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Significant  progress  has  been  achieved  under  grant  AFOSR-86-0031  in  understand¬ 
ing  the  physics  underlying  ultrafast  transient  phenomena  that  occurs  in  the  sem¬ 
iconductor  microstructures.  This  understanding  is  essential  for  making  the  neces¬ 
sary  advances  for  the  future  generation  of  ultrafast  microelectronic  devices. 

Microstructure  samples  have  been  obtained  primarily  from  H.  Morkoc  of  University 
of  Illinois  and  other  groups  such  as  M.  Niigaki  of  Hamamatsu  Photonic  KK,rR(dl32s 
Bertaska  of  McDonnei  Douglas,  and  Emil  Koteles  of  GTE. 

The  following  summarizes  our  accomplishments  in  the  major  effort  on  III- V  sem¬ 
iconductor  microstructures  and  a  minor  effort  on  II-VI  compound. 

1.  Nonequiiibrium  Phonon  Effects  on  the  Energy  Relaxation  and  Lifetime 
of  Photogenerated  Carriers  in  GaAs  MQW 

We  have  carried  out  time-resolved  photoluminescence  with  the  finest  time- 
resolution  of  2ps  performed  to  date  to  investigate  the  ultrafast  processes  which 
occurs  in  undoped  GaAs/AlGaAs  microstructures. 

The  existence  of  a  large  population  of  nonequilibrium  (hot)  phonons  after  an  initial 

o 

rapid  carrier  cooling  in  an  undoped  multiple  GaAs  quantum  well  (55A)  structure 
excited  by  500fs  laser  pulses  is  experimentally  verified  by  time-resolved  photo- 
luminescence  studies  using  a  2ps-time-rcsolution  streak  camera.  The  hot  phonon 
population  is  directly  monitored  by  measuring  time-resolved  phonon  replica 
luminescence  below  n=l  electron  to  heavy-hole  transition  energy.  It  reaches  to  a 
maximum  in  ~3ps  after  the  end  of  laser  pulse  and  decays  with  a  JOps  time  con¬ 
stant.  Its  effects  on  carrier  dynamics  are  the  following  two  aspects.  First,  the 
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energy  relaxation  of  carriers  after  an  initial  rapid  cooling  is  substantially 
suppressed  due  to  re-absorption  of  hot  phonons  by  carriers.  Second,  the  carrier 
density  decay  is  strongly  influenced  by  the  hot  phonons  in  the  system  which 
stimulate  phonon  replica  emission.  This  stimulated  phonon  replica  emission  is 


mainly  radiated  and  highly  polarized  in  the  well  plane  (x-y  plane).  An  effective 
carrier  depletion  time  [density  decreases  by  a  factor  of  c-1  from  ne  it  =0)]  is 
obtained  to  be  as  short  as  lOps. 

This  work  is  of  particular  importance  to  the  semiconductor  community  because  it 
solves  the  long  standing  controversial  issue  regarding  the  phonon  assisled  transi¬ 


tions  below  n=l  electron  to  hole  transition  in  GaAs  quantum  wells  that  was  ini¬ 
tially  suggested  by  Professor  Holonyak’s  group  (see  for  example,  Phys.  Today,  40, 
Oct.  1980).  We  have  explicitly  demonstrated  that  the  hot  phonons  emitted  by  hot 
electrons  can  stimulate  phonon-replica-emission. 

Part  of  this  work  was  accepted  to  be  published  in  Phys.  Rev.  B. 

2.  Observation  of  Ultrafast  Lateral  Diffusion  of  Photogenerated  Carriers  in 
An  Asymmetric  GaAs  Quantum  Well 


i 

I 

i 

e 


An  experimental  investigation  of  in-plane  expansion  of  photogene  rated  carriers  in  a 
single  asymmetric  GaAs  quantum  well  is  carried  out  using  a  streak  camera  and  a 
subpicosecond  laser  source  at  4.3K.  It  is  found  that  the  carrier  expansion  process  is 
enhanced  in  the  well  plane  by  restricting  carrier  expansion  in  the  well  direction. 
The  effective  diffusivity  D'  of  photoexcited  carrier  in  the  asymmetric  well  has 
been  directly  determined  to  be  l()6c/n 2 1  s  at  4.3K,  which  is  about  four  orders 
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magnitude  larger  than  the  thermodynamic  value  in  bulk  GaAs.  The  expansion 
velocity  is  determined  to  be  3.2xl08cm  /  s  ,  which  is  four  times  larger  than  the 
initial  Fermi-velocity  of  the  electrons. 

The  observation  of  an  anomalously  fast  carrier  expansion  process  in  an  asymmetric 
quantum  well,  but  not  in  the  symmetric  well  is  most  important.  It  provides  a 
new  insight  to  the  problem  of  how  to  design  semiconductor  structures  in 
order  to  obtain  the  desired  carrier  transport. 

This  work  was  submitted  to  Phvs.  Rev.  Lett,  for  publication. 

3.  Determination  of  Band  Offsets  via  Optical  transitions  in  Ultrathin 
Quantum  Wells 

This  work  demonstrates  what  experimental  data  from  optical  transition  should  be 
used  and  how  it  is  sensitive  to  the  Q  -value  of  band  offset.  We  have  shown  that 
the  precise  Q  -value  can  be  determined  by  systematically  measuring  the  n=l  light- 
and  heavy-hole  energy  spacing  (A£  )  in  the  sensitive  zone  of  well  width  (15  to 

o 

80A).  Furthermore,  it  is  also  shown  that  different  connection  rules  and  hole 
masses  yield  different  Q  -value  and  the  appropriate  connection  rule  can  be  discrim¬ 
inated  also  by  systematically  measuring  A£  as  well  as  the  energy  spacing  between 
ground  and  first  electron  eigen  state  (AE  J2  )  in  the  sensitive  zone  of  L. . 

This  work  was  published  in  Phys.  Rev.  B,  37,  1048  (1988). 

Recently,  we  are  investigating  the  influence  of  different  connection  rules  on  the 
exciton  binding  energy  in  1II-V  semiconductor  quantum  wells,  and  hence  on  the 


determination  of  band  offset  based  on  the  data  of  A E . 


4.  Physics  in  semiconductor  GaAs  and  GaSe  under  laser-driven  shock-wave 
compression 
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Using  a  pump-probe  technique,  the  photoluminescence  peak  of  bulk  GaAs  was 
observed  to  blue  shift  and  split  into  two  components  under  the  laser-induced  shock 
wave  loading  condition.  The  two  components  correspond  to  the  transitions  from 
the  T6  conduction  band  to  the  valence  heavy-  and  light -hole  subbands,  because  of 
symmetry  breaking  by  the  uniaxial  shock  compression  along  the  [001]  direction. 
From  the  blue  shift,  we  deduced  our  picosecond- laser-induced  shock  pressure  of  ~ 
10k  bar. 


In  GaSe,  we  observed  a  broadening  of  the  spontaneous  emission  which  is  propor¬ 
tional  to  the  shock  pressure  and  attributed  to  a  shock-wave-induced  exciton  colli¬ 
sion  mechanism  due  to  the  directional  motion  of  particles  in  the  shocked  region. 
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These  works  were  published  in  Appl.  Phys.  Lett.  51,  1780  (1987)  and  Appl.  Phys. 
Lett  52,  93  (1987). 


5.  Optical  transition  and  recombination  lifetime  in  quasi-zero  dimensional 
electron  system  in  CdSx  Sc  X_x 
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We  report  on  the  observation  of  optical  transitions  between  quantized  levels 
(IS  ,\P  )  in  the  conduction  and  valence  bands  in  quasi-zero  dimensional  electron 
systems  in  CdSx  Sc  !_x  by  4K  and  300K  photoluminescence  measurements.  From 
the  relative  luminescence  intensities  arising  from  the  IS  —IS  and  IP  — 1 P 


transitions,  it  is  shown  that  the  photoexcited  electrons  (holes)  occupied  in  the 
different  quantized  levels  in  conduction  (valence)  band  are  not  coupled,  or  only 
weakly  coupled.  Picosecond  luminescence  studies  of  the  optical  transitions  using  a 
streak  camera  reveal  a  significant  influence  of  three-dimensional  confinement  on  the 
recombination  probabilities  of  photogenerated  electrons  and  holes.  The  recombina¬ 
tion  life  time  of  electrons  and  holes  in  quasi-zero  dimensional  system  decrease  with 
the  diameter  (d)  of  system  from  210ps  for  d=10.2nm  to  70ps  for  d=7.4nm  at  4K. 
The  ratio  of  recombination  life  time  for  IS  —IS  transition  and  IP  —IP  transition 
seems  to  be  independent  of  the  diameter  of  system,  and  is  measured  to  be  about 
3.5.  This  ratio  can  be  interpreted  by  the  calculated  squared  transition  matrix  ele¬ 
ments  and  the  nonradiative  processes. 

This  work  was  published  in  Appl.  Phys.  Lett.  51,  1839  (1987). 
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Effects  of  nonequilibrium  phonons  on  the  energy 
relaxation  and  recombination  lifetime  of  photogenerated 
carriers  in  undoped  GaAs  quantum  wells 

Kai  Shum,  M.  R.  J unnarkar ,  H.  S.  Chao,  and  R.  R.  Alfano 

Institute  for  L'ltrafast  Spectroscopy  and  Lasers 
Departments  of  Electrical  Engineering  and  Physics 
The  City  College  of  New  York,  New  York,  New  York,  10031 

and 

H.  Morkoc 

Coordinated  Science  Laboratory,  University  of  Illinois  at 
Urbana-Champaign,  Urbana,  Illinois  61801 

ABSTRACT 

Time-resolved  and  time- integrated  photoluminescence  studies  of 
an  undoped  multiple  GaAs  quantum  well  structure  excited  by  0.5ps 
laser  pulses  have  revealed  several  important  experimental  observations 
on  the  behavior  of  photoexcited  carriers.  A  large  population  of  none¬ 
quilibrium  longitudinal  optical  phonons  produced  in  the  energy  relaxa¬ 
tion  process  of  hot  carriers  manifests  itself  by  the  nonequilibrium  pho¬ 
non  stimulated  phonon-replica  which  located  at  ~30meV  below  n=l 
electron-hole  transition.  The  energy  relaxation  is  substantially 
suppressed  due  to  the  existence  of  nonequilibrium  phonons  after  an  ini¬ 
tial  rapid  cooling  (0-5ps).  The  number  of  photoexcited  carriers 


decreases  anomalously  fast  within  the  first  30ps  after  the  excitation  by 
laser  pulse.  An  effective  carrier  depletion  time  is  determined  to  be  as 
short  as  lOps.  A  mechanism  which  leads  to  such  a  short  carrier  deple¬ 
tion  time  is  associated  with  the  nonequilibrium  phonon  stimulated 
phonon- replica. 


PACS  numbers:  63.20.Kr,  72.20.Jv,  73.40.Lq. 
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I  Introduction 


Recently,  there  have  been  growing  interest  in  studying  the  dynamics  of  photo¬ 


generated  carriers  in  semiconductor  quantum  well  structures  since  it  reflects  the  fun¬ 


damental  interactions  among  electrons,  holes  and  phonons.  These  interactions  deter¬ 


mine  the  performance  of  ultrahigh  speed  electronic  and  optoelectronic  devices.  In  the 


case  of  bulk  semiconductors,  the  slowing  of  hot  carrier  relaxation  is  attributed  to 


either  screening1  of  electron-phonon  interaction  or  re-heating  of  carriers  by  nonequili¬ 
brium  (NE)  phonons  produced  in  the  relaxation  process  of  hot  carriers.”’3  A  compara¬ 


tive  study4  of  the  hot  electron  cooling  rates  in  undoped  multiple  quantum  well 


(MQW)  structures  and  in  bulk  GaAs  by  time-resolved  measurements  of  optical 


absorption  and  gain  concluded  that  the  rates  were  approximately  same  at  a  photogen¬ 
erated  carrier  density  of  2.5xl017  cm"3.  This  is  expected  from  a  simple  theory5  in 


which  lattice  is  treated  as  a  heat  bath  for  quasi-equilibrium  carriers.  In  the  case  of 


modulation  doped  MQW  structures,  quasi-steady-state  experiments  have  been  carried 
out  from  two  groups6,7  to  investigate  the  interaction  of  electrons  and  holes  with  pho¬ 
nons.  These  experiments  have  generated  conflicting  results  regarding  the  presence6  or 
absence7  of  NE  phonons.  Time-resolved  PL  measurements  with  ~20ps  time  resolution 


were  reported  by  Ryan  et  al .  who  found  that  the  cooling  of  hot  carriers  was 


anomalously  slow  after  40ps. 


Reviewing  the  aforementioned  previous  works  two  important  points  were 


noticed.  First,  one  can  not  extract  information  on  the  initial  carrier  relaxation  process 


either  from  quasi-steady-state  experiments  or  from  time-resolved  measurements  with 


time  resolution  greater  than  lOps.  To  assess  the  importance  of  NE  phonon  effect  on  the 


energy  relaxation  of  hot  carriers  one  uses  a  corresponding  theory  to  match  the  experi¬ 


mental  carrier  temperature  cooling  curve  or  the  logarithm  of  power  loss  as  function 


of  the  inverse  of  the  carrier  temperature.  There  has  been  no  report  on  the  direct  evi¬ 


dence  for  existence  of  NE  phonons  in  quantum  well  structures  under  high  photoexci- 
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tation.  It  is  necessary  to  study  the  initial  carrier  relaxation  process  to  substantiate  the 
existence  of  NE  phonons  and  its  importance  on  the  hot-carrier  relaxation  since  the 
phonon  lifetime  is  about  7ps.  Second,  one  can  not  easily  obtain  information  on  photo¬ 
generated  carrier  lifetime  in  modulation  doped  MQW  structures  because  the  large 
doped-in  carrier  density  is  comparable  to  the  photoexcited  carrier  density.  Further¬ 
more,  carrier  lifetime  can  influence  the  cooling  of  hot  carriers.9  Therefore,  it  is  neces¬ 
sary  to  study  photogenerated  carrier  dynamics  in  undoped  MQW  structures  in  order 
to  obtain  information  both  on  the  energy  loss  rates  and  on  the  carrier  lifetime. 

In  this  paper,  experimental  observations  are  reported  from  the  measurements  of 
time-resolved  photoluminescence  (PL)  with  a  2ps  time  resolution  and  time-integrated 
luminescence  spectra  from  an  undoped  GaAs  MQW  structure  with  dependences  on 
lattice  temperature,  excitation  intensity,  and  polarization.  The  NE  phonons  emitted 
by  hot  electrons  are  directly  observed  by  measuring  the  time-integrated  as  well  as 
time-resolved  phonon  replica  luminescence  lying  below  the  n=l  electron  to  heavy- 
hole  transition  energy.  The  energy  relaxation  of  hot  electrons  is  found  to  be  substan¬ 
tially  suppressed  when  a  large  population  of  NE  phonons  is  present  after  an  initial 
rapid  cooling.  The  photoexcited  carrier  density  extracted  from  a  fitting  of  time- 
resolved  PL  profiles  at  different  emitted  photon  energies  decreases  non-exponential ly 
and  very  rapidly  within  the  first  30ps.  An  effective  carrier  depletion  time  is  deter¬ 
mined  to  be  as  short  as  lOps.  The  mechanism  leading  to  such  short  carrier  depletion 
time  is  associated  with  NE  phonon  stimulated  phonon-replica. 

11  Sample 
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The  undoped  GaAs/Al^Ga^As  MQW  structure  investigated  was  grown  by 
molecular-beam  epitaxy  on  a  (001)-oriented  undoped  GaAs  substrate.  The  MQW  con- 

©  o 

sists  of  50  periods  of  55A  thick  GaAs  and  100A  thick  Al03GaQ7As  layers,  and  fol¬ 
lowed  by  a  1.2 um  GaAs  buffer  layer.  The  lateral  size  of  the  structure  is  about 
2x4mm .  The  good  quality  of  sample  is  confirmed  by  well  resolved  heavy-hole  and 
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light-hole  excitonic  structures  of  room  temperature  PL  spectra.  The  FWHM  for  the 
heavy-hole  excitons  is  about  13meV.  The  sample  was  mounted  on  the  cold  finger  in 
an  optical  helium  cryostat. 

Ill  Experiment 

An  ultrashort  laser  pulse  of  0.5ps  duration  at  620nm  was  used  to  excite  the 
electron-hole  pairs  in  the  sample.  Most  of  light  is  absorbed  in  both  the  GaAs  wells 
and  the  AlGaAs  barriers.  Photogenerated  carriers  in  the  barriers  will  either  diffuse  to 
the  GaAs  buffer  layer  or  be  captured  into  the  wells.  The  latter  will  result  in  uncer¬ 
tainty  about  carrier  density  in  the  wells  by  a  fact  of  ~J.  Residual  transmitted  light 
is  absorbed  in  the  GaAs  buffer  layer.  The  diameter  of  excitation  area  is  about  160wm . 
To  avoid  damaging  of  the  sample  surface,  the  maximum  excitation  power  is  adjusted 
such  that  photogenerated  carrier  density  is  on  order  of  I019cm'3.  The  0.5ps  laser 
pulse  was  generated  from  a  colliding-pulse  passive  mode-locked  dye  laser  and 
amplified  by  a  four-stage  dye  amplifier  pumped  with  a  frequency  doubled  Nd:  YAG 
laser  operated  at  20Hz.  In  order  to  keep  the  time  resolution  of  Hamamatsu  streak 
camera  system  within  2ps,  PL  was  spectrally  resolved  using  different  narrow  band 
filters.  The  luminescence  intensities  were  corrected  for  the  nonlinearity  of  streak 
rate,  the  spectral  response,  and  the  transmission  of  each  narrow  band.  For  the  time- 
integrated  PL  measurements,  a  Stanford  boxcar  and  a  GaAs  photomultiplier  were 


h 


used. 

IV  Results  and  Discussion 
1.  Steady-state 
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Time-integrated  luminescence  spectra  from  the  MQW  structure  excited  by  the 
0.5ps  laser  pulse  at  various  lattice  temperatures  (7 L  )  are  shown  in  Fig.l.  These  spec¬ 
tra  were  taken  in  a  conventional  backward  Raman  configuration  r  (y  ,y  )z  ,  where  z  is 
the  growth  direction.  In  order  to  eliminate  radiation  from  the  edges  of  the  sample. 
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the  luminescence  spot  was  first  imaged  on  an  aperture  with  aid  of  a  streak  camera 
and  then  re-focused  onto  a  vertical  slit  of.a  grating  spectrometer.  Several  features  are 
displayed  in  the  data  shown  in  Fig.l.  (1 )  The  emission  band  A  on  high  energy  side  of 
the  spectra  arises  from  the  recombination  of  photogenerated  n=l  electron-hole.  The 
peak  of  A  shifts  towards  the  low  energy  side  as  the  TL  increases.  An  high  energy 
tail  on  the  peak  A  develops  with  increasing  TL  . 

(2)  A  broad  emission  band  C  at  low  energy  side  of  the  spectra  arises  from  the 
GaAs  buffer  layer.  The  total  emission  intensity  of  C  decreases  as  TL  increases. 

(3)  The  most  interesting  feature  of  the  spectral  data  is  the  appearance  of  an 
emission  band  B,  located  about  30meV  below  the  n=l  electron-heavy-hole  transition. 
This  B  emission  band  is  attributed  to  the  NE  phonon  stimulated  phonon-replica.  The 
following  four  reasons  support  this  assignment  to  the  B  band: 

(0  This  B  emission  band  does  not  appear  in  the  PL  spectra  taken  at  low  power 
excitation  about  lW/cm2  at  4.3K  using  a  cw  488nm  line  of  an  argon-ion  laser.  How¬ 
ever,  the  PL  spectra  using  a  very  weak  train  of  laser  pulses  (l20fs)  directly  from  the 
colliding-pulse  passive  mode-locked  dye  laser  as  the  excitation  source,  with  an  excita¬ 
tion  power  density  in  range  of  10-5— 10_3W  /  cm2,  shows  a  weak  electron  to  acceptor 
emission  band  separated  by  17meV  from  the  n=l  electron-heavy-hole  transition 
accompanied  by  its  phonon-replica10  at  low  temperatures.  This  extrinsic  emission 
band  disappeared  completely  when  the  sample  temperature  was  raised  to  ~  80K, 
whereas  the  B  band  displayed  in  Fig.l  exists  up  to  room  temperature.  Moreover,  the 
relative  time-integrated  PL  intensity  of  the  B  to  A  band  under  0.5ps  light  pulse  exci¬ 
tation  shown  in  the  Fig.l  decreases  as  the  excitation  power  density  decreases.  The  B 
band  in  Fig.l  is  hardly  visible  when  the  excitation  power  density  is  lower  than 
lO-3/**, ,  where  Pm  —1012W  /  cm2  is  the  maximum  value  of  the  excitation  power  den¬ 
sity.  This  contradicts  what  is  generally  expected  for  impurity  emission.  Since  the 
concentration  of  acceptors  is  low  in  our  sample  as  confirmed  by  the  luminescence  stu- 
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dies  using  the  weak  train  of  120fs  light  pulses,  an  impurity  emission  band  should  be 
more  readily  apparent  and  pronounced  at  lower  excitation.  Therefore,  the  B  emission 
band  shown  in  Fig.  1  can  not  be  attributed  to  electron  to  acceptor  luminescence  under 
high  excitation. 

(ii)  The  intensity  ratio  between  the  B  and  A  emission  bands  as  indicated  by  tri¬ 
angle  in  Fig.2  increases  as  lattice  temperature  increases.  To  convince  the  emission 
band  B  arises  from  NE  phonon  stimulated  recombination  process,  we  calculate  the 
sum  of  occupation  numbers  both  for  the  equilibrium  phonons  (  lattice  temperature 
T i  )  and  for  the  NE  phonons  (  carrier  temperature  Tc  )  by  assuming  the  carrier  tem¬ 
perature  is  same  as  the  effective  temperature  for  the  NE  phonons  after  30ps.  This 
will  be  discussed  later  with  the  carrier  temperature.  The  calculated  result  is  shown 
by  the  solid  line  in  Fig.2.  The  fitting  of  the  total  phonon  occupation  number  to  the 
obtained  intensity  ratio  from  the  spectra  is  impressively  good.  This  implies  that  the 
intensity  of  band  B  is  well  correlated  with  the  LO-phonon  population.  Because  of  the 
participation  of  NE  phonon  the  intensity  of  B  band  (see  also  time-resolved  lumines¬ 
cence  at  780nm  in  Fig.4)  is  strongly  enhanced  especially  at  low  temperatures.  This 
further  supports  our  assignment  of  the  B  band  to  the  NE  phonon  stimulated  phonon 
replica. 

(iii)  It  is  also  not  possible  to  attribute  the  emission  band  B  to  the  emission  from 
renormalized  band-band  transition  because  the  peak  position  of  the  A  band  (766.7nm) 
and  the  B  band  at  4.3K  do  not  change  with  variation  of  excitation  intensity  by  a  fac¬ 
tor  of  ~  103.  A  further  support  to  above  statement  is  the  fact  that  the  spectral  posi¬ 
tion  of  peak  A  exactly  coincides  with  each  other  for  the  two  different  luminescence 
studies  using  different  light  excitation  sources:  one  is  the  weak  120fs-pulse  train  with 
a  repetition  rate  of  125MHz,  and  the  other  is  the  amplified  0.5ps-pulse  with  a  20Hz 
repetition  rate. 

(iv)  We  have  also  studied  the  polarization  of  the  emission  bands  A,  B.  and  C  by 


J£JiV>V|Vj 


"mcFonCTiiTW!  jutkj*  nr  urainwmflrjv  wv  jti  *mm  **  u*  - 


-  8  - 

measuring  "right-angle"  11,1“  z  (y  ,  x  )y  and  z  (y  ,  z  )y  integrated  luminescence  spectra. 
It  was  found  that  the  B  band  emitted  at  the  sample  edge  along  the  y-direction  (see 
inset  of  Fig.3)  was  highly  polarized  in  the  x-direction.  The  A  and  C  emission  bands 
were  depolarized,  independent  of  excitation  power  density  and  lattice  temperature. 
The  intensity  ratio  between  the  x-  and  z-direction  for  the  B  band  is  about  20  in  the 
lattice  temperature  range  from  4.3  to  300K  at  full  excitation  power  density  P„, .  But, 
the  ratio  is  strongly  dependent  on  the  excitation  power  density  at  a  given  lattice  tem¬ 
perature.  For  example,  at  TL  =  20A' ,  the  ratio  decreases  from  20  to  —6  as  the  excita¬ 
tion  power  density  decreases  from  P,„  to  0.017/5,,.-  Fig.3  shows  the  spectra  for 
T i  =  100 A'  at  two  different  excitation  intensities.  The  lower  solid  trace  was  taken 
in  z  (y  ,  x  )y  configuration  (TE  polarization)  at  the  excitation  density  of  0.0l7/’m  .  It  is 
identical  to  the  spectrum  taken  in  z  (y  ,  z  )y  configuration  (TM  polarization)  which  is 
not  shown.  The  B  band  disappears  due  to  the  absence  of  large  population  of  NH  pho¬ 
nons  at  the  low  excitation  and  high  lattice  temperature.  The  upper  solid  and  broken 
traces  are  the  TM  and  TE  spectra  at  the  full  excitation  Pn, ,  respectively.  The  inten¬ 
sity  ratio  of  TE  and  TM  remains  20.  excitation  Pm .  It  should  be  pointed  out  that  the 
energy  position  of  A  band  in  Fig.3  is  lower  than  that  in  Fig.l  by  ~4meV  while  the 
position  of  B  band  remains  unchanged.  This  low  energy  shift  of  A  band  is  due  to  the 
self-absorption  effect.  The  polarization  behavior  of  B  band  is  consistent  with  the 
results  from  Raman  scattering  experiments  reported  by  Zucker  et  al.]~  where  the 
pump  photons  are  provided  by  the  external  laser  source.  In  our  case  the  pump  pho¬ 
tons  arise  from  the  recombination  of  n=l  electron  and  n=l  heavy-holc  at  the  subband 
edges.  This  process  is  strongly  enhanced  by  the  presence  of  a  large  number  of  NE 
phonons  emitted  by  hot-electrons  at  high  excitation.  The  depolarized  feature  of  the 
band  C  is  expected  since  the  emission  is  from  the  GaAs  buffer  layer.  However,  the 
polarization  behavior  of  the  emission  band  A  is  not  understood.  In  terms  of  the  selec¬ 
tion  rule  of  dipole  recombination  described  by  Iwamnra  ct  a/.13,  the  emission  intensity 


of  n=l  electron- light-hole  for  the  TE  polarization  should  be  four  times  larger  than 
that  for  the  TM  polarization.  There  should  be  no  TM  emission  arising  from  n=l 
electron-heavy-hole  recombination.  It  should  be  mentioned  here  that  the  B  band  in 
the  Fig.l  which  was  detected  along  the  z-direction  is  depolarized  and  much  weaker, 
at  least  a  fact  of  20,  than  detected  along  the  y-direction.  The  former  is  expected  due 
to  the  symmetry  and  the  latter  provides  an  important  rule  for  understanding  a  fast 
carrier  density  decreasing  process  which  will  be  discussed  in  the  timc-rcsolved  section. 

In  order  to  further  substantiate  our  assignment  of  the  B  band,  two  additional 
questions  must  be  addressed: 

Cl)  Why  no  phonon  replica  appears  on  the  high  energy  side  of  the  A  band,  if  a 
large  NE  phonon  population  really  arises?  A  high  energy  replica  would  arise  from  a 
recombination  of  an  energetically  elevated  electron  (hole)  by  absorbing  a  phonon  and 
hole  (electron).  However,  the  elevated  electron  (hole)  will  be  very  quickly  scattered 
by  the  other  hot  electrons  (holes)  through  a  strong  electron-electron,  electron-hole, 
and  hole-hole  interaction  before  it  recombines  with  a  hole  (an  electron).  Therefore, 
the  high  energy  replica  can  not  be  observed  in  PL  spectra.  But.  this  re-heating  process 
should  result  in  a  retardation  of  hot  carrier  cooling  which  is  indeed  consistent  with 
the  time-resolved  PL  data  described  in  the  next  section. 


(2)  Why  is  the  energy  separation  between  the  peaks  of  A  and  B  bands  (/; Ag  ),  on 
order  of  ~  30meV  which  is  smaller  than  expected  (3GmeV)7  The  30m  eV  energy 
separation  in  fact  gives  a  further  strong  support  of  our  assignment  of  the  B  band.  In 
an  equilibrium  state,  the  interaction  strength  of  electron  and  LO  phonon  in  a  bulk 


GaAs  is  proportional  to  -iy,  where  Q  is  the  phonon  wave  vector.  Therefore,  the  pho¬ 


nons  in  vicinity  of  Q  =  0  would  be  expected  to  more  strongly  couple  with  electrons 
giving  rise  to  a  phonon  replica  at  just  36meV  below  its  primary  emission  band.  In 
our  highly  photoexcited  MQW  sample,  LO  phonons  are  driven  to  a  NE  state  by  a  rapid 
initial  energy  relaxation  process  of  photoexcited  hot  electrons.  The  develop  of  these 


NE  both  in  time  domain  and  in  two-dimensional  wave  vector  space  (</*,  parallel  to  the 
well  plane)  has  been  recently  studied  by  several  groups.14,15  Two  important  results 
on  the  number  of  NE  LO  phonon  population  are  :  (i)  At  given  q ,  it  can  reach  to  a 
maximum  in  about  0.5  -  3ps  and  then  decrease,  (ii)  At  a  given  time,  it  increases  stee¬ 
ply  from  a  minimum  of  q  mm,  reaches  to  a  maximum  at  qc ,  and  then  decreases.  What 
is  essential  concerning  the  question  is  that  there  is  a  maximum  NE  LO  pho¬ 
non  population  which  is  located  both  in  a  narrow  wave  vector  space  of  the 
vicinity  of  qc  and  in  the  well  plane  and  lasts  in  a  short  time  period.  This 
portion  of  NE  phonons  behaves  like  coherent  Bosons.  The  existence  of  the  "coherent 
Bosons"  will  tend  to  increase  the  number  of  the  Bosons  (near  qc  )  in  the  system  with 
a  rate  proportional  to  the  present  number  of  Bosons  resulting  in  a  stimulated  phonon 
emission  process.  This  process  can  give  rise  to  a  phonon  replica  emissions,  when  an 
electron  in  the  conduction  subband  recombines  with  a  hole  in  the  valence  subband, 
not  only  a  photon  but  also  a  phonon  or  more  phonons  will  be  emitted  to  join  the 
Boson  system.  In  order  to  conserve  the  momentum  of  electron-phonon  system  in  this 

recombination,  the  electron  must  possess  a  momentum  of  qr  (neglecting  the  momen- 

2  2 
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turn  of  the  emitted  photon)  which  is  associated  with  a  kinetic  energy  of  E„  =  — - 

2  mr 

Therefore,  the  energy  EAB  should  be  equal  to  LO  phonon  energy  minus  E  which  is 
about  6meV  corresponding  to  the  values  of  1.03xl0bem_l  and  0.067m 0  for  qt  and  elec¬ 
tron  effective  mass  m, .  Thus,  this  gives  30meV  for  EAB  . 

Over  years,  there  has  been  much  debate13,16"18  over  the  interpretation  of  the 
spectral  features  locate  below  the  n-1  electron  to  heavy-hole  transition  energy, 
mainly  because  many  species  exist  in  this  energy  region:  LO-phonon  replica,  impurity 
states,  and  many  body  band  gap  renormalization.  Holonyak  and  co-workers17  have 
demonstrated  the  LO-phonon  participation  in  GaAs-AIGaAs  based  QW  laser  emission 
by  observing  mre  than  one  LO-phonon  sidebands  in  the  laser  operation  below  the 
n-1  confined-particle  transition.  Recently,  Skolnick  ct  al ,18  confirmed  that  a  peak 
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below  the  n=l  electron- heavy-hole  transition  is  phonon  replica  by  observing  the  cou¬ 
pling  at  both  the  GaAs-  and  InAs-like  bulk  LO-phonon  energies  of  InGaAs  QW.  Our 
measurements  support  these  arguments. 

2.  Time-resolved 

Time-resolved  PL  intensities  of  the  MQW  structure  at  4.dK  were  measured  over 
the  spectral  range  from  720  lo  780nm.  Lour  representative  streaks  at  given  selected 
wavelengths  are  shown  in  Fig.4.  The  emission  centered  at  770nm  arises  from  the 
recombination  of  band  edge  electron  and  heavy-hole.  The  radiation  with  wavelengths 
less  than  770nm  are  from  the  recombination  of  hot  carriers.  The  time-resolved 
luminescence  intensity  of  the  phonon  replica  is  centered  at  780nm.  Each  lumines¬ 
cence  profile  is  an  average  from  10  individual  shots  using  a  prepulse19  for  averaging. 
The  left  peak  of  the  dotted  curve  is  the  prepulse  which  reflects  the  2ps  temporal  reso¬ 
lution.  The  right  peak  on  the  same  curve  is  the  Rayleigh  scattering  light  from  the 
sample  surface  which  defines  the  "zero"  time  for  our  analysis.  Several  features 
appear  in  the  data  displayed  in  the  Fig.4.  (i)  The  rise  time  for  all  the  luminescence 
profiles  are  instrumental;  a  2ps  up-limit  reflects  the  rapid  thermalization,  capturing, 
and  initial  cooling,  (ii)  The  shape  of  rise-part  of  luminescence  profile  at  780nm  is 
similar  to  that  at  770nm,  but  delayed  by  '"'dps.  This  implies  that  the  emission  at 
780nm  does  not  originate  from  the  same  band  as  the  emission  at  770nm;  otherwise, 
the  rise-part  of  the  luminescence  at  780  should  start  at  the  same  zero  point  as  that  at 
770nm.  The  delay  of  '''dps  is  consistent  with  our  assignment  of  the  emission  band 
center  at  780nm  to  the  stimulated  phonon-replica.  This  dps  is  the  time  required  to 
establish  a  large  population  of  NE  phonons  at  qr .  (iii)  The  luminescence  decay  time20 
of  dOps  (broken  curve)  for  the  phonon  replica  emission  at  780nm  is  shorter  than  the 
luminescence  decay  time  of  60ps  at  770nm.  This  supports  that  the  emission  at 
780nm  can  not  be  due  to  impurity  emission  at  high  excitation. 
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To  study  the  NE  phonon  effect  on  energy  relaxation  process  quantitatively,  car¬ 
rier  temperature  and  density  as  function  of  time  should  be  determined  simultane¬ 
ously.  An  expression  in  the  time  domain  is  introduced  to  fit  the  time-resolved  PL 
data  by  using  two  adjustable  parameters,  namely,  the  carrier  density21  (nc )  and  the 
electron  temperature  (Tc ).  For  the  direct  optical  transition  in  MQW,  the  lumines¬ 
cence  intensity  is  given  by: 

I(E,,t  )  =  C,(\-e~"T')[\M,-hh  I  2p,  f  ,  Phh  f  h  +  I  Me  -ih  Ml) 

where  the  prjlhji,  ~~mcMhJh  are  the  density  of  state  for  electron,  heavy-hole,  and 
light-hole;  the  A/r  (ih  >  is  the  matrix  element22  for  electron  to  hhUh)  transition;  C, 
absorbed  all  the  constant  factors  including  the  corrections  for  detector  response  and 
the  transmission  of  each  narrow  band  filters  used;  r,  is  the  rise  time  of  luminescence 
and  is  set  to  be  ~  lps  which  is  an  up-limit  of  thermal izalion  time  of  electron-hole 
system; 


is  the  Fermi-Dirac  distributions  for  the  electrons  in  conduction  band  [with  subscript 
e)  and  for  the  holes  in  valence  band  (with  subscript  h). 

An  unique  set  parameters  of  Tc(t  )  and  u,  (r  )  was  determined  by  consistently 
fitting  all  the  luminescence  profiles  detected  at  different  photon  energies  E, .  Three 
calculated  luminescence  profiles  corresponding  the  experimental  data  are  shown  in 
Fig.4  by  the  thick  solid  curves. 

The  determined  Tc  as  function  of  time  is  plotted  as  the  solid  curve  in  Fig.5.  The 
shaded  area  reflects  the  extent  of  the  uncertainty  in  deducing  the  carrier  temperature 
within  the  first  4ps  due  to  our  limited  time  resolution.  The  plotted  data  in  Fig.5  is 
interpreted  as  follows.  Carrier-carrier  collisions  quickly  leads  to  a  thcrmalizcd  distri¬ 
bution  at  very  high  Tc  within  lps  due  to  a  large  number  of  hot  carriers  excited.23 
Although  an  initial  thermal  ization  process  of  carriers  can  only  be  probed  by  fern- 
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tosecond  spectroscopy.  An  initial  cooling  of  the  thermal  i/.ed  distribution  is  stu¬ 
died  with  our  present  time  resolution  of  2ps,  providing  information  about 
the  energy  relaxation  and  the  density  decay  of  photogenerated  carriers.  The 
initial  cooling  within  first  5ps  is  very  fast  (250K/ps)  reflecting  both  the  screening1 
and  NE  phonon  effects  are  small.  After  the  initial  cooling,  a  large  number  of  longitu¬ 
dinal  optical  phonons  in  a  finite  wavevector  space  accumulates  due  lo  the  finite  pho¬ 
non  lifetime.24  The  phonons  emitted  by  the  hot  electrons  are  reabsorbed  by  the  elec¬ 
trons  as  a  reverse  process  of  the  emission  giving  rise  to  a  slower  cooling  for  the  hot 
carriers.'5  It  should  be  emphasized  here  that  the  rapid  initial  cooling  rules  out  the 
importance  of  significant  screening  of  electron-phonon  interaction  in  the  present 
study.  If  the  initial  screening  was  important,1  the  initial  rapid  cooling  would  be 
slowed  and  a  NE  phonon  population  would  not  be  built  up. 

Another  important  observation  displayed  in  the  Eig.5  is  that  the  time  constant 
for  the  slow  decay  component  of  carrier  temperature  cooling  curve  is  ~ 30 ps  which  is 
same  as  the  decay  time  of  the  NE  phonon  stimulated  phonon  replica  emission  at 
780nm  (see  Eig.4).  This  suggests  that  ~5ps  is  required  for  the  electron  sys’em  and 
the  NE  phonon  system  to  be  essentially  equilibrated  with  each  other.  The  coupled 
electron-phonon  system  decays  with  a  common  decay  constant  of  about  30ps.  It 
should  be  emphasized  that  this  time  constant  (30ps)  is  the  lifetime  of  NE  phonons 
which  differ  from  the  equilibrium  phonon  lifetime  of  ~7ps.  The  longer  lifetime  for 
the  NE  phonons  than  equilibrium  phonon  is  due  to  the  coupling  between  the  hot  elec¬ 
trons  and  NE  phonons. 

The  determined  quasi-Fermi-encrgies  for  electrons  ar  '  holes  are  plotted  in  Fig.6. 
The  changes  of  u,  and  uh  are  very  rapid  with  the  first  lOps  reflecting  a  rapid  decrease 
of  the  carrier  density.  The  behavior  of  the  degeneracy  for  electrons  and  holes  with 
function  of  time  is  reversed  due  the  difference  in  effective  masses  for  the  electrons 


and  holes. 
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The  quasi-equilibrium  distribution  functions  for  the  electrons  and  the  holes  can 
be  experimentally  determined  from  Tc(t  )  and  ).  Using  these  distributions  the 

energy  loss  rates  for  the  electrons  and  holes  can  be  obtained.  The  energy  loss  rate 
( Pcm  )  is  defined  as  follows: 


P c  ,  hh  — 


d  <E  > 


e  Jih 


d  r  /o  €/  (T‘  M  ,e)de  ] 
~dt  r  "  '  :  J 


(3) 


dt  dl  "  f,  f  (Tc,u,M,e)de 

When  the  electron-hole  system  is  treated  as  the  Maxwcll-Boltzmann  gas.  the  energy 

loss  rate  for  the  carriers  (c=e,hh)  is  given  by: 


Pr  =  K 


.  dTc  it  ) 


dt 


(4) 


It  should  be  pointed  that  the  energy  loss  rates  for  the  electrons  and  heavy-hole 
obtained  by  Eq.(3)  are  the  net  loss  rates.  These  include  the  energy  loss  due  to 
electron(hole)-phonon  interactions  and  the  energy  exchange  between  electrons  and 
holes  by  carrier-carrier  scattering.  The  hole-phonon  scattering  rate  was  calculated25 
and  measured6  to  be  2.5  to  3  times  larger  than  the  rate  for  the  electrons  due  to  the 
additional  coupling  through  the  deformation-potential.  The  energy  transfer  from 
electrons  to  holes  should  be  expected  in  our  photogenc rated  carrier  system  as  long  as 
the  characteristic  time  for  the  electron-hole  scattering  is  shorter  than  the  electron- 
phonon  scattering  time. 


The  determined  energy  loss  rates  using  Eqs.  3  or  4  are  shown  in  Fig.7.  The  solid 
and  dotted  curves  were  obtained  from  Eq.3  for  the  electrons  and  heavy-holes,  respec¬ 
tively.  The  dot-dashed  curve  was  obtained  from  Eq.(4).  It  should  be  pointed  out  that 
the  curves  obtained  from  either  Eqs-?  or  4  are  the  experimental  data  because  the  dis¬ 
tribution  function  f  (,Tc,ur/,)  was  experimentally  determined.  The  broken  line  in 
the  Fig.7  is  calculated  based  on  the  simple  theory5  where  the  lattice  is  treated  as  a 
heat  bath  for  the  quasi-equilibrium  carriers. 

The  salient  features  of  the  data  displayed  in  the  Eig.7  are: 


(i)  The  energy  loss  rate  for  the  electrons  is  more  than  two  orders  smaller  than 
predicted  by  the  simple  theory5  (broken  line).  The  reason  for  the  difference  arises 
from  the  re-heating  of  electrons  by  the  NE  phonons.  The  energy  loss  rate  determined 
for  the  electrons  by  scattering  with  LO-phonons  via  the  wave- vector  dependent 
Frohlich  interaction  may  be  even  smaller  than  what  we  have  determined  because  the 
net  energy  loss  rate  determined  from  Eq.3  is  sum  of  rates  due  to  electron-phonon 
interactions  and  electron-hole  interactions.  The  smaller  Phh  than  Pc  is  because  of  the 
presence  of  light-hole  population  and  the  different  masses  for  the  electron  and  hh 
which  leads  to  different  behavior  of  quasi-Fermi  energies  for  the  electron  and  hh  as 
function  of  time. 

(ii)  The  net  energy  loss  rates  for  the  electrons  and  heavy-holes  almost  follow  the 
rate  determined  from  Eq.(4).  The  reasons  for  this  are  the  rapid  decrease  of  ur  as 
shown  in  the  Fig.6  as  well  as  the  slow  cooling  of  the  carrier  temperature  which 
result  in  uc  /  KTC  <0.5,  making  the  electron-hole  system  behave  more  like  a 
Boltzmann  gas. 

(iii)  The  energy  loss  rates  increase  as  Tc  increases.  The  curves  in  the  Fig.7  even 
bends  over  and  approaches  the  maximum  value  of  5.3xlO_8\Vr  for  bulk  GaAs.  This 
"bending"  behavior  reflects  the  initial  rapid  cooling  of  hot  carriers  and  was  also 
observed  by  Shah  et  aL 6  This  may  be  the  indication  for  the  presence  of  NE  phonons  in 
the  system. 

The  average  phonon  emission  time  '  for  hot  carrier  to  emit  a  LO  phonon  is 
defined  by: 

f arg  =  e ~^L0 1  KTc  .  (5) 

The  value  of  rol,?  as  function  of  the  carrier  temperature  is  plotted  in  Fig.8.  This  is 
just  another  way  to  describe  the  energy  relaxation  process.  As  can  be  seen  the  Tavg  is 
not  a  constant  but  carrier  temperature  dependent  and  hence  time  dependent. 
For  Tc  >\200K  or  t  <  5ps  the  Tavg  remains  at  a  value  of  about  lps.  At  the  initial 
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stage  of  carrier  cooling  only  an  up-limit  of  0.9ps  was  determined  from  the  data 
shown  in  the  Fig.4.  The  theoretical  predicted  value  should  be  ~-0.16/>5  7  which  is 
consistent  with  our  value  within  the  limited  time  resolution.  As  time  t  >  5ps  (or 
Tc  <  1200AT  )  the  Ta,.s  increases  quickly  with  decrease  of  the  carrier  temperature 
reflecting  that  a  large  NE  phonon  population  is  built  up  and  the  carrier  cooling  is 
suppressed.  At  Tc  =  AWK  ,  r,Il?  is  lOps.  Ryan  ct  al .8  obtained  a  constant  value  of 
7ps  for  t ave  by  studies  of  time-resolved  PL  of  modulation  doped-MQW.  The  reason 
for  this  is  simply  because  their  time  resolution  of  '"'20ps  did  not  allow  them  to  mon¬ 
itor  the  entire  carrier  cooling  process. 

Information  about  the  time  dependence  of  initial  carrier  population  can  be 
obtained  using  the  experimentally  determined  distribution  function  f  .  The  densi¬ 
ties  are  obtained  by  the  following  expressions: 


f  °° 

nr(t  )  =  J()  P,-  f  cd  e  ; 

(6) 

nM  0  )  =  fft  Phh  f  h  d  e  ; 

(7) 

"m  ^  Pin  /WL 

(8) 

where  A E  is  the  energy  separation  between  the  hh-  and  lh-subbands  at  the  zone 
center.  The  carrier  densities  for  the  electrons,  heavy-holes,  and  light-holes  as  function 
of  time  are  plotted  in  Fig.9.  The  quantity  of  In  [n,.  (r  )]  is  also  plotted  in  this  figure  as 
a  broken  curve. 

The  salient  feature  of  the  data  in  the  Fig.9  is  that  the  carrier  density  decreases 
non-exponentially  and  varies  (by  a  factor  of  10)  rapidly  within  the  first  30ps  after 
the  excitation  by  laser  pulse.  One  can  not  define  a  single  carrier  lifetime  due  to  the 
non-exponential  nature  of  the  density  curves.  However,  an  effective  carrier  depletion 
time  [density  decreases  by  a  factor  of  e~‘  from  nr(t  =0)]  is  deduced  to  be  as  short  as 
lOps.  The  slow  component  of  density  decay  curve  reflects  part  of  bimolecular  recom¬ 
bination  process.  The  fast  decrease  of  carrier  density  occurring  in  such  a  short  time 
period  cannot  be  accounted  for  by  the  usual  bimolecular  recombination,  which  is  a 
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much  slower  process,  on  a  nanosecond  time  scale  in  a  bulk  GaAs  and  about  350ps  in 
quantum  well  structure  with  well  thickness  of  ~5nm.26  Several  other  possibilities 
should  be  pointed  out  and  discussed: 

(i)  Rapid  carrier  diffusion  process  may  be  utilized27  to  shorten  the  lifetime  of  the 
carriers  within  an  excited  spot.  Assuming  the  velocity  of  photoexcited  carriers  mov¬ 
ing  in  the  directions  perpendicular  to  the  z-axis  is  on  an  order  of  107cm  /  s  ,  an 
extremely  small  radius  of  about  1  um  for  the  excitation  spot  would  be  required  to 
decrease  carrier  density  by  a  factor  10  in  the  excited  region.  Since  the  radius  of  the 
excitation  spot  was  abcut  160um  in  our  experiment,  the  carrier  diffusion  out  of  the 
excited  region  can  not  result  in  the  ultrafast  depletion  time  of  lOps. 

(ii)  Auger  recombination  resulting  from  strong  electron-electron  Coulomb 
interaction  may  be  very  efficient  and  causes  a  substantial  shortening  of  carrier  life¬ 
time  at  high  carrier  density  in  a  quasi-two-dimensional  carrier  system.  But  from  a 
comparative  study,28  Auger  recombination  in  semiconductor  MQW  structure  and 
bulk  were  found  to  be  close  to  each  other.  Using  a  value  of  1.5xl0_29cmbr  for  the 
nonradiative  Auger  rate  in  bulk  GaAs  at  room  temperature.29  Auger  recombination 
lifetime  at  carrier  density  of  1019cm'3  should  be  b70ps  which  is  about  100  times  larger 
than  the  carrier  depletion  time.  However,  a  detailed  study  of  Auger  recombination  in 
semiconductor  MWQ  structures  under  the  condition  of  high  carrier  density,  high  car¬ 
rier  temperature,  low  lattice  temperature,  and  presence  of  a  large  population  of  Mi 
LO  phonons  does  not  exist  in  the  literature. 

(iii)  Finally,  let  us  consider  what  role  the  NF  phonons  can  play  in  reducing  the 
initial  carrier  depletion  time  to  lOps.  As  discussed  in  the  steady-state  section,  there  is 
a  stimulated  phonon  emission  process  accompanied  by  an  intense  stimulated  phonon 
replica  which  is  mainly  radiated  at  the  sample  edge  and  propagates  parallel  to  the 
well  plane.  We  found  that  the  stimulated  phonon  replica  detected  along  the  y- 
direction  also  had  a  decay  time  constant  about  30ps  which  is  same  as  that  detected 
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along  the  z-direction.  Therefore,  it  is  most  likely  that  the  initial  rapid  carrier  deple¬ 
tion  is  due  to  the  participation  of  the  NT  LO  phonons  which  stimulates  the  emission 
of  phonon  replica.  It  should  be  noted  that  this  stimulated  process  is  remarkablely 
different  from  the  regular  stimulated  emission  process  in  which  phonons  do  not  parti¬ 
cipate.  As  in  semiconductor  lasers,  the  stimulated  emission  is  less  effective  at  room 
temperature.  Moreover,  a  longitudinal  optical  phonon  (LO)-phonon  replica  plays  no 
significant  role  in  the  laser  operation  of  bulk  GaAs.  which  is  expected  because  of  the 
weak  electron- phonon  coupling.  Therefore,  the  prominence  of  phonon-assistant 
recombination  in  our  quantum  well  structure  is  due  to  the  presence  of  a  large 
number  of  NH  LO- phonons  produced  in  the  relaxation  process  of  hot  carriers.  By 
employing  the  N E- phonon  model  we  able  to  explain  not  only  the  slow  car  rier  cixUing 
hut  also  the  ultrashort  carrier  depletion  time  m  highly  photoexcited  undoped  quantum 
well  structur  es. 

V  Conclusion 

With  2ps  time  resolution  we  are  able  to  study  the  initial  energy  relaxation  of 
the  hot  carriers  and  the  decrease  of  carrier  density  simultaneously.  I  he  existence  of  a 
large  population  of  NE  phonons  in  highly  excited  semiconductor  quantum  well  struc¬ 
tures  is  experimentally  verified.  Its  effect  on  the  energy  relaxation  is  to  slow  down 
the  cooling  rate  after  an  initial  rapid  cooling  (0-5ps).  A  new  mechanism  used  to 
explain  the  ultrashort  carrier  depletion  time  —  \()ps  deduced  from  the  fitting  of 
time-resolved  PL  profiles  at  different,  emitted  photon  energies  is  proposed  to  be  associ¬ 
ated  with  NE  phonon  stimulated  phonon  replica  emission  rather  than  other  nonradia- 
tive  processes. 
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Figure  Caption 

Fig.  1  Tiine-integrated  luminescence  spectra  at  various  lattice  temperatures.  The  peaks 
A,  B,  and  C  are  explained  in  the  text. 

Fig.2Triangles  are  the  ratios  of  the  peak  intensity  of  B  and  A  indicated  in  the  Fig.  1 . 
The  solid  curve  is  explained  in  the  text. 

Fig.JTiine-integrated  PI.  spectra  taken  in  the  geometry  as  shown  in  the  inset  at  100K. 

Fig.TTime-resolved  PL  profiles  (thin  solid  curves)  from  the  symmetric  MQW  at  4K  at 
various  wavelengths.  The  narrow  peaks  within  0-1  Ops  are  the  pre-pulses  used 
for  averaging.  The  next  peak  of  dotted  curve  is  the  Rayleigh  scattering  light 
from  the  sample  surface.  The  thick  broken  curve  is  generated  by  an  expression 
of  I(Ejt)=I<{ l-e  ,/T'V  ' ' Tj  and  the  thick  solid  curves  are  generated  by 
Eq.(  1 ). 

Fig.5 Experimentally  determined  carrier  temperature  as  function  of  time.  The  shaded 
area  indicates  the  extent  of  uncertainty  in  deducing  carrier  temperature  within 
first  4ps. 

Fig.6Experimental  determined  quasi-Fermi  energies  for  electrons  (up  solid  curve)  and 
for  holes  (down  solid  curve). 

Fig.7Experimentally  determined  net  energy  loss  rates  for  electrons  (e)  and  heavy- 
holes  (hh)  obtained  from  two  quasi-equilibrium  Fermi-Dirac  distribution  func¬ 
tions.  The  dot-dashed  curve  is  obtained  from  Maxwell-Boltzmann  distribution 
function.  The  broken  line  is  theoretical  result  for  bulk  (iaAs  without  the  effects 
of  screening  and  nonequilibrium  phonons.  The  power  unit  is  watts. 

Fig.SExpcrirnentally  determined  average  phonon  emission  time  as  function  of  carrier 
temperature. 

Fig.9Experimentally  determined  carrier  densities  as  function  of  time,  e  for  electron, 
hh  for  heavy-hole,  and  Ih  for  light-hole.  The  broken  curve  is  the  nature  loga- 
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ABSTRACT 

Lateral  expansion  of  energetic  carriers  excited  by  an  intense  0.5-ps 
laser  pulse  in  a  single  asymmetric  GaAs  quantum  well  is  investigated 
at  4.3K  using  a  streak  earn  ,ra.  On  a  picosecond  time  scale,  it  is 
observed  that  the  carriers  expand  from  the  excitation  spot  at  an 
ultrafast  speed  of  ~  10 *cm  /  s.  A  comparison,  under  identical  condi¬ 
tions,  of  the  asymmetric  and  symmetric  quantum  well  results  sheds 
light  on  the  mechanism  that  leads  to  such  an  extraordinary  large  car¬ 
rier  expansion  velocity. 

PACS  numbers:  66.30.Dn,  72.20.  Jv.,  73.40.  Lq. 


Over  several  years,  spatial  expansion  of  photoexcited  carriers  in  semiconductor 
bulk  and  microstructures  has  been  a  subject  of  great  interest  Based  on  the  photo- 
luminescence  (PL)  data,  Forchel  ct  al .'  have  suggested  that  the  fast  ambipolar 
diffusion  of  optically  excited  carriers  in  a  wide  varity  of  semiconductors  (e.g.,  Si,  Ge, 
GaAs,  CdS,  and  ZeSe)  at  a  speed  comparable  to  Fermi  velocity  is  essential  for  the 
understanding  the  broad  line  shape  of  the  luminescence  spectra.  Later,  Steranka  and 
Wolfe2  pointed  out  that  the  fast-diffusing  model1  overestimated  carrier  expansion 
velocity  by  an  order  of  magnitude  on  nanosecond  time  scale  in  indirect-gap  Si  due  to 
neglecting  the  local  lattice  heating  effect.  Subsequently,  Tsen3  confirmed  the  results 
of  Steranka  and  Wolfe"  on  a  picosecond  time  scale  by  time-resolved  Raman  scattering. 
Using  excitc-probe  technique,  Majumder  ct  al4  found  the  carrier  expansion  velocity 
was  smaller  than  Fermi  velocity  in  the  direct-gap  CdS  and  CdSe.  In  contrast,  many 
groups  deduce  from  their  experiments  carrier  expansion  velocity  vr  ^  107cm  /  s  in 
GaAs5-6  and  C'dSe.7,8  This  value  is  roughly  equal  to  the  Fermi  velocity  of  carriers  and 
exceeds  group  velocities  of  longitudinal  and  transverse  acoustic  phonons  by  about  2 
orders  of  magnitude,  thus,  supporting  the  fast-diffusing  model  of  Forchel  ct  al}  in 
these  direct-gap  polar  materials.  Theoretically,  Wautelet  and  Van  Vechten9 
have  inferred  that  hot  carriers  generated  by  nanosecond  laser  pulses  will  be 
confined  near  the  excitation  spot;  while  a  spatial  expansion  of  photoexcited  car¬ 
riers  is  expected  from  both  hydrodynamic5’10  and  thermodiffusive11  points  of 
view.  Recently,  Steranka  and  Wolfe2  have  suggested  that  nonequilibrium 
(NE)  phonon  effects12  (phonon  wind)  must  be  taken  into  account  in  a  more 
realistic  description  of  the  dynamic  carrier  expansion  process.  The  NE  phonons 
created  in  hot -carrier  relaxation  process  provide  a  large  driving  force  for  car¬ 
rier  expansion. 

In  GaAs/AlGaAs  system,  the  band-discontinuities  at  the  interfaces  of  GaAs  and 
AlGaAs  leads  to  the  formation  of  two-dimensional  electronic  energy  subbands  in  the 
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GaAs  region.  Recent  Raman-scattering  studies  by  Sood  ct  al.  3  provided  evidence  for 
the  presence  of  interface  phonon  modes  and  confined  longitudinal  optical  (LO)  and 
transverse  optical  (TO)  phonons  in  GaAs  well  layer.  One  would  like  to  kown  how 
carrier  confinement,  modified  phonon  modes,  and  well  shape  influence  the  dynamic 
process  expansion  of  photoexcited  carriers  in  a  quantum  well  structure.  Most 
recently,  Tsen  and  Morkoc14  have  shown  no  influence  of  the  effects  of  reduced  dimen¬ 
sionality  for  carriers  and  phonons  on  photogenerated  carrier  expansion  in 
GaAs/AlGaAs  symmetric  MQW  structures  using  picosecond  time-resolved  Raman 
scattering  at  a  low  carrier  density  of  ~-10!(jcm_i.  Their  measured  expansion  velocity 
of  8.6x1 06cm  /  s  is  consistent  with  the  fast-diffusing  model  of  Forchel  ct  alJ 

In  this  letter,  we  report  on  the  first  real  time  investigation  of  ultrafast  hot  car¬ 
rier  expansion  process  associated  with  an  ultrashort  carrier  depletion  time  of  ~ 
lOps  in  an  asymmetric  GaAs  quantum  well  at  4.TK  using  a  2-ps-resolution  streak 
camera  system.  By  studies  of  time-resolved  PL  arising  from  the  recombination  of 
energetic  electron  and  hole  detected  from  a  lateral  displaced  region  ~  MtOum  away 
from  the  center  of  the  excitation  spot  which  is  ~  160 um  in  diameter,  it  is  observed 


that  the  hot  carriers  expand  from  the  excitation  spot  at  an  ultrafast  speed  of  ~ 
108cm  /  s.  Through  a  comparison  of  the  symmetric  and  asymmetric  quantum 
well  results  under  identical  conditions,  it  is  concluded  that  the  observation  of  such  a 
fast  hot  carrier  expansion  process  is  made  possible  only  in  the  asymmetric  well 
with  the  higher  energy  barrier  such  as  to  prevent  from  carrier  expansion  in  the  well 
direction.  Density  and  cooling  of  photogene  rated  carriers  are  independently  deter¬ 
mined  by  analysis  of  time-resolved  PL  from  the  excitation  spoL  We  have  demon¬ 
strated  that  cooling  of  hot  carriers  is  severely  slowed  by  the  NT  phonon  heating  effect 
and  the  carrier  density  decays  very  rapidly  within  the  first  TOps  after  the  excitation 
by  0.5ps  laser  pulse  due  to  the  NL  phonon  stimulated  phonon  replica. 


In  our  experiments,  an  intense  laser  pulse  of  0.5ps  at  (>20nm  was  used  to  excite 
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the  electron-hole  pair  density  of  ~  1014cm~3.  The  excitation  source  is  from  a 
colliding-pulse  passive- mode- locked  dye  laser  amplified  by  a  four-stage  dye  amplifier 
pumped  with  a  frequency  doubled  Nd:  YAG  laser  operated  at  20Hz.  The  undoped 
GaAs/Al^Gaj  xAs  asymmetric  and  symmetric  quantum  wells  investigated  in  this 
work  were  grown  by  molecular-beam  epitaxy  on  a  (001  Loriented  undoped  GaAs  sub¬ 
strate.  The  asymmetric  quantum  well  consists  of  layers  from  the  front  surface  to 

o  o  o 

substrate  150A  Al03Ga07As,  50A  GaAs  well,  200A  AlAs  (asymmetric  barrier),  and 
followed  by  1  um  Al03Ga07As,  and  0.2 um  GaAs  buffer.  The  symmetric  well  struc- 

o  o 

ture  used  for  the  comparative  studies  consists  of  50  periods  of  55A  GaAs  and  100A 
Alo3GaQ7As  layers  followed  by  a  1.2 um  GaAs  buffer  layer.  The  excellent  quality  of 
the  samples  are  confirmed  by  well  resolved  heavy-hole  and  light-hole  excitonic  struc¬ 
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tures  in  room  temperature  PI  spectra. 

First,  we  discuss  the  measurements  of  the  differences  in  the  spatial  widths  of 
time-integrated  laser  and  PL.  intensity  profiles  on  the  sample  surface.  I  he  hoi  PI. 
profiles  (solid  curves)  detected  within  a  spectral  window  of  35nm  (‘JOmeV  )  centered 
at  700nm  (1.771eV)  for  the  asymmetric  and  symmetric  wells  along  with  the  profiles 
(broken)  of  corresponding  exciting  laser  at  a  fixed  intensity  are  plotted  in  Fig. la  and 
Fig. lb,  respectively.  I  he  gain  of  detecting  system  was  adjusted  such  that  the  max¬ 
imum  intensities  recorded  for  the  PL  and  laser  wrere  roughly  same.  The  insets  show 
the  potential  profiles  of  the  asymmetric  and  symmetric  wells.  It  is  clear  that  the 
half-width  of  PL  profile  for  the  asymmetric  well  is  much  broader  than  for  the  sym¬ 
metric  wells.  This  implies  that  the  energetic  carriers  (about  3  LO- phonon  energies 
above  the  subband  edge)  in  the  asymmetric  well  is  much  more  mobile  than  in  the 
symmetric  well. 


< 

eS 


f 


v.V 


By  making  use  of  the  two  electronic  windows  of  the  streak  camera,  time- 
resolved  PL  intensities  of  the  asymmetric  well,  Iw j(t)  and  Iw.2(t),  from  the  excited 
and  unexcited  regions,  respectively,  were  measured.  The  subscripts  wl  and  w2 
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(r  i~*r  t)  stand  for  the  two  windows  as  being  positioned  in  the  Fig.  1.  The  lu.j(t)  and 
Iw2(t)  (thin  solid  curves)  measured  at  4.4K  with  ~2ps  time  resolution  are  displayed 
in  F;ig.2.  The  peak  of  the  broken  curve  is  the  Rayleigh  scattering  light  ((>20nm)  from 
the  sample  surface  which  defines  the  "zero"  time  for  our  analysis.  The  thick  solid 
curves  are  the  theoretical  fits  which  will  be  discussed  below.  As  can  tx-  seen  the  rise 
of  the  lu.j  is  time-: 'solution  limited,  while  the  rise  of  1^.,  is  much  slower  than  that 
of  the  Iu. j.  The  slow  rise  time15  of  1^.,  (r,“  2  ^  lOOps)  is  essential  to  suppoort  that 
the  luminescence  detected  from  the  laterally  displaced  w2  indeed  arises  from  the 
recombination  of  electron  and  hole  whioh  come  from  the  wl.  It  eliminates  one  possi¬ 
ble  artifact  that  may  lead  to  anomalously  strong  luminescence  detected  from  the  w2. 
Such  a  possibility  arises  from  the  1  um  layer  of  AI(nGa()7As  which  is  next  to  the 
AlAs  barrier.  This  layer  is  also  highly  excited  by  the  (>2()nm  laser  pulse.  Possible 
recombination  radiation  (although  we  did  not  detect  for  the  asymmetric  well,  neither 
for  the  symmetric  well)  from  Al()^Ga()7As  layer  will  excite  carriers  within  the  well. 
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This  secondary  excitation  would  not  be  limited  to  the  originally  excited  region  giving 
rise  secondary  luminescence.  Since  the  speed  of  light  in  GaAs  is  about  lOOum  /ps,  the 
luminescence  intensity  Iw2  would  rise  quite  similarly  with  the  lwJ.  This  is  contrary 
to  our  experimental  observation. 

Having  known  the  expansion  distance  ( 294i im  )  from  the  center  of  the  w  1  to  the 
inner  edge  of  the  w2  and  assuming  the  time  it  takes  for  carriers  to  expand  from  the 
wl  to  w2  where  they  recombine  is  equal  to  the  rise  time  of  1  v  an  up-limit  of  the 
expansion  velocity  v,  is  estimated  to  be  about  2.94xl()scm  /  j  (294um  /  lOOps). 

The  decay  of  the  lu  ]  is  dependent  upon  a  transient  carrier  distribution  in  the  wl 
which  determines  the  number  of  carriers  in  the  given  energy  space  as  function  of 
time.  In  order  to  obtain  the  carrier  distribution  (assumed  to  be  lermi-Dirac  distribu¬ 
tion),  the  time- resolved  PI.  profiles  detected  in  narrow  (20meV)  spectral  regions 


which  cover  the  whole  emission  band  were  measured  from  the  region  wl  at  4.3K. 
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We  have  introduced  an  expression16  using  the  carrier  temperature  [7V  it  )]  and  the 
quasi-Fermi-energy  for  electrons  [i^c  (r  )]  as  a  set  of  parameters  to  consistently  fit  these 
measured  profiles,  three  typical  measured  profiles  as  well  as  theoretically  fitted 
curves  (thick  solid)  are  plotted  in  Fig.Ta.  The  subband  edge  luminescence  is  centered 
~760nm.  The  cooling  of  hot  carriers  is  substantially  slowed  down  after  ~5ps  due 
to  reabsorption  of  NH  phonons  by  the  carriers.16"18  The  electron  density  nc  (r  )  and 
In  [n,.  (r  )]  are  plotted  in  Fig.Jb.  The  salient  feature  of  the  plotted  data  is  that  the  car¬ 
rier  density  decreases  nonexponent ially  and  very  rapidly  form  0  to  BOps.  One  can  not 
define  a  single  lifetime  to  describe  the  carrier  density  decay.  However,  an  effective 
carrier  depletion  time  [the  time  for  n,  to  decrease  by  a  factor  of  e_1  from  n,At  =(>)] 
can  be  deduced  to  be  as  short  as  lOps.  Such  an  extremely  rapid  carrier  density 
decrease  can  not  lx-  accounted  for  by  an  usual  bimolecular  recombination  which  is  in 
nanosecond  time  scale  for  a  bulk  Ha  As  and  about  3.M)ps  lor  .1  quantum  well  structure 
with  well  thickness  of  ~5()A.1<y  But  it  can  be  explained  consistently  b\  the  participa¬ 
tion  of  NT  phonons  in  the  highly  excited  system.  Since  NT  phonons  aie  local  1/ed  both 
in  a  wavevector  space  and  in  a  real  physical  space  (in  the  well  plane)  they  will 
behave  like  coherent  Bosons.  The  existence  of  the  coherent  liosons  will  tend  to 
increase  the  number  of  the  Bosons  in  the  system  with  a  rate  proportional  to  the 
present  number  of  Bosons.  This  is  a  stimulated  phonon  emission  process  and  can  be 
accomplished  by  a  phonon  replica,  when  an  electron  recombines  with  a  hole,  not  only 
a  photon  but  also  a  phonon  or  more  phonons  will  be  emitted  to  join  Boson  system. 
Such  a  NT  phonon  stimulated  phonon  replica  has  been  reported  previously  by 
Holoayak  ct  aL20  We  have  recently  confirmed16  this  novel  process  in  the  symmetric 
wells  by  monitoring  the  phonon  replica  which  is  mainly  radiated  in  the  lateral  x-y 
plane  in  time  domain  as  well  as  detailed  polarization  studies  of  integrated  lumines¬ 
cence  spectra. 

Using  the  determined  carrier  distribution,  the  solid  curve  as  shown  in  Tig. 2a  is 


theoretically  generated  to  fit  the  measured  Iu.r  Good  fitting  indicates  that  the  photo- 
generated  carrier  system  is  not  a  purely  diffusive  system. 


We  now  discuss  possible  origins  for  the  observed  anomalously  large  expansion 
velocity  of  the  energetic  carriers  in  the  asymmetric  quantum  well  but  not  in  the 
symmetric  quantum  well. 

According  to  the  thermodiffusive  model,1,11  a  classical  diffusion  equation7  can  be 
employed  to  fit  I  ,,(t).  Assuming  the  initial  photogenerated  carrier  profile  is  a  sym¬ 
metric  bell-shaped  function7  the  carrier  density  as  function  of  time  t  and  distance  r 
can  lx-  written  as 


n  (  r  ,  t  )  =  nlte  ' 1  ?(l  +  ]U°g  \  D  ’  )-1/  2  e  "  2>  «'<*  2  +  •»/>  ’ '  , 


(1) 


where  D  is  an  effective  diffusion  constant  for  the  energetic  carriers  in  the  given 
spectral  region  and  time  scale  of  20ps,  n„  is  the  carrier  density  at  the  center  of  laser 
spot,  and  a  is  a  width-parameter  for  the  initial  carrier  profile  which  is  di  ternuned  to 
be  12 Oum .  When  electrons  and  holes  move  together  in  the  lateral  x-y  plane  of  the 
asymmetric  quantum  well,  i.e.  ambipolar  diffusion  applies  to  the  present  case,  the 
luminescence  detected  from  the  w2  should  be  given  by: 


/  u  2  U  )  ~~  f  n‘(  r  ,t  )  dr  .  (2) 

Using  the  alxive  expression  to  fit  the  measured  I  v  two  parameters  of  /)  and  rff  2 
have  Ixcn  precisely  determined  to  be  0.7x1  (/’cm  -/  a  and  13ps,  rcspc'ctively. 

It  should  be  investigated  whether  the  determined  values  for  D'  and  t -  can  be 
interpreted  from  the  thermodiffusive  model.1,1 1  Several  aspects  are  discussed  below: 

(1)  The  value  of  /)'  is  four  orders  magnitude  greater  than  the  regular  impurity 
scattering  limited  diffusivity.  The  enhancement  of  /)'  may  arise  from  high  carrier 

o 

temperature  and  the  screening  effects  at  high  carrier  density.  However,  the  screen¬ 
ing  effects  do  not  seem  to  be  effective  since  the  initial  carrier  cooling  is  very  rapid.16 
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(2)  The  decay  rates  of  carriers  which  are  located  in  the  same  energy  space  but  in 

the  different  spatial  regions  of  the  w  1  and  w2  are  comparable  reflecting  the  dominant 

process  for  carrier  density  decay  is  not  due  to  the  carrier  expansion  out  from  the  wl 

but  due  to  the  NE  phonon  stimulated  phonon  replica.  It  is  further  inferred  that  a 

large  number  of  the  NT  phonons  emitted  by  hot  electrons  may  exist  in  an  area  of  at 

least  ~  300 um  in  diameter.  These  phonons  which  have  momenta  parallel  to  the  x-y 

plane  relax  to  acoustic  phonons  of  large  momentum  and  then  these  acoustic  phonons 
-■) 

may  drive-  the  carriers  outside  the  excited  region.  A  rough  estimation  seems  support 
this  argument.  Let  us  take  a  value  of  10~8  dynes  per  pair2  for  acoustic  phonon-wind 
force,  an  excess  energy  of  140meV  per  electron-hole  pair,  and  an  acceleration  time  of 
Lips.  Then,  a  value  of  2.8xlOscm  /  s  for  carrier  expansion  velocity  is  predicted  by 
phonon-wind  model.2,12 

(3)  From  the  determined  D  and  tJ  2,  the  expansion  velocity  is  calculated  by 
-  —  =  2.3xl(),cm  /  \  which  is  consistent  to  the  thermocliffusiw  model1,  1 


VI 


and  the  estimated  up- limit.  In  a  two-dimensional  electron  system  quasi-Fcrmi-energy 
is  proportional  to  carrier  density.  Therefore,  the  Fermi- velocity  vF  and  hence  v,  is 
proportional  to  the  square  root  of  the  carrier  density.  Then,  our  observed 
v,  =  2~-3xlOscm  /  s  at  the  high  carrier  density  of  ~  10iqcm-3  is  consistent  with  the 
measured  v,  =  S.ftxH/’cm  /  s  by  Tsen  and  Morkoc14  at  the  !o\e  carrier  density  of  ~ 
10  lhcm  -3. 

(4)  Using  the  determined  I)'  and  r/2,  an  effective  "diffusion  length"  is  calcu¬ 
lated  to  be  30 um  .  According  to  this  value,  there  would  be  no  significant  carriers  in 
the  region  of  w2  which  is  ~3 OOum  away  from  the  excitation  center.  But  the  total 
detected  photons  from  the  wl  and  w2  wuthin  the  first  20ps  after  the  excitation  by 
the  laser  pulse  is  comparable.  This  might  be  again  interpreted  by  the  phonon-wind 
model2  in  which  the  concept  of  conventional  diffusion  length  may  not  apply  since 
the  collisions  between  carriers  and  phonons  no  longer  limit  how  far  carriers  can  be 
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drifted. 

(5)  Finally,  it  must  be  noted  that  the  anomalously  fast  expansion  is  only 
observed  in  the  asymmetric  well.  This  is  very  interesting  because  such  a  simple 
modification  of  quantum  well  shape  in  comparison  with  the  symmetric  well  makes  a 
big  difference  on  the  dynamic  carrier  expansion  process.  Due  to  the  asymmetric  well 
shape  there  are  no  transmission  states  (virtual  states  or  unconfined  states)  above  the 
lower  energy  barrier  unlike  the  symmetric  well  case.  This  is  an  intrinsic  property  of 
the  asymmetric  quantum  well.  The  photogenerated  carriers  may  be  immediately 
scattered  by  the  higher  energy  harrier  and  obtain  momenta  from  the  interface  which 
assists  carriers  to  expand  rapidly  in  x-y  plane. 

In  summary,  we  have  observed  for  the  first  time  the  ultrafast  expansion  of  the 
energetic  carriers  in  the  asymmetric  quantum  well.  It  is  found  that  the  unique  well 
shape  of  the  asymmetric  quantum  well  is  responsible  lor  the  observation  of  the 
extraordinary  large  expansion  velocity  of  the  energetic  carriers.  However,  no  distinc¬ 
tion  has  been  made  among  various  models  to  explain  the  observed  ultrafast  carrier 
expansion  process.  Both  the  rapid  carrier  density  decay  and  the  slow  cooling  of  hot 
carriers  in  the  asymmetric  and  symmetric  wells  arc  attributed  to  the  presence  of  the 
nonequilibrium  phonons  generated  by  the  energy  relaxation  of  photoexcited  carriers. 

The  research  was  supported  by  the  U.S.  Air  Force  Office  of  Scientific  Research 
under  Grant  No.  AFOSR-86-OOT1.  We  thank  K.  Bajaj  for  helpful  discussions. 
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Figure  Ca  ption 

Fig.  1  Time-integrated  spatial  profiles  of  laser  and  the  luminescence  in  a  spectral  win¬ 
dow  of  35nm  centered  at  700nm  at  4.3K  from  the  asymmetric  well  (a)  and 
symmetric  wells  (b).  The  wl  ad  w2  are  regions  where  the  time-resolved 
luminescence  profiles  are  obtained.  The  insets  show  the  potential  profiles  of  the 
asymmetric  and  symmetric  wells. 

Fig.2 Time-resolved  photoluminescence  of  the  asymmetric  well  at  4.3K  from  the  spa¬ 
tial  region  of  the  wl  (a)  and  w2  (b).  The  solid  curves  are  the  theoretical  fits  as 
discussed  in  the  text. 

Fig.3(a)  Time-resolvecl  photoluminescence  from  the  asymmetric  quantum  well  at 
4.3K  at  various  wavelengths,  (b)  The  electron  density  decay  curves. 
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Electron-hole  recombination  lifetimes  in  a  quasi-zero-dimensional  electron 
system  in  CdS^Se^, 
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The  recombination  lifetimes  for  the  radial  and  angular  quantum  number  conserved  15-15  and 
\P-\P  transitions  from  three-dimensionally  confined  electrons  in  CdS,  Se,  _  x  were  measured 
by  time-resolved  photoluminescence  (PL).  The  assignment  of  the  observed  transitions  was 
supported  by  calculations  of  eigen  energy  levels  and  squared  matrix  element  ratio  for  these 
transitions  as  well  as  well-resolved  PL  peaks  arising  from  15-15  and  \P-\P  transitions. 


Recently,  carriers  localized  in  semiconductor  micro- 
structures  have  attracted  much  attention  because  of  their 
novel  optical  properties  for  potential  device  applications. 
Electrons  localized  in  semiconductor  crystallites  of  diame¬ 
ters  ranging  from  30  to  800  A  embedded  in  a  transparent 
insulating  matrix  are  confined  in  three  dimensions.1'3  Such 
quasi-zero-dimensional  (quasi-OD)  electron  systems  also 
exist  in  bulk  alloy  semiconductors  due  to  compositional  fluc¬ 
tuations.4  Electronic  motion  in  these  systems  no  longer  fol¬ 
lows  a  well-defined  energy  momentum  relation  because  the 
Hamiltonian  of  the  system  does  not  commute  with  the  mo¬ 
mentum  operator  due  to  potential  discontinuity  at  the  ciys- 
tallite  surface.  The  envelope  wave  functions  in  spherical  co¬ 
ordinates  [  r,0,d>  ]  and  eigen  energies  of  conduction  electrons 
and  valence  holes  localized  in  an  infinite  spherical  well  with¬ 
in  the  effective  mass  approximation  are  given  by 

\nLm)  =  C„JL  (*„tr)  Yl„  (0,«>)  ( 1 ) 

and 


=*Tw./2  mf,ha2. 


(2) 


respectively,  where  the  subscripts  n,  L,  and  m  are  effective 
radial,  angular,  and  magnetic  quantum  numbers,  respective¬ 
ly;  the  j  and  Y  are  the  spherical  Bessel  and  spherical  Har¬ 
monic  functions,  respectively;  a  is  the  radius  of  crystallite; 
and  mrh  is  the  effective  mass  of  electron  or  the  isotropic  hole 
mass.3  The  values  of  x*l  for  the  lowest  two  states  of  either 
conduction  electron  or  valence  hole  15  and  IP  are  rr  and 
4.49,  respectively.  5  stands  for  L  =  0  and  P  for  L  =  1.  The 
allowed  transitions  which  conserve  angular  and  radial  quan¬ 
tum  numbers  are  15-15, 1P-1P,  and  higher  transitions.  The 
physical  picture  of  quasi-OD  electron  system  described 
above  has  been  experimentally  verified  by  several  groups.1'3 
Large  and  fast  optical  nonlinearities  arising  from  photogen¬ 
erated  electron  hole  pairs  have  been  reported  in  quasi-OD 
electron  systems.6  7 

In  this  letter,  we  report  on  measurements  of  ultrashort 
recombination  lifetimes  of  the  radial  and  angular  quantum 
number  conserved  15-15and  \P-\P  transitions  in  quasi-OD 
electron  systems  in  CdS,  Se,  _  ,  at  4.3  K  using  a  streak  cam¬ 
era  detection  system.  The  time-resolved  photoluminescence 
(PL)  detected  at  various  emission  energies  allows  us  to  un¬ 
ambiguously  identify  the  15- 15  and  1P-1P  transitions.  The 
recombination  lifetime  of  \P-\P  transition  was  measured  to 
be  3.5  times  shorter  than  the  15- 15  transition.  The  ultrafast 
decay  of  the  IPexcitation  may  have  practical  importance  for 
the  construction  of  ultrafast  reversible  optical  switches.  The 


assignment  of  the  observed  5  and  P  transitions  was  support¬ 
ed  by  calculations  of  eigen  energy  levels  and  squared  matrix 
element  ratio  for  these  transitions,  as  well  as  the  observation 
of  optical  transitions  in  steady-state  PL. 

The  samples  investigated  were  four  optical  glass  filters 
Coming  2-61,  2-59,  2-58,  and  2-64.  The  value  of  x  for  each 
sample  was  accurately  obtained  from  chemical  analysis.8  A 
second  harmonic  ( 530  nm)  of  a  Nd-glass  laser  pulse  of  8  ps 
duration  was  used  to  excite  the  samples  on  the  front  surface. 
The  maximum  optical  energy  incident  onto  the  front  surface 
was  about  40 /zJ.  The  spot  size  was  about  8x  10~3  cm:.  For 
steady-state  PL  experiments,  an  argon-ion  laser,  a  Spex  dou¬ 
ble  grating  spectrometer,  an  S-20  photomultiplier,  and  a 
lock-in  amplifier  were  used. 

The  time-resolved  PL  profiles  obtained  at  various  emis¬ 
sion  energies  at  4.2  K  for  sample  2-58  are  shown  in  Figs. 
1(a)— 1  (e).  The  left  curve  in  Fig.  1(a)  shows  the  temporal 
profile  of  the  exciting  laser  pulse  which  reflects  the  time  reso¬ 
lution  of  the  detection  system  ( 10  ps).  The  dotted  curves  in 
Fig.  1  are  the  double-exponential  fits  to  the  data  with  a  value 
of  12  ps  for  the  rise  time  r, .  The  decay  times  of  these  time- 
resolved  luminescence  profiles  show  only  two  distinct  val¬ 
ues.  For  emission  energies  ranging  from  2.213  to  2.175  eV 
and  from  2.1  to  1.967  eV  the  decay  times  are  29  and  100  ps, 
respectively. 

The  decay  time  as  function  of  emitted  photon  energy  E 
for  four  quasi-OD  electron  systems  at  4.2  K  is  summarized  in 
Fig.  2.  For  a  comparison,  the  exciton  lifetime  versus  exciton 
energy  in  a  quasi-OD  system  in  a  bulk  CdSosjSe^  alloy 
compound  studied  by  Kash  el  a l.9  is  also  included  in  this 
figure.  The  most  remarkable  feature  of  interest  in  the  data  is 
the  appearance  of  two  distinct  energy  regions  in  which  the 
decay  time  is  different  by  a  factor  of  about  3.5  for  the  samples 
2-64,  2-58  and  the  bulk  CdSojjSeo,,  while  the  other  two 
samples  do  not  show  this  feature  because  either  IP  electron 
and  hole  were  not  excited  or  the  luminescence  from  the  1P- 
1P  electron-hole  recombination  was  not  detected.  It  is  ex¬ 
pected  that  the  decay  time  should  be  nearly  invariable  over 
the  energy  range  of  emission  if  only  the  lowest  confined  state 
15  for  electrons  and  holes  is  occupied.  The  wide  spectral 
range  most  likely  reflects  the  crystallite  size  distribution  and 
the  fluctuations  in  the  value  of  x  from  crystallite  to  crystal¬ 
lite.  When  the  two  lowest  states  15  and  IP  are  substantially 
occupied  and  the  recombination  lifetimes  associated  with 
these  two  states  are  considerably  different,  a  steplike  change 
of  decay  time  will  be  observed  in  two  distinct  energy  regions. 
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FIG  1.  Time-resolved  PL  profiles  de¬ 
tected  at  different  energies  (a)-(e)  for 
the  sample  2-58  at  4  2  K.  The  dotted 
curves  are  the  theoretical  fit  with  1 2  ps 
for  rise  time  and  decay  times  are  indi¬ 
cated  on  the  corresponding  profiles. 
The  left  curve  of  (a)  is  the  temporal 
response  profile  to  the  laser  pulse  used 
to  excite  the  samples. 
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The  exact  lifetime  ratio  in  the  two  energy  regions  depends  on 
the  transition  matrix  elements.  The  above  argument  should 
in  principle  explain  what  we  observed.  However,  the  energy 
dependence  of  localized  exciton  lifetime  in  the  bulk  com¬ 
pound  CdSo  jj  Se<j47  as  displayed  in  Fig.  2(e)  also  exhibits  a 
steplike  feature.  The  explanation  given  by  Kash  et  al ,9  was 
based  on  the  model  suggested  by  Cohen  and  Sturge4  where 
the  exciton  migration  goes  from  a  site  with  higher  energy 
into  another  site  with  lower  energy.  This  raises  a  question 
whether  the  exciton  migration  mechanism  can  apply  to  the 
present  case.  Since  exciton  migration  must  involve  a  trans¬ 
port  over  a  relatively  large  distance  due  to  the  nature  of  the 
exciton-phonon  interaction,4  such  migration  may  occur  in 
relatively  large  crystallites.  This  apparently  may  explain  the 
results  observed  in  samples  2-58  and  2-64.  In  smaller  crystal¬ 
lites  the  migration  does  not  occur  in  correspondence  with 
the  samples  2-61.  However,  the  exciton  migration  picture 
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FIG.  2.  Lifetimes  »re  plotted  as  function  of  emitted  photon  energy.  The 
sample  temperature  for  (a),  (b),  (c),  (d)  is  4.2  K  and  for  (e)  is  2  K.  The 
arrows  indicate  the  energy  positions  for  the  IS- 15  transition  (5)  and  the 
\P-\P  transition  (P). 
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cannot  explain  the  steplike  energy  dependence  oflifetime  as 
well  as  the  ratio  with  the  value  of  about  3.5  because  the  depth 
of  a  potential  well  caused  by  compositional  fluctuations  is 
entirely  random  in  nature.  Furthermore,  the  exciton  lifetime 
would  increase  monotonically  with  a  decrease  of  exciton  en¬ 
ergy. 

A  two-state  (15,  IF)  model  is  introduced  here  which 
consistently  interprets  all  the  results  shown  in  Fig.  2  as  well 
as  other  results  observed  by  Cohen  and  Sturge4  and  reported 
by  Kash  et  al .’  It  is  possible  that  the  excitons  are  confined  by 
the  potential  wells  in  the  bulk  alloy  compound  which  have 
two  lowest  states  designated  by  15  and  IF.  The  exciton  in 
higher  lying  state  lFcan  be  scattered  to  the  lower  lying  state 
15  by  emission  of  phonons  instead  of  migrating  from  site  to 
site. 

In  order  to  verify  the  two-state  model,  we  have  calculat¬ 
ed  the  eigen  energy  levels  of  the  quasi-OD  electron  systems 
and  the  recombination  lifetime  ratio  for  15-15  and  1F-1F 
transitions.  The  confinement  energy  for  the  ground  states 
(£  is  )  is  equal  to  E  +  Eh[S  for  the  quasi-OD  electron  sys¬ 
tems,  where  E's  and  EhiS  are  the  lowest  confinement  ener¬ 
gies  for  electrons  and  holes,  respectively.  These  energies 
were  calculated  from  the  bulk  band  gaps10  at  300  K  and  the 
measured  peak  energies  of  the  first  derivative  of  room-tem¬ 
perature  reflectance  ( dR/dk )  by  subtraction.  Using  the 
electron  effective  mass,  the  isotropic  hole  effective  mass  giv¬ 
en  in  Ref.  4,  and  the  measured  £,  9,  the  effective  diameter  of 
crystallite  for  each  sample  and  eigen  energies  were  calculat¬ 
ed  based  on  Eq.  ( 2 )  and  listed  in  Table  I.  The  emitted  photon 
energies  at  4.2  K  for  15-15and  1F-1F  transitions  in  Table  I 
were  obtained  by  adding  corresponding  confinement  ener¬ 
gies  to  the  values  of  bulk  gap10  at  4.2  K.  These  calculated 
energy  positions  are  located  in  Fig.  2  as  arrows  labeled  by  5 
and  F  for  clarity  and  are  in  reasonable  good  agreement  with 
the  two-state  model. 

The  radiative  lifetime  ratio  for  15-15  and  IF- IF  can  be 
predicted  following  the  approach  given  by  Casey  and  Pan- 
ish."  The  matrix  element  for  a  transition  from  localized 
state  in  conduction  band  to  the  localized  state  in  valence 
band  is  M  =  MbM„y,  where  Mb  is  the  average  matrix  ele¬ 
ment  for  the  Bloch  states  for  bands  in  absence  of  new  eigen 
states  due  to  the  confinement  and  Mmy  is  the  envelope  part 
of  the  matrix  element  which  is  Menv  =  ( 1Z.  | \L  ).  Since  the 
IF  state  is  a  threefold  degeneracy  state  with  magnetic  quan¬ 
tum  number  of  —  1,0,  and  1,  the  envelope  matrix  element 
ratio  of  \Mmy  |?P_  1P  to  |A/CT>  |}s_  1S  can  be  readily  calculat¬ 
ed  to  be  5.  When  the  4  K  recombination  lifetimes  are  purely 
radiative,  then  the  ratio  for  r,s  /r,P  should  be  equal  to  5.  The 
discrepancy  between  the  predicted  value  of  5  and  the  experi¬ 
mentally  measured  value  of  3.5  for  rls/rXP  arises  from  the 
existence  of  nonradiative  transitions.  A  simple  rate  equation 
analysis  helps  to  substantiate  this  statement.  Suppose  1/ 
T\,  =  l/rn  +  1  /rnl,  where/'  =  5,F;  r„(pl  and  r „fp)  are  the 
radiative  and  nonradiative  recombination  lifetimes  for  the 
15(F)-15(F)  transition,  respectively.  The  ratio  is 
I'ls/i’ie  =  (r„/rrp)l(r,p  +rv)/(r„  +rnj)].  When 
~r„p  and  •< r„ .r,,, ,  then  rls/rIP  - 1;  hence,  there  should  be 
only  one  lifetime.  This  is  opposite  to  the  experimental  re¬ 
sults.  On  the  other  hand,  the  ratio  rxs/rxr  is  about  3.5  when 
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TABLE  I  Measured  energy  peak  of  first  derivative  of  reflection  at  300  K,  bulk  energy  gaps  at  4  2  and  300  K,  the  calculated  confinement  energies,  emission 
energies  for  IS-1S  and  \P-\P  transitions  and  effective  diameters  Id),  the  measured  lifetimes  for  I5-I5and  \P-\P  transitions  as  well  as  their  ratio  are 
displayed 


(300  K) 

Energy 

4  K 

peak 

(300  K) 

emission 

4  K. 

Sample  for 

Bulk 

Confined  energy 

energy* 

lifetime 

CdS,  Se,  ,  dR  /dX 

g»P‘ 

(meV) 

(eV) 

(ps) 

No. 

x  d  (eV) 

(eV)  £(,+5;, 

Eu  E-u  E  f. 

F * 

15-15  \P-\P 

15  IP  r„/r„ 

(A) 

*  From  Ref  10 

"Measured  energy  separation  (meV)  between  5  and  P  peaks  at  300  K  from  Fig.  3(a). 
‘  Peaks  are  broad  and  unresolved. 

‘‘The  anticipated  values  assuming  a  same  ratio  r,,/r„  for  the  samples  2-58  and  2-64 


nonradiative  recombination  lifetimes  (r„,  =  r  )  are  about 
5  times  larger  than  the  radiative  lifetime  t,p  of  the  excited 
state  IP  for  the  assumption  of  t[P  ~t  and  r,s~r„. 

Steady-state  PL  studies  at  4  and  300  K  were  performed 
to  substantiate  our  findings  from  the  picosecond  studies.  In 
Fig.  3(a),  300  K  PL  spectra  of  sample  2-6 1  at  various  excita¬ 
tion  levels  ( see  figure  caption )  are  plotted.  The  broken  verti¬ 
cal  line  indicates  the  peak  position  of  dR  /dX.  As  expected 
the  salient  feature  in  the  luminescence  spectra  is  the  appear¬ 
ance  of  a  two-peak  structure  separated  by  220  meV  on  the 
high-energy  side  of  the  broken  vertical  line.  Based  on  the 
two-state  model  we  attribute  these  two  peaks  to  the  IS- 15 
and  1 P- 1 P  optical  transitions’  as  indicated  by  S  and  P  in  the 
figure. 

The  4  K  PL  spectrum  of  2-61  excited  by  the  457.9-nm 
line  of  the  argon-ion  laser  is  depicted  in  Fig.  3(b).  As  can  be 
seen  from  the  spectrum,  the  position  of  the  15-15  peak  shifts 
to  higher  energy  and  its  FWHM  (58  meV)  decreases  by  a 


WAVELENGTH (nm) 

FIG.  3.  (a )  Photoluminescence  spectra  of  sample  2-6 1  at  various  excitation 
levels  taken  at  room  temperature  using  the  488  nm  line  of  an  argon-ion 
laser.  The  excitation  levels  (intensity  scale)  were  2.5  (2)  and  10(10)  times 
higher  for  the  broken  curve  and  the  dot -dashed  curve  than  the  solid  curve, 
respectively.  The  broken  vertical  line  indicates  the  peak  position  of  the  first 
derivative  of  reflectance  at  room  temperature  5  and  P  stand  for  the  15-15 
and  IP-IP  transitions,  respectively  (b)  Photoluminescence  spectra  of  2-6 1 
at  4.3  K.  The  inset  shows  that  the  15  peak  broadening  is  dominated  by  the 
diameter  fluctuation  of  microcrystallites.  The  solid  dots  are  calculated  by 
assuming  a  Gaussian  shape  of  the  probability  density  function  for  the  ran¬ 
dom  variable  d  with  a  variance  of  25  A:. 
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factor  of  2  in  comparison  with  room-temperature  data.  The 
latter  indicates  that  the  broadening  of  the  15-15  transition  at 
room  temperature  is  not  entirely  dominated  by  the  size  fluc¬ 
tuation.  The  inset  of  Fig.  3(b)  shows  a  spectrum  fit  by  as¬ 
suming  a  Gaussian  shape  of  the  probability  density  function 
for  the  random  variable  d.  The  fit  gives  a  mean  value  of  74  A 
for  the  effective  diameter  and  a  full  width  at  half-maximum 
of  15  A  which  characterizes  the  size  fluctuation  of. micro- 
crystallites.  The  structure  next  to  the  main  15  peak  is  from 
the  \P~\P  transition  since  its  energy  separation  with  the  1 P 
peak  is  about  2 1 0  meV. 

We  have  also  measured  4  and  300  K  luminescence  spec¬ 
tra  of  2-59,  2-58,  and  2-64  samples.  The  two-peak  structure 
was  observed  for  2-59  at  both  4  and  300  K.  The  peaks  were 
broader  for  2-58  and  2-64  than  for  the  other  two  samples  and 
the  IP  structure  was  not  clearly  identified. 

In  summary,  we  have  reported  on  the  measurements  of 
recombination  lifetimes  of  15-15  and  IP-IP  transitions  in 
quasi-OD  electron  systems  in  CdS.Se,  _  „ .  The  assignment 
to  these  observed  transitions  is  supported  by  calculations  of 
eigen  energy  levels  and  the  corresponding  squared  matrix 
element  ratio  as  well  as  the  direct  observation  of  optical  tran¬ 
sitions  between  quantized  energy  levels  15  and  IP  in  the 
conduction  band  and  the  valence  band. 
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Abstract 


Subpicosecond  time-resolved  pnotoluminescence  studies  of  a  symmetric  multiple  quantum  well  structure  and 
an  asymmetric  quantum  well  with  a  similar  well  thickness  have  revealed  several  important  experimental 
results  on  the  behavior  of  pnotoexcited  carriers  in  these  microstructures:  (1)  The  energy  relaxation  process 
is  substantially  suppressed  m-  to  the  existence  of  a  large  population  of  ncnequilibrium  phonons  after  an 
initial  rapid  cooling.  (2)  The  two  different  well  structures  do  not  show  any  difference  in  the  energy 
relaxation  process  under  the  same  experimental  conditions.  <  3 )  The  photoexoitec  carrier  density  deduced 
from  exact  and  consistent  fittings  of  the  time-resolved  pnotoluminescence  profiles  at  various  emitted  pnoton 
energies  decreases  nonexpor.enti  al  ly  and  very  rapidly  within  the  first  30  ps  after  the  end  of  a  0.5  ps  pulse. 
An  effective  carrier  depletion  time  is  determined  to  be  as  short  as  10  ps  in  either  well.  A  mecnanism  which 
leads  to  such  short  carrier  depletion  time  is  found  to  be  associated  with  ?  he  nonequi librium  phonon  enhanced 
p  no  non  replica  emission.  Unlike  conventional  stimulated  emission  ser.avicr,  this  mechanism  is  even  more 
effective  in  exhausting  carriers  from  the  system  at  higher  temperatures.  (•*)  The  difference  in  potential 
well  profile  corresponds  to  a  very  different  behavior  cf  carrier  diffusion  process  in  the  lateral  well  plane. 
The  carrier  diffusion  is  enhanced  in  the  well  plane  of  the  asymmetric  well  by  restricting  carrier  diffusion 
in  the  growth  direction.  Diffusivity  0  of  the  pnotoexcited  carriers  in  the  asymmetric  well  has  been  directly 
determined  to  be  1 2 ‘cm Vs  at  u.JK,  which,  is  about  four  criers  of  magnitude  larger  than  the  value  in  bulk 
GaAs . 
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Over  the  last  few  years,  there  r.as  Seen  a  growing  interest  to  understand  the  fundamental  interactions  of 
confined  and  free  electrons  and  holes  with  the  phonons  m  bulk  semiconductor  and  its  multiple  quantum  well 
structures.  This  understanding  is  crucial  for  the  operatic.",  of  ultr.anigh  speed  devices  based  on  the 
principles  of  ultrafast  carrier  dynamics.  The  slowing  cf  carrier  coding  in  bulk  semiconductors  was 

attributed  to  the  screening1  of  the  electron-phonon  interaction  or  the  :■■■■- "e  at  1  rg  of  carriers  by  hot- 
phonons. i>1  The  first  comparison*  of  tr.e  hot-electron  ocolir.g  rates  m  ur. .doped  multiple  quantum  well  (MQW) 
structures  and  in  the  bulk  GaAs  was  achieved  by  time-resolved  measurements  of  optical  absorption  and  gain. 
The  cooling  rates  were  approximately  the  same  for  bulk  and  MQW  at  a  carrier  density  of  2.5xl0‘7  cm-’.  This 

is  expected  from  a  simple  theory5’  where  the  lattice  is  as  a  heat  bath  for  the  quasi-equilibrium 

carriers.  Quasi-steady  state  experiments  have  teen  carried  out  by  two  groups6' 7  to  investigate  the 

interaction  of  electrons  and  holes  with  phonons  in  modulation-doped  MQW  structures.  These  experiments  have 
generated  conflicting  results  regarding  the  presence  of  nonequi 1 ibrium  phonons.  The  first  measurements 
of  time-resolved  photoluminescence  from  modulation  doped  semiconductor  MQW  with  greater  than  20  ps  time 

resolution  were  reported  by  Ryan  et  al.“  who  found  that  the  cooling  of  hot  carriers  was  anomalously  slow 

after  R0  ps.  By  reviewing  these  and  ether  previous  works  we  noticed  three  important  points:  First,  by 
studying  the  carrier  dynamics  in  modulation-doped  MQW  structures  which  are  designed  for  high  mobility 
devices,  one  can  obtain  the  energy  loss  rates  of  the  photo-injected  majority  carriers  by  extracting  the 
carrier  temperature  as  a  function  of  time  from  the  high  energy  tails  of  time-resolved  pnotoluminescence 
spectra.  Due  to  a  large  extrinsic  carrier  density,  at  least  comparable  to  the  pnotoexcited  carrier  density, 
one  cannot  obtain  information  on  photo-injected  carrier  lifetimes.  In  these  experiments  the  carrier  lifetime 
is  essentially  close  to  the  value  in  bulk  which  is  about  1  ns. 6  It  13  necessary  to  study  photo- injected 

car-ier  dynamics  in  undoped  MQW  structures  in  order  to  obtain  information  on  both  the  energy  loss  rates  and 
carrier  lifetimes.  The  latter  is  cruci.I  for  applications  1  r:  high  speed  optoelectronic  devices.  Several 

attempts  have  been  made  to  reduce  the  carrier  lifetime  through  r.onradiativa  recombination  pr  ccesses . 10 
Second,  one  cannot  extract  information  on  the  initial  carrier  relaxation  process  from  both  quasi-steady  state 
experiments  and  time-resolved  measurements  with  time  resolution  greater  than  11  ps.  Study  of  the  initial 
carrier  relaxation  process  is  crucial  to  substantiate  the  existence  of  nonequilibrium  phonons  and  Its 
importance  to  the  hot-carripr  relaxation  since  the  phor.on  lifetime  is  about  5-7  ps.  Third,  the  rapid 
expansion  of  photoexcited  electron-hole  plasma  at  high  density  in  quantum  well  structures  has  not  been 
experimentally  investigated  on  a  picosecond  time  scale.  Only  an  indirect  measurement  of  the  expansion  of 
electron-hole  plasma  in  G3A3  MQW  structures  at  low  carrier  density  has  been  reported  where  the  drift 

velocity  of  the  electrons  was  comparable  to  the  Fermi- vel  oci  t  y,  ind  therefore,  the  transport  of  the 
photoexcited  carriers  In  GaAs  quantum  wells  was  thermodi f f uni ve . * 1 

In  this  paper,  several  experimental  results  are  reported  from  the  measurements  of  time-resolved 
pnotoluminescence  with  2  ps^time  resolution  as  well  as  t. C“  studies  of  lattice  temperature  dependence, 

excitation  dependence,  and  polarization  of  Integrated  luminescence  spectra  from  an  undoped  GaAs  symmetric 
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MQW  structure  and  an  asymmetric  quantum  well.  We  have  found:  (1)  The  energy  relaxation  process  is 

substantially  suppressed  due  to  the  existence  of  a  large  population  of  nonequilibrium  phonons  after  an 

initial  rapid  cooling.  (2)  The  two  different  well  structures  do  not  show  any  difference  in  the  energy 

relaxation  process  under  the  same  experimental  conditions.  (3)  The  photoexcited  carrier  density  -^enacted 
from  fittings  of  time-resolved  phctoluminescence  profiles  at  various  emitted  pnonon  energies  decreases  r.on- 
exponentially  and  very  rapidly  within  the  first  30  ps  after  the  (*nd  of  a  0,5  p3  laser  pulse.  An  effective 

carrier  depletion  time  is  determined  to  be  as  siic,  t  as  10  ps  in  either  type  of  wells.  A  mechanism  wnich 

leads  to  such  a  short  carrier  depletion  time  is  associated  with  nonequiiibrium  phonon  ennanced  phonon- 
replica-emission.  Unlike  the  conventional  stimulated  emission  process  this  mechanism  is  even  more  effective 

in  exhausting  carriers  from  the  system  at  higher  temperatures.  The  bandgap  renormalization  in  the 

photoexcited  region  is  not  observed  from  integrated  luminescence  spectra  as  a  consequence  of  ul 
carrier  depletion  time  and  hot  carrier  temperature.  (*0  There  is  no  or  little  carrier  diffusioi 
lateral  well  planei  of  the  symmetric  wells  except  for  the  carriers  at  the  subband  edge,  while  in 


in  the  growth  direction, 
determined  to  be  106cm2/s 
bulk  GaAs. 


Sambles 


The  undoped  GaAs/AlxGai_xAs  symmetric  MQW  structure  and  the  asymmetric  quantum  well  investigated  in 
work  were  grown  by  molecular-beam  epitaxy  cr.  a  (001  )-oriented  undoped  CaAs  substrate.  The  symmetric 
consists  of  50  periods  of  55&  thick  GaAs  and  100X  thick  Al0_  ,Ga0. 7A3  layers  followed  by  a  1.2  urn  buffer 
layer.  The  asymmetric  quantum  well  consists  of  layers  frem  the  front  surface  to  substrate  150(1  Al0  ,Ca, 
50a  GaAs  well,  200A  AlAs  (asymmetric  barrier),  followed  by  1  urn  A  10_  ,Ca^  ,As,  and  0.2  urn  CaAs  buffer"  lay 
Si-doped  GaAs  substrate.  The  lateral  sizes  are  2x4m m,  2x6mm  for  the  symmetric  MQW  structure  am 
asymmetric  well,  respectively.  The  samples  were  mounted  cn  a  cold  finger  in  a  helium  cryostat. 
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An  ultrashort  light  pulse  of  500  fs  duration  at  620  r.m  was  used  to  excite  the  electron-hole  pairs  with  a 
carrier  density  of  10‘5crrT3.  This  light  source  was  generated  from  a  colliding- pulse  passive1'  mode-locked  dye 
laser  and  amplified  by  a  four-stage  dye  amplifier  pum pea  with  a  frequency  doubled  Nd: YAG  laser  at  20  Hz. 
The  photoluminescence  was  spectrally  filtered  using  various  narrow  Dana  filters,  ana  temporally  dispersed  by 
a  2  ps  time  resolution  Hamamatsu  streak  camera  system.  The  luminescence  intensities  were  corrected  for 
nonlinearity  of  the  streak  rate,  the  spectral  response,  and  the  transmission  of  each  narrow  band  filter.  The 
time  resolution  for  the  present  work  is  about  2  ps.  The  experimental  apparatus  used  in  this  research  has 
been  described  in  detail  by  Ho  et  al.’2 


Schematic  energy  band  diagrams  are  shown  in  Fig.  1  for  a  symmetric  and  an  asymnetric  quantum  well.  Using 
a  finite-depth  square  well  model  and  material  parameters  given  in  Ref.  13  the  subband  energy  positions  in  the 
wells  wpre  calculated.  There  are  two  subbands  for  electrons  (e)  and  heavy-holes  (hh)  and  only  one  subband 
for  light-holes  (lh).  The  good  quality  of  the  samples  is  confirmed  by  spectrally  well  resolved  heavy-hole 
and  light-hole  excitonic  structures  at  room  temperature.  The  FWHMs  for  the  heavy-hole  excitons  are  about  15 
meV  and  12.7  meV  for  the  symmetric  and  asymmetric  quantum  wells,  respectively. 


The  electron-hole  pairs  were  generated  by  photo-excitation.  These  photogenerated  carriers  therm3lize 
within  themselves  via  carrier-carrier  interaction  to  roach  a  quasi-equilibrium  state  in  less  than  1  ps. 
Simultaneously,  three  other  processes  take  place.  The  first  one  is  the  loss  of  energy  of  the  electron-hole 
system  through  the  electron-phonon  interaction.  The  second  process  is  the  carrier  diffusion  process  due  to 
the  gradient  of  carrier  density  between  photoexcited  and  non-excited  areas.  The  third  process  is 
recombination  of  electron-hole  pair  through  either  nonradiative  or  radiative  mechanisms.  The  work  reported 
here  on  the  photoluminescence  experiments  uses  the  last  process  to  probe  the  first  two  processes. 


Energy  Relaxation 


a)  Steady  State  Luminescence 


Time-integrated  luminescence  spectra  of  thp  symmetric  MQW  structure  excited  by  the  subpicosecond  laser 
pulse  at  various  lattice  temperatures  (T^)  are  shown  in  Fig.  2.  These  spectra  were  taken  in  the  z(y,y)z 
configuration.  The  conventional  Raman  notation  (x,y,z  are  the  crystal  axes,  and  z  is  the  growth  direction)  is 
used.  In  order  to  eliminate  radiation  from  the  edges  of  the  sample  the  luminescence  spot  wa3  imaged  on  an 
aperture  with  the  aid  of  a  3treak  camera  and  then  focused  onto  a  vertical  slit  of  a  grating  spectrometer. 
Several  features  are  displayed  in  the  data. 


(1)  The  emission  peak  (A)  on  the  high  energy  stdp  of  the  spectra  is  from  the  recrmbinatlon  of 
photogenerated  n-1  electrons  and  n-1  heavy-holes.  This  peak  shirts  towards  the  low  energy  side  as  TL 
Increases.  A  high  energy  tail  on  the  main  peak  (A)  develops  with  the  increase  of  TL.  For  TL>'50K,  a 
shoulder  on  the  high  energy  side  of  the  peak  A  arising  from  the  recombination  of  n-1  electrons  and  n»i  light- 
holes  can  be  identified. 
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(2)  A  broad  emission  band  (C)  at  the  low  energy  side  of  the  spectra  arises  from  the  1.2  um  GaAs  buffer 
layer.  The  total  emission  intensity  decreases  as  increases. 

(3)  The  most  interesting  feature  of  the  data  is  the  appearance  of  an  emission  band  (3)  below  the  n-1 
electron  to  n-1  heavy-hole  transition  which  is  attributed  to  the  emission  from  the  nonequilibrium  phonon 
enhanced  phonon  replica. 

The  following  reasons  support  this  assignment: 

(i)  The  B  emission  band  does  not  appear  in  the  steady-state  photoluminescence  spectra  taken  at  low 
power  excitation  of  about  1W/cm2  at  4.3K  using  a  488  nm  line  of  an  argon-ion  laser.  However,  the 
photoluminescence  spectra  using  a  weak  train  of  laser  pulses  (-  200  fs)  directly  from  the  colliding-pulse 
passively  mode-locked  (CPM)  dye  laser  oscillator  as  the  excitation  source,  with  an  excitation  power  density 
in  the  range  of  1 0"5-1  0"3W/cm2,  show  a  weak  electron-acceptor  emission  band  separated  by  17  meV  from  the  n-1 
electron-heavy-hole  transition  accompanied  by  its  phonon-replica1'  at  low  temperatures.  This  extrinsic 
emission  band  disappeared  completely  when  the  sample  temperature  was  raised  to  -  80K,  whereas  the  3  band  in 
Fig.  2  exists  up  to  room  temperature.  Moreover,  the  relative  time-integrated  intensity  ratio  between  B  and  A 
bands  in  Fig.  2  detected  under  amplified  0.5  ps  light  pulse  excitation  decreases  as  the  excitation  power 
density  decreases  and  the  B  band  is  not  visible  when  the  excitation  power  density  is  lower  than  10"’Pm,  where 
Pm  -  1012W/cm2  is  the  maximum  value  of  the  excitation  power  density.  This  contradicts  what  is  generally 
expected  for  an  impurity  emission.  Since  the  concentration  of  acceptors  is  low  in  our  sample  as  confirmed  by 
the  luminescence  studies  using  the  weak  train  of  200  fs  light  pulses,  the  impurity  emission  should  be  readily 
apparent  and  more  pronounced  at  lower  excitation.  Therefore,  the  B  emission  band  shown  in  Fig.  2  cannot  be 
attributed  to  electron  to  acceptor  luminescence. 

(ii)  The  ratio  of  emission  intensities  of  B  and  A  emission  bands  increases  as  T^  increases.  The  values 
of  the  ratio  as  a  function  of  T^  are  shown  in  Fig.  3  as  triangles.  This  behavior  is  remarkably  different 
from  the  conventional  stimulated  emission  process  in  which  phonons  do  not  participate.  As  in  semiconductor 
lasers,  the  stimulated  emission  is  less  effective  when  the  sample  is  at  room  temperature.  To  convince  us 
that  the  emission  band  B  is  due  to  a  phonon-assisted  process,  we  calculated  the  sum  of  occupation  numbers 
for  both  equilibrium  phonons  (lattice  temperature  T[J  and  nonequilibrium  phonons  (carrier  temperature  Tc)  by 
assuming  the  carrier  temperature  is  the  same  as  the  effective  temperature  for  the  nonequilibrium  phonons 
after  30  p3  (see  later  discussion  on  the  carrier  temperature).  The  calculated  result  is  shown  by  the  solid 
line  in  Fig.  3.  The  fitting  of  the  total  phonon  occupation  number  to  the  obtained  intensity  ratio  from  the 
spectra  is  impressively  good.  This  implies  that  the  intensity  of  peak  B  is  well  correlated  with  the  L0- 
phonon  population.  Because  of  participation  of  nonequilibrium  phonons  the  intensity  of  the  B  band  (see  also 
time-resolved  luminescence  at  780  nm  in  Fig.  4)  is  strongly  enhanced  especially  at  low  temperatures.  This 
further  supports  our  assignment  of  the  B  band  to  the  nonequilibrium  phonon  enhanced  phonon  replica. 


(ill)  It  is  not  possible  to  attribute  the  emission  band  B  to  the  emission  from  the  renormalized  band-band 
transition  because  the  peak  positions  of  the  A  band  (766. 7nm)  and  B  band  as  well  as  the  energy  separation 
between  them  at  4.3K  do  not  change  with  variation  of  the  excitation  intensity  by  a  factor  of  9x10"'  (2ND3). 
Further  support  to  the  above  statement  is  the  fact  that  the  spectral  position  of  peak-A  exactly  coincides 
with  each  other  for  two  different  luminescence  studies  using  different  light  excitation  sources;  the  one  is 
the  weak  200  ts-pulse  train  with  a  repetition  rate  of  114  MHz,  while  the  other  is  the  amplified  0.5  ps-pulse 
with  a  20Hz  repetition  rate.  As  the  observed  band  edge  transition  energy  does  not  change  over  a  range  of 
photogenerated  carrier  densities  (101,-101 5/cmJ)  the  bandgap  renormalization  in  this  study  must  be  excluded. 

(Iv)  We  have  also  measured  the  polarization  of  the  emission  bands  A,  B,  and  C  by  adopting  a  "right- 
angle"ls>“  geometry  for  obtaining  integrated  luminescence  spectra.  It  was  found  that  the  emission  of  band  B 
detected  along  the  y-direction  (emitted  at  the  sample  edge)  was  strongly  polarized  in  the  x-direction,  while 
the  A  and  C  emission  bands  are  depolarized  independent  of  excitation  power  density  and  lattice  temperature. 
The  intensity  ratio  between  the  x-  and  z-polarization  for  the  B  band  was  about  20  and  independent  of  lattice 
temperatures  from  4.3  to  300K  at  full  excitation  power  density  Pm.  But  it  is  strongly  dependent  on  the 
excitation  power  density  at  a  given  lattice  temperature.  For  example  at  TL-100K,  the  B  band  cannot  be 
clearly  identified  from  the  A  band  and  shows  no  difference  between  the  two  polarizations  at  the  excitation 
power  density  0.017  Pm  (2ND13),  while  the  ratio  at  full  excitation  Pm  is  20.  The  polarization  behavior  of  the 
B  band  is  consistent  with  the  Raman  scattering  measurements  reported  by  Zucker  et  al.11  where  the  initial 
photons  are  provided  by  the  external  laser  source.  In  our  case,  the  initial  photons  are  provided  by  the 
luminescence  from  the  recombination  of  n-1  electrons  and  n-1  heavy-holes  at  the  subband  edges.  This  process 
is  3trongly  enhanced  by  the  presence  of  a  large  number  of  nonequilibrium  phonons  emitted  by  the  hot- 
electrons  at  high  excitation.  The  polarization  behavior  of  emission  band  C  is  expected  since  the  emission  is 
from  the  1.2pm  GaAs  buffer  layer.  This,  in  fact,  provides  a  good  check  of  our  system.  However,  the 
polarization  behavior  of  the  emission  band  A  is  not  understood  in  terms  of  the  selection  rules  of  dipole 
recombination  described  by  Iwamnra  et  al.1T,  where  the  emission  intensity  for  x-polarization  (TE  polarization) 
should  be  at  least  four  times  larger  than  that  for  z-polarization  (TM  polarization)  due  to  different  density 
of  states  for  heavy-holes  and  light-hole3. 

There  has  been  much  debate1’"1’  on  the  interpretation  of  the  spectral  feature  below  the  n-1  electron 
heavy-hole  transition  energy.  This  is  because  the  energy  levels  for  Impurity  states,  longitudinal  optical 
phonon  (LO)-phonon  replica,  and  Intrinsic  bandgap-reduction  due  to  many-body  interactions  roughly  coincide 
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uith  each  other.  Holonyak  and  co-workers1*  have  demonstrated  the  LO-phonon  participation  In  GaAs-AlGaAs 
based  QW  laser  emission  by  the  observation  of  more  than  one  LO-phonon  sideband  of  laser  operation  below  the 
n*1  confined-particle  transition.  Also,  the  phonon  replica  of  n*l  electron  to  heavy-hole  emission  has  been 
observed  in  the  InGaAs  based  systems.  The  observation  of  the  coupling  at  both  the  GaAs-  and  Ir.As-like  bulk 
LO-phonon  energies  of  InGaAs  QW  by  Skolnick  et  al.2°  has  confirmed  that  the  peak  below  the  n  =  l  electron- 
heavy-hole  transition  is  a  phonon  replica.  Furthermore,  it  should  be  pointed  out  that  our  assignment  for  the 
B  emission  band  is  consistent  with  the  time-resolved  pnotoluminescence  data  described  in  the  next  section. 

b)  Time-Resolved  Luminescence 

Time-resolved  photoluminescence  emitted  from  the  symmetric  MQW  structure  at  n.  3K  are  shown  in  Fig.  U  for 
various  emitted  photon  energies.  The  emission  centered  at  770  nm  arises  from  the  recombination  of  electrons 
and  heavy-holes  near  the  respective  first  subband-edges.  The  emission  with  wavelengths  less  than  770  nm  is 
from  the  recombinations  of  hot  carriers.  The  time-resolved  luminescence  intensity  assigned  to  the  phonon 
replica  3t  780  nm  is  shown  in  Fig.  k.  Each  luminescence  profile  was  the  average  of  10  individual  shots  using 
a  prepulse21  for  averaging.  The  left  peak  of  the  dotted  curve  is  the  prepulse  which  reflects  the  2  ps 
temporal  resolution.  The  right  peak  on  the  same  curve  is  the  Rayleigh  scattering  light  from  the  sample 
surface  which  defines  the  "zero"  time  for  our  analysis.  Several  features  appear  in  the  data  displayed  in  Fig. 

(i)  Although  the  risetime  for  all  the  luminescence  profiles  is  instrumental,  the  2  ps  uppei — limit 
reflects  the  rapid  thermalization  and  initial  cooling.  (ii)  The  shape  of  the  rise-time  of  the  luminescence 
profile  at  780  nm  is  similar  to  that  at  770  nm,  but  delayed  by  3  ps.  This  implies  that  the  emission  at  780 
nm  does  not  originate  from  the  same  band  as  the  emission  at  770  nm;  otherwise,  the  rise-time  of  the 
luminescence  at  780  should  start  at  the  same  zero  point  as  that  at  770  nm.  The  3  ps  delay  is  consistent 
with  our  assignment  of  the  emission  band  center  at  790  nm  to  the  phonon  band.  This  time  is  just  the 
effective  time  required  to  establish  a  large  nonequilibrium  hot  phonon  population.  (iii)  The  luminescence 
decay  time  of  30  ps  (broken  curve)  for  the  phonon  replica  emission  at  780  nm  cannot  be  due  to  impurity 
emission  at  high  excitation. 

The  decay  times  were  obtained  by  fitting  the  data  in  Fig.  k  to  the  expression  given  by: 

"t/Td 


I(E,t)-I0[  1  Tp  je 


(!) 


where  the  parameters  (I0,tr  Td^  are  the  proportional  constant,  the  risetime,  and  the  decay  time,  respectively. 

To  study  the  energy  relaxation  process  quantitatively,  the  carrier  temperature  as  a  function  of  time  must 
be  determined.  There  are  many  ways  to  obtain  the  time  evolution  of  carrier  temperature.  One  way  is  to 
analyze  the  high  energy  tail  of  time-resolved  spectra  using  a  Maxwell-Boltzmann  distribution.6-*  Another  way 
is  to  measure  the  ratio  of  luminescence  intensity  profiles  as  a  function  or  time  at  two  energy  positions 
within  the  exponential  energy  tail  by  also  assuming  a  Maxwell-Boltzmann  distribution.  The  above  two  methods 
cannot  extract  information  on  the  photogenerated  carrier  density.  An  expression  in  the  time  domain  is 
introduced  to  fit  the  experimental  luminescence  temporal  profiles  by  using  two  adjustable  parameters, 
namely,  the  carrier  densities  (ne,  n^p,  n^p)  and  the  electron  temperature  (Tc). 

Assuming  the  direct  optical  transition,  the  luminescence  intensity  is  given  by: 


I(Ei,t)=Ci(1  e  ^U^e-hh  l*Pe  ^e  Phh  ^h  *  l^e-lhlJ(>e  ^e  Plh  ^h^* 


(2) 


where  pe  hh  -  me  lh  is  the  density  of  state  for  electrons,  heavy-holes,  and  light-holes;  Mp-hh(lh)  is 
the  matrix  'element  'for2  electron  to  hh(lh)  transitions;  Cj  absorbs  all  the  constant  factors  including  the 
corrections  for  detector  response  and  the  transmission  of  each  narrow  band  filter  used;  tr  is  the  risetime  of 
luminescence  and  is  set  to  be  1  ps  (approximately  the  response  time  of  the  system)  which  is  an  upper-limit 
of  thermalization  time  of  the  electron-hole  system; 


1  /(etci"ue,h 


)/KTr 


1) 


(3) 


is  the  Fermi-Dirac  distribution  for  electrons  in  the  conduction  band  (with  subscript  e)  and  for  holes  in  the 
valence  band  (with  subscript  h).  The  quasi-Fermi  energy  of  holes  (uh)  is  related  to  the  quasi-Fermi  energy 
of  electrons  (u^)  and  the  carrier  temperature  by  the  relation. 


ne’nhh  *  nlh 


This  value  13  given  by: 


with 


uh  >  KTclnCe31^1  ♦eup/KTc)-l  ] 


w 


(5) 


mhh  6  mlh  (6> 

An  unique  set  of  parameters  of  Tc(t)  and  ue(t)  was  used  to  consistently  fit  all  the  luminescence  profiles 
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detected  at  different  photon  energies  Ej.  Three  typical  calculated  luminescence  profiles  are  shown  in  Fig.  <4 
by  the  thick  solid  curves. 


In  Fig.  5  the  time-resolved  photoluminescence  profiles  from  the  asymmetric  quantum  well  are  shown  at 
three  typical  energies.  The  emission  centered  at  760  nn  arises  from  the  recombinations  of  electrons  and 
heavy-holes  near  the  respective  first  subband-edges.  The  theoretical  solid  curves  in  Fig.  5  were  generated 
using  Eq.(2).  The  same  set  of  values  for  Tc  and  up  used  to  fit  the  data  for  the  symmetric  quantum  well  were 
used  to  fit  the  luminescence  prefiles  of  the  asymmetric  quantum  well.  The  energy  relaxation  process  is 
similar  for  both  wells  despite  the  different  shape  of  the  well.  This  is  as  expected! 


The  experimentally  determined  v; 
The  shaded  area  reflects  the  extent 
due  to  our  limited  time  rescluti 
collisions  quickly  lead  to  a  therm 
excited  carriers.22  The  therma 
spectroscopy.  The  Initial  cooli 
resolution  of  2  pa,  providing  in 
lifetime.  The  initial  cooling 
distribution  for  electrons.  After 
finite  wave  vector  space  ar»  uccur.u 
electrons  may  be  reabsorbed  by  Vr 
cool  \r.s  for  * h**  r.o T.-C3  priors.  It 
importance  of  significant  screening 
screening  was  important,  the  mitia 
not  te  built  up.  Another  important 
slow  component  of  the  carrier  t-m: 
phonon  enhanced  phonon  replica  er.i 
and  tre  nonequi librium  phor.or,  syst' 
a  single  decay  constant  of  37  ps  a: 


aiues  of  Tc  as  a  function  of  time  are  plotted  as  the  solid  curve  in  Fig.  6. 
t  of  uncertainty  in  deducing  the  cannier  temperature  within  the  first  H  ps 
ion.  The  plot  of  Tc  vs  t  is  interpreted  as  follows:  carrier-carrier 

alized  distribution  at  very  nigh  Tc  within  1  ps  due  to  a  large  number  of 
lization  process  of  carriers  can  only  be  probed  by  femtosecond 

ng  of  the  thermalized  distribution  is  studied  with  our  present  time 
(formation  about  both  the  electron-phonon  interaction  and  the  carrier 
within  tne  first  ‘3  ps  is  very  fast  (25CK/ps)  assuming  a  Fermi-Dirac 
r  the  initial  cooling,  a  large  number  of  longitudinal  optical  phonons  in 
il  iteq  u-‘  to  the  finite  pnonon  lifetime.23  The  phonons  emitted  by  the  hot 
he  electrons  as  the  reverse  process  of  emission  giving  rise  to  a  slower 
sr.onl  5  b"  -mpr.  is: zee  r.»r«  tr.at  the  rapid  initial  cooling  rules  out  the 
•:  of  the  tror-nr.onor.  interaction  in  the  present  study.  If  the  initial 
*1  ’"in?  ceding  would  be  slow  and  a  nonequilibrium  phonon  population  would 
i  otser.Mt; -r.  displayed  m  Fig.  o  is  that  the  time  constant  -  30  ps  for  the 

per:' are  cooling  curve  is  the  same  as  the  decay  time  of  the  nonequilibrium 

ission  at  TOO  nm.  This  suggests  that  tne  time  required  for  the  electron 
e-.s  to  “quilitrate  with  ••  a cn  otr.er  is  *»-3  PS.  The  system  then  decays  with 
3  a  coupled  hot  electron-pnonon  system. 


Th“  experimentally  determined  quati-Fermi-erergies  for  electrons  and  holes  are  plotted  in  Fig.  7.  The 
changes  of  up  and  u^  are  very  rapid  within  the  first  10  ps  and  the  behavior  of  the  degeneracy  for  electrons 
and  holes  as  a  function  of  time  are  reversed. 

Using  the  experimentally  determined  distributions  for  electrons  and  holes  in  terms  of  Tc(t)  and  up  y/t) 
the  energy  loss  rates  for  electrons  and  holes  as  function  of  Tc  or  np  can  be  obtained.  The  energy  loss  rate 
(?*,r.n*  is  defined  as  follows: 


d<E>„  rr 


C  f  l  Tp  , Ur 


f  (Tp.up^yifj.C/dC 


However,  if  the  electron-hole  system  is  treated  as  a  Maxwell -Boltzmann  gas,  i.e.,  the  electrons  and  holes 
have  the  same  distribution  function  which  is  a  Maxwell-Boltzmann  distribution,  the  energy  loss  rate  for  the 
carriers  (c-e,hh)  is  given  by: 


It  should  be  pointed  out  that  the  energy  loss  rates  for  electrons  and  heavy-holes  obtained  by  Eq.  (7)  are 
the  net  loss  rate3.  It  includes  the  energy  loss  due  to  electron(hoie)-phonon  interactions  and  the  energy 
exchange  between  electrons  and  holes  by  carriei — carrier  scattering.  The  hole-phonon  scattering  rates  were 
calculated2"  and  measured1  to  be  2.5  to  3  times  larger  than  the  rate  for  electrons  due  to  the  additional 
coupling  through  the  deformation-potential.  Therefore,  an  energy  transfer  process  from  electrons  to  holes 
should  be  expected  in  our  photogenerated  carrier  system.  The  experimentally  determined  energy  loss  rates 
using  Eq.  (7)  or  Eq.  (8)  are  shown  in  Fig.  8.  The  solid  and  dotted  curves  were  obtained  from  Eq.  (7)  for 
electrons  and  heavy-holes,  respectively.  The  dot-dashed  curve  was  obtained  from  Eq.  (8).  It  should  be 
pointed  out  that  the  curves  obtained  from  either  Eq.  (7)  or  Eq.  (8)  are  the  experimental  data  because  the 
distribution  f(Tc,uPih)  was  experimentally  determined.  The  broken  line  is  calculated  based  on  the  simple 
theory1*  where  lattice  can  be  treated  a3  a  heat  bath  for  quasi-equilibrium  carriers.  This  theoretical  curve 
does  not  agree  with  experimentally  determined  curves.  The  smaller  Phh  and  Pp  is  because  of  the  presence  of 
light-hole  population  and  the  different  masses  for  electron  and  hh  which  lead  to  different  behaviors  of 
qua3i-Fermi  energies  for  electrons  and  heavy  hole3  as  a  function  of  time. 
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The  salient  features  of  the  experimental  data  displayed  in  Fig.  8  are: 


i 

£ 
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(i)  The  experimentally  obtained  energy  loss  rate  for  electrons  is  more  than  two  orders  smaller  than 
predicted  by  the  simple  theory5  (broken  line).  The  reasons  for  the  differences  arise  from  the  reabsorption 
of  nonequilibrium  phonons  by  electrons.  As  stated  above,  the  energy  loss  rate  determined  for  electrons  by 
scattering  alone  with  LO-phonons  via  the  wave-vector  dependent  Frohiich  interaction  may  be  even  smaller  than 
our  determination  because  the  net  energy  loss  rate  determined  from  Eq.  (7)  is  the  sum  of  rates  caused  by 
electron-phonon  and  electron-hole  interactions. 

(ii)  The  net  energy  loss  rates  for  electrons  and  heavy-holes  almost  follow  the  rate  obtained  from  Eq.  (8). 
The  reasons  for  this  are  the  rapid  decrease  of  Up  as  shown  in  Fig.  (8)  and  the  slow  cooling  of  the  carrier 
temperature  which  results  in  ue/KTc  <  0.5.  This  makes  the  electron-hole  system  behave  more  like  a  Boltzmann 
gas. 


£ 


(iii)  The  energy  loss  rates  increase  as  Tc  increase.  The  dot-dashed  curve  even  bends  over  and  approaches 
the  maximum  value  of  5.3xlO"‘W  for  bulk  GaAs.  This  "bending"  behavior  reflects  the  initial  rapid  cooling  of 
hot  carriers  and  was  also  observed  by  Shah  et  al.J  This  may  be  the  indication  for  the  presence  of 
nonequilibrium  phonons  in  the  system. 
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In  Fig.  9,  we  have  plotted  the  average  phonon  emission  time  for  hot  carriers  defined  by:’ 

_  .  -Elo'KTc 

Tavg  pe 

as  a  function  of  the  carrier  temperature.  This  is  just  another  way  to  describe  the  energy  relaxation 

process.  As  can  be  seen  the  tavK  is  carrier  temperature  dependent  and  therefore  time  dependent.  For 

TC>1  200K  or  t<5  ps  the  Tav^  remains  at  a  value  of  about  !  ps.  As  time  t>5  ps  (or  Tc<12CCK)  the  taVg 
increases  quickly  with  decrease  of  the  carrier  temperature.  At  Tc  =  980K,  i2Vg  is  10  P3-  Ryan  et  al . 
obtained  a  constant  value  of  7  ps  for  taVg  by  studies  of  time-resolved  photoluminescence  of  modulation 
doped-uQW.  The  reason  for  the  difference  is  simply  pecause  their  time  resolution  did  not  allow  them  to 
monitor  the  entire  carrier  cooling  process.  Gur  upper  limit  of  G.9  ps  for  the  initial  carrier  cooling  agrees 
with  the  theoretically  predicted’  value  of  0.16  ps  within  our  resolution.  As  soon  as  a  large 

nonequilibrium  phonon  population  is  built  up  the  carrier  cooling  is  suppressed. 


We  have  extensively  discussed  the  energy  relaxation  process  based  on  the  experimentally  determined 
distribution  function  f  e  h(Tc,ue  ,  c ) .  In  the  same  manner  information  about  the  carrier  lifetime  can  be  also 
obtained  using  determined  fe  '  In  Fig.  10  th?  carrier  densities  for  electrons,  heavy-holes,  and  light-roles 
as  a  function  of  time  are'  determined  from  the  experimental  data.  The  densities  are  obtained  by  the 
following  expressions: 


(10) 
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where  AE  is  the  energy  separation  between  the  hh-  and  lh-subtands  at  the  zone  center.  The  relation  ne  -  n 
♦  njh  must  be  satisfied.  The  quantity  of  In  [ne(t)j  is  also  plotted  in  Fig.  10  as  a  broken  curve. 


hh 


The  salient  feature  of  the  data  in  Fig.  10  is  that  the  carrier  density  decreases  nonexponentially  and  very 
rapidly  within  the  first  30  ps  after  the  end  of  a  0.5  ps  pulse  excitation.  One  cannot  define  a  single  carrier 
lifetime  due  to  the  nonexponential  nature  of  the  density  decay  curves.  However,  an  effective  carrier 
depletion  time  (time  for  the  density  to  decrease  by  a  factor  of  e-'  from  ne(t-0))  can  be  deduced  to  be  as 
3hort  as  10  ps.  This  decrease  of  carrier  density  occurring  in  such  a  short  time  cannot  be  accounted  for  by 
the  usual  blmolecular  recombination,  which  is  a  much  slower  process,  on  the  nanosecond  time  scale  in  the 
bulk  GaAs  and  about  350  ps  In  quantum  well  structure  with  well  thickness  of  -  5  nm.”  This  short  carrier 
depletion  time  is  not  caused  by  the  rapid  carrier  diffusion  which  will  be  discussed  in  the  diffusion  section. 

Such  a  short  carrier  depletion  time  can  be  explained  consistently  by  the  participation  of  nonequilibrium 
phonons.  The  nonequilibrium  L0  phonons  In  the  quantum  well  system  can  play  two  essential  roles.  First, 
they  substantially  suppress  the  en  rgy  relaxation  process  of  carriers  due  to  reabsorption  of  these  phonons 
by  carriers.  Second,  since  these  phonons  are  localized  in  both  wavevt"tor  space  and  the  real  physical  space 
(In  the  well  plane)  they  will  behave  like  coherent  Bosons.  The  existence  of  the  coherent  Bosons  will  tend  to 
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increase  the  number  of  the  bosons  in  the  system  with  a  rate  proportional  to  the  present  number  of  Bosons 
causing  a  rapid  carrier  relaxation.  This  is  a  stimulated  phonon  emission  process  and  can  only  be 
accomplished  by  the  phonon  replica  emission,  i.e.,  when  an  electron  in  the  conduction  subband  recombines  with 
a  hole  in  the  valence  subband  not  only  a  photon  but  also  a  phonon  or  more  phonons  will  be  emitted  to  Join 
the  Boson  system.  Since  this  is  a  stimulated  process,  it  is  very  effective  in  exhausting  carriers  from  the 
system.  Unlike  the  conventional  stimulated  emission  in  a  semiconductor  laser,  this  process  is  expected  to  be 
even  more  effective  at  room  temperature  since  equilibrium  phonon  occupation  will  also  increase  as  lattice 
temperature  increases.  We  have  confirmed  the  nonequilibrium  phonon  induced  stimulated-phonon-emission  by  the 
detailed  polarization  studies  of  integrated  luminescence  spectra  discussed  in  the  last  section.  By  employing 
the  nonequilibrium-phonon  model  we  are  able  to  explain  not  only  the  slow  carrier  cooling,  but  also  the 
ultrashort  carrier  depletion  time  in  highly  photoexcited  undoped  quantum  well  structures. 

Carrier  Diffusion 

Electron-hole  plasma  expansion  at  high  photoexcitation  power  density  in  semiconductor  bulk  and  quantum 
well  structures  is  a  subject  of  great  interest  from  both  the  points  of  view  of  theory  and  device  application. 
Various  methods  have  been  employed  to  study  the  dynamic  process  of  carrier  diffusion  such  as  spatially  and 
spectrally  time-resolved  luminescence  spectroscopy20-27  and  Raman  scattering. “i2*  Different  mechanisms  have 
been  suggested  to  explain  the  rapid  carrier  diffusion  in  bulk  semiconductors.  Rapid  electron-hole  diffusion 
has  been  observed27  in  CdSe  at  RK  but  not  at  room  temperature.  The  absence  of  a  larger  number  of  thermal 
phonons  at  low  temperature  made  it  possible  to  observe  the  rapid  carrier  diffusion  process.  Using  a  streak 
camera,  Junnarkar  and  fllfano2’  were  able  to  measure  the  difference  of  spatial  profiles  of  the  pump  laser  and 
the  photoluminescence  at  various  excitation  levels  by  which  the  diffusion  velocity  was  estimated.  Junnarkar 
and  Alfano2’  obtained  a  diffusion  velocity  on  the  order  of  10‘cm/s  at  high  photoexcitation  even  at  room 
temperature.  Such  a  large  diffusion  velocity  was  attributed  to  screening  of  polar  electron-phonon 
interactions  and  partially  screened  nonpolar  electron-phonon  interactions.  High  internal  pressure  of  a  very 
dense  electron-hole  plasma  has  been  proposed20  to  account  for  the  rapid  expansion  of  plasma.  The  pioneering 
work  of  Tsen  and  Morkoc11  using  the  time-resolved  Raman  scattering  technique  has  shown  no  Influence  on 
photogenerated  carrier  diffusion  by  carrier  confinement,  zone  folding  of  phonon  modes,  and  phonon  confinement 
in  symmetric  GaAs/AlGaAs  MQW  structure  at  low  carrier  density.  In  view  of  their  works  we  should  point  out 
that  the  well  thickness  of  their  MQW  structure  was  so  thick  that  any  confinement  effects  would  be  smeared 
out;  and  the  carrier  density  deduced  from  Raman  data  was  the  total  carrier  density  within  whole  energy 
distribution.  Therefore,  they  could  not  obtain  the  information  about  dynamic  diffusion  of  carriers  with 
different  kinetic  energies.  In  the  following  subsection  the  results  on  dynamics  of  carrier  diffusion 
processes  in  the  x-y  plane  (lateral)  will  be  presented. 

A  schematic  diagram  of  the  photoluminescence  image  on  the  screen  of  the  streak  camera  is  depicted  in  Pig. 
1  1  to  show  how  we  study  dynamics  of  carrier  diffusion.  The  area  inside  the  filled-llne  circle  is  the  laser 
excited  area.  If  there  is  no  carrier  diffusion  photoluminescence  should  be  restricted  to  inside  this  circle. 
Suppose  sufficient  electrons  and  holes  can  diffuse  far  enough  to  the  unexcited  area  where  they  can  recombine 
one  should  be  able  to  detect  the  luminescence  from  the  unexcited  area  (outside  the  filled-line  circle). 
Using  the  streak  camera  we  were  able  to  measure  the  laser  spatial  profile  and  luminescence  profile  at  the 
sample  surface.  Those  spat.ul  intensity  profiles  were  integrated  over  a  time  duration  of  1  ns.  The  width 
(FWHM)  of  spatial  profile  should  give  information  on  the  carrier  diffusion  length  within  1  ns.  However,  this 
information  is  not  adequate  to  deduce  accurate  values  of  diffusion  constant  and  diffusion  velocity.  To  study 
the  dynamic  diffusion  process  quantitatively  we  make  use  of  two  electronic  windows  to  obtain  time-resolved 
photoluminescence  from  the  excited  (I3C)  and  unexcited  (Ise)  areas  at  the  same  time  (and  thus  laser  intensity 
fluctuation  was  completely  avoided). 

The  time-integrated  laser  and  luminescence  profiles  at  the  surface  of  the  55 X  symmetric  MQW  structure  are 
plotted  in  Fig.  12.  The  dotted-curve  is  the  laser  profile.  The  inside  and  outside  solid  curves  are  the 
luminescence  profiles  at  7A0  nm  and  770  nm,  respectively.  The  broken-curve  is  the  luminescence  profile  at 
830  nm,  which  comes  from  the  1.2  urn  buffer  layer.  The  gain  of  the  detecting  system  was  adjusted  such  that 
the  maximum  Intensity  recorded  was  roughly  the  same.  The  widths  of  luminescence  profiles  above  the  subband 
edge  (760-700  nm)  are  exactly  the  same  as  the  width  of  the  laser  profile  (dotted- line)  reflecting  no  or 
little  diffusion  for  high  energy  carriers  within  our  spatial  resolution.  The  width  of  the  luminescence 
profile  at  the  subband  edge  (770  nm)  is  about  2  times  wider  than  the  width  of  the  laser  profile.  This  energy 
dependence  of  the  carrier  diffusion  process  is  not  usually  expected  because  a  larger  kinetic  energy  is 
associated  with  larger  carrier  velocities.  However,  the  nonequilibrium  phonons  in  the  system  during  the 
first  30  ps  after  the  laser  excitation  may  inhibit  carrier  diffusion,  especially  for  the  carriers  with  high 
energy.  The  extent  of  carrier  diffusion  at  the  subband  edge  seems  to  be  consistent  with  the  results  given  by 
Tsen  and  Morkoc.11  The  diffusion  velocity  determined  by  them11  is  8.6x10‘cm/s.  According  to  this  value  the 
carriers  will  travel  about  100  pm  within  1  ns  which  will  broaden  the  width  of  luminescence  profile  from  200 
tun  to  about  A00  pm.  This  agrees  with  the  measured  lateral  width  or  the  luminescence  profile  at  770  nm  in 
Fig.  12.  It  should  be  pointed  out  that  the  initial  carrier  density  in  our  experiment  was  about  100  times 
larger  than  in  their  experiment.  Due  to  the  short  carrier  depletion  time  (10  ps)  in  our  case,  the  average 
carrier  density  within  1  ns  would  be  similar  in  both  cases. 

The  time-integrated  laser  and  luminescence  profiles  at  the  surface  or  the  asymmetric  quantum  well  are 
plotted  in  Fig.  1  3.  The  Innermost  solid  curve  and  dotted  curve  are  the  laser  profile  and  the  Caussian-line- 
flt  of  the  laser  profile,  respectively.  The  center  wavelengths  of  the  luminescence  are  indicated  on  the 
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curves.  The  key  differences  of  the  spatial  luminescence  intensity  profiles  of  the  asymmetric  well  from  that 
of  the  symmetric  wells  are:  the  widths  of  luminescence  profiles  are  much  wider  than  the  width  of  the  laser 
profile;  and  the  width  of  the  luminescence  profiles  do  not  depend  much  on  the  wavelength  of  luminescence. 
The  possiole  mechanisms  are  given  below. 

In  Fig.  1H,  we  have  plotted  the  time-resolved  photoluminescence  from  I3e  and  Isc  measured  at  4K  with  2  ps 
time  resolution  for  asymmetric  well.  Based  on  the  fact  that  the  width  of  the  spatial  luminescence  profile 
does  not  depend  on  the  photon  energy  detected  within  our  spatial  resolution,  w»  placed  a  wide-band  filter 
centered  at  700  nm  in  the  luminescence  path  to  increase  the  strength  of  the  luminescence  intensity.  The 
broken  curve  is  the  Rayleigh  scattering  light  from  the  sample  surface.  As  can  be  seen,  the  risetime  of  the 
Isc  is  time-resolution  limited,  while  the  rise  of  I3e  is  much  slower  than  that  of  Is;;  3nd  also  well  resolved. 
The  fast  decay  of  Iso  has  been  discussed  in  the  energy  relaxation  section.  It  should  Be  emphasized  here  that 
using  the  same  set  of  experimentally  determined  Tc(t)  and  ue(t)  we  have  determined  the  carrier  depletion  time 
is  10  ps.  The  only  choice  for  theoretical  fittings  to  Ise  is  the  use  of  standard  classical  diffusion  equation 
to  calculate  the  carrier  density  as  a  function  of  time  t  and  distance  r.  The  exact  values  of  separation 
between  the  two  windows  and  the  widths  of  the  windows  are  given  in  Fig.  II.  This  method  should  be  adequate 
as  long  as  electrons  act  more  like  particles  rather  than  waves. 

The  diffusion  equation  is  given  by: 
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where  t  is  the  carrier  lifetime,  and  D  is  the  diffusion  constant  of  carriers.  The  initial  condition  for 
carrier  density  is  assumed27  to  be  the  same  as  the  laser  profile  which  is: 
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where  n0  is  the  carrier  density  at  the  center  of  the  laser  spot  at  t=0,  and  a  is  1  2C  pm.  obtained  by  fitting 
the  laser  profile  to  Eq.  (14).  The  solution  of  the  second  order  differential  Eq .  (13)  with  the  given  initial 
condition  Eq.  (14)  is  given  by:27 
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Assuming  ambipolar  diffusion,  i.e.  electrons  and  holes  have  tn<  same  diffus  on  velocities, 
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Using  Eq.  ( 1 6)  the  parameters  or  D  and  t  have  teen  determined  to  be  li‘cs2/s  and  12  ps,  respectively. 

The  results  are  very  surprising  in  two  respects.  First,  the  value  fcr  'J  is  more  than  four  orders  of 
magnitude  greater  than  the  conventional  diffusivity.  The  diffusion  velocity  v  •  »'(r/t)  .  j.JxIO’cm/s  is  four 
times  larger  than  the  initial  Fermi-velo"ity.  It  shoulJ  Be  pointed  out  ‘rat  the  values  for  D  and  t  deduced 
from  the  time-resolved  luminescence  profile  Ise  with  the  theoretic, al  expression  Eq.  (lb)  are  consistent  with 
time-integrated  spatial  luminescence  profile.  Second,  the  carrier  lifetime  at  t<25  ps  deduced  from  the 
diffusion  Eq.  (13)  is  the  same  as  the  carrier  depletion  time  determined  from  the  tixe-rescl ved  luminescence 
data  shown  in  Fig.  4.  Thl3  demonstrates  that  our  analysis  is  sel  r-consistent .  Using  the  determined 
n(r,t)  we  can  3how  the  contribution  of  the  second  term  m  the  right-side  of  diffusion  Eq.  (13)  in  losing  the 
carrier  at  a  given  position  r  is  always  smaller  than  the  first  tern  and  thus  the  ultrashort  carrier  depletion 
time  cannot  be  attributed  to  the  carrier  diffusion  tut  to  the  nonequiliorium  phonon  enhanced  phonon  replica 
emission. 

Several  possible  origins  of  the  extraordinary  rapid  carrier  diffusion  in  the  x-y  pi  me  of  the  asym, metric 
quantum  well  can  be  given  at  this  stage.  (1)  The  internal  pressure10  of  high  electron-hole  pairs  may  be 
extremely  large  when  the  density  is  well  3bove  the  electron-hole  liquid  density.  (2)  There  are  a  large 
number  of  the  coherent  nonequilibrium  phonons  emitted  By  hot  electrons  in  the  system.  These  phonons  have 
momenta  parallel  to  the  x-y  plane  and  may  also  drive  the  carriers  outside  the  excited  region  by  adding  the 
phonon  momenta  to  the  carriers.  (3)  Due  to  the  asymmetric  potential  well  profile  there  are  no  transmission 
states  (virtual  3tates  or  unconfined  states)  above  the  energy  barriers  unlike  the  symmetric  well  case.  This 
13  the  intrinsic  property  of  asymmetric  quantum  wells.  The  photogenrrated  carriers  may  be  Immediately 
scattered  by  the  higher  energy  barrier  and  obtain  the  momentum  from  the  interface  which  assists  carriers  to 
expand  in  the  x-y  plane.  Due  to  a  large  kinetic  energy  of  -300  meV,  the  photogenerated  carriers  acquire  a 
la  ge  initial  diffusion  velocity.  The  origins  (1)  and  (2)  seem  to  be  inconsistent  with  the  symmetric  well 
case.  However,  whether  the  mechanism  described  in  (3)  is  responsible  Tor  the  observed  rapid  diffusion  is 
still  an  open  question  worthy  of  further  study. 
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Conclusion 

In  this  paper  extensive  data  have  been  presented  on  time-resolved  photoluminescence  studies  as  well  as 
other  related  investigations  of  a  symmetric  MQW  structure  and  an  asymmetric  quantum  well  with  similar  well 
thickness.  With  2  ps  time  resolution  we  were  able  to  study  the  initial  energy  relaxation  and  carrier 
diffusion  processes  simultaneously  under  the  same  conditions. 

The  existence  of  a  large  population  of  nonequilibrium  phonons  in  highly  excited  semiconductor  quantum  well 
structures  is  experimentally  verified.  Its  effect  on  the  energy  relaxation  is  to  slew  down  the  cooling  rate 
after  an  initial  rapid  cooling  (0-5  ps).  The  symmetric  and  asymmetric  potential  well  profiles  do  not  show 
any  difference  in  the  energy  relaxation  process  under  the  same  conditions. 

A  new  mechanism  to  explain  the  ultrashort  carrier  depletion  time  deduced  from  the  consistent  fitting  of 
time-resolved  photoluminescence  profiles  at  various  emitted  photon  energies,  is  proposed  to  be  nonequilibrium 
phonon  enhanced  phonon  replica  emission  rather  tnan  other  nonradiative  processes  suen  as  Auger  recombination, 
carrier  diffusion,  and  the  conventional  radiative  stimulated  emission. 

The  differences  in  carrier  diffusion  for  two  different  potential  well  shapes  are  remarkable.  For  the 
first  time  we  have  observed  the  ultrafast  carrier  diffusion  in  the  x-y  plane  of  the  asymmetric  quantum  well, 
while  there  is  little  or  no  carrier  diffusion  except  for  carriers  at  the  subband  edge  in  the  x-y  plane  of  the 
symmetric  quantum  well  under  similar  experimental  conditions. 
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Fig.  1  The  schematic  potential  profiles  of  a 
symmetric  and  an  asymmetric  quantum  wells.  The 
excitation  photons  are  absorbed  in  both  of  the 
barriers  and  the  wells  for  symmetric  wells, 
while  for  the  asymmetric  well  the  absorbtion 
occurs  in  one  of  the  barriers  and  the  well  . 
The  photoluminescence  originates  from  the 
recombination  of  n*l  electrons  and  n  =  l  heavy- 
hole  as  well  as  light-hole. 
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Fig.  3  The  triangles  are  the  ratios  of  the  peak 
intensity  of  B  and  A  indicated  in  Fig.  2.  The 
solid  curve  is  explained  in  the  text. 
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Fig.  2  The  time-integrated  luminescence  spectra 
at  various  lattice  temperatures.  The  peaks  A, 
8,  and  C  are  explained  in  the  text. 


f'K.  9.  Time-resolved  photclumineseence 

profiles  (solid  curves)  from  the  symmetric  MQW 
at  9K  at  various  wavelengths.  The  narrow  peaks 
within  0-10  ps  are  the  prepulses  used  for 
average.  The  next  peak  of  the  aottra  curve  is 
the  Rayleigh  scattering  light  from  the  sample 
surface 
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Fir.  5  Time-resolved  pnotolur.ir.escpr.ee  from  the 
asymmetric  quantum  well  at  «K  at  various 
wavelengths. 
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Fir,.  '!  Experimentally  determined  quasi 
energies  for  electrons  (up  solid  curve)  a 
holes  idown  suliu  curve). 


Fig.  6  Experimentally  determined  carrier 
temperature  a3  a  function  of  time.  The  shaded 
area  indicates  the  extent  of  uncertainty  in 
deducing  carrier  temperature  within  the  first  k 
ps. 


Fig.  o  Experimentally  determined  net  energy 
loss  rates  for  electrons  (e)  and  heavy-holes 
(hh)  obtained  from  Eq.  (7).  The  Cot-dashed  curve 
is  obtained  from  Eq.  (8).  The  broken  line  is  the 
theoretical  result  for  bulk  GaAs. 


'  SPIE  Vol  793  Ultralast  Lasers  Probe  Phenomena  in  Bulk  and  Microstructure  Semiconductors  1 1987 1 


.'•'Ok-o.- v.v.v.v. 


■450 


150 


4  50 


i? 


T  (om ) 


Fig.  13  Tine-intefrat.ee  spatial  profiles  cf  r:c.  '  ■  Tire-resoivea  snetei  unirescer.ee  cf  • 

laser  and  the  luminescence  fron  the  asymmetric  asymmetric  weli  at  -K  from  different  spat 

well  at  various  wavelengths  at  vk.  r-mcr.s  as  eistir.ruishe :  ay  ,  -r.d  lr.,. 

solid  „rves  are  tr.e  theoretic  a-  fit  ty  using 
ff)  me  i-d.  ilo)  for  :5,  and  1 ,  respectively. 
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Abstract 

The  non-commutation  relationship  between  the  Hamiltonian  of  the  system  and  the  electron  momentum 
operator  in  the  confinement  --direction  was  taken  into  account  using  random  variables  in  the  calculation  of 
electron  scattering  by  polar  optical  phonons  in  quantum  wells.  The  rate  at  the  onset  of  phonon  emission  is 
found  to  Pe  i  times  smaller  tnan  previously  reported.  Furthermore,  this  rate  becomes  on  increasing  function 
of  initial  electron  energy  instead  of  a  decreasing  function. 

It  is  important  to  take  into  account  tr.e  effect  of  electron  confinement  on  the  scattering  of  electrons  by 
polar  optical  phonons  (PC?)  in  quasi-two-dimensional  (Q2D)  quantum  wells.1-*  The  Hamiltonian  H  of  the 
electron  system  and  the  oper  itor  k.  for  the  electron  momentum  in  z-d;roction  do  not  commute  with  eacn  ''then 
because  of  tne  potential  profile  experienced  by  electrons  in  the  confinement  z-directior..  That  is, 


where  vTz;  is  tne  potential 


m  z-directior.. 


The  usual  conservation  law 


momentum  in 


(1) 

-direction: 


are 


must  be  modified  m  the  sense  that  I -. ,  k1-  and  k-.  are  r.o  longer  deterministic  otservab 
random  var.abies  ir.v.:.  Here,  qz  is  the  z-compcnent  of  tne  pnon-n  wave  vector,  k.  ir.d 
wave  vectors  before  and  after  scattering,  respectively.  The  statistics  of  each  r.v.  is  characterize; 
probability  density  function.  The  pioneering  won-  of  'Hess',  Price2  and  Leourton*  essentially  adopted 
momentum-selection  rule”  m  z-direction  for  calculation  of  phonon  emission  and  phonon  absorption  rate: 
treat  qz  as  a  regular  det  -  rm.in.c  1 1  c  variable*.  Althouvn  -  i-deiy'  realized  "fuzzy  “amentum  conservatic 
direction,  he  used  instead  a  "momentum  conservation  approximation".  Furthermore,  the  value  of  k, 
function  in  the  coordinate  representation  was  considered  is  tne  -lectror.  momentum  in  z-directior 
analysis.  This  approicn  is  not  valid  since  is  equal  to  u  rr-  use  value  of  n/u  in  an  "infmn 
approximation  and  for  a  "finite  well"  kz  is  just  a  parameter  that  related  to  the  eigen- value 
Hamiltonian  of  tne  electron  system. 


(2) 

but  are 
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The  thesis  of  this  paper  is  to  present  a  theoretical  analysis  of  electron  scattering  by  POP  in  Q2D 
structures  placing  emphasis  on  confinement  effects  by  taxing  into  account  the  non-commutation  relationship 
between  the  Hamiltonian  -per .at or  H  and  the  electron  momentum  opera*  or  k..  We  sh"w  new  "fuzzy  mementum 
conservation"  in  z-dirertion  is  exactly  treated  using  n-;:*ner  the  "non-momentum-seiection  rule"  n*“T  the 
"momentum  conservation  approximation"  but  by  ccm.side-.ng  ~t/  is  i  random,  variable.  A.  key  result  that 


distir.guisnes  this  work  from,  tne  previc 
reduced  by  a  factor  of  3. 


-ru,!-* 


rat  e 


"jno 


emission  is 


There  are  two  important  aspects  of  tne  non-c'mmut  a?  ion  relatinshsp  given  by  Eq.  (1).  First,  the 
corresponding  observables  cannot  be  measured  simultaneously  with  arbitrary  accuracy.  The  product  of  their 
quantum  root-nean-sqjare  deviation  must  rot  Oe  less  than  toe  half  of  tr.e  quantum  expectation  value  the 
d  V  ( z ) 

operator  i  in  a  given  state.  Second,  the  electron  energy  due  to  the  motion  in  z-direction  has  no 

correlation  with  the  corresponding  electron  moment  j-  in  this  direction.  The  eigen-value  of  momentum 
operator  kz  cannot  precisely  be  determined  due  to  the  electron  confinement  by  the  potential  well  and  must  be 

ted  m  the  f irsl  bound  state  the  quantum 
ii,kzj  ia  zero  resulting  in  j  precise  det^r r.inat . on  of  trie  lowest  eigenvalue  of  operator 

hose  two  points  play  the  key  role  in 


considered  as  a  random,  variable  (r.v.).  When  an  electron  is 
expectation  v ^ i ue  or  l H , k ^  j  1 3  zu  r  o  rosuitir. ^  in  a  pro c  *  3  c  C • 
H  and  root-mean-square  deviation  for  kz  car.  b 
calculating  the  electron-POP  scattering  rate. 


•ter-med. 


In  our  analysis,  the  unscreened  electron-.":?  interaction,  equilibrium  bulk  phonon  spectrum  and  an  infinite 
depth  square  well  with  the  lowest  subband  occupied  are  assumed. 


Within  the  framework  of  the  ef  f  ect  i  ve-mass-approacn , 
representation  is  then  given  by 


| k  >* 


ai  kr 
/  A 


0(3 


the  wave  function  for  electron  in  coordinate 


(3) 
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where 


iU)-\l^-  sin  ~ 


W 


is  the  wave  function  in  z-direction  describing  the  bound  state.  The  vectors  r=(x,y)  and  k  =  (kx,ky)  denote  the 
two  dimensional  electron  position  and  the  wavevector  parallel  to  the  interfaces,  respectively.  The 
corresponding  electron  energy  can  be  expressed  as 


where 


h2k2 

E=  +  F 

2m  tz’ 


F  M  " 

2  2mL  2' 


(5) 


(6) 


Here  m  is  the  isotropic  electron  effective  mass,  Lz  is  the  well  width  and  A  is  the  layer  area. 

The  probability  density  function  for  the  r.v.  kz  can  be  obtained  from  the  wave  function  given  by  Eq.  ( U )  by 
a  representation  change  to  the  wavevector  space  using 


♦(kz)  =  - 


/2i 


t>(z)e'il<zZ  dz. 


This  yields 


«LZ  eos2(kzLz  2) 

P(k2)  -  |<t><k,)|2  =  Va-k  2l,2  ^yr 


(7) 


(8) 


Since  there  exists  a  "fuzzy  momentum  conservation  law"'  in  z-direction  given  by  Eq.  (2),  the  probability 
density  function  for  the  r.v.  qz,  p(qz),  can  be  obtained  from  the  convolution  expression: 


p(qz)  1  p(  ikz)  *  p(  *k  ’z ). 


(9) 


The  electron  scattering  rate  (wk.k)  per  unit  area  in  k-space  in  Q2D  structures  can  be  expressed  in  Born 
approximation  by  Fermi's  golden  rule  by: 


2k'k 


U  -n  A 


h(2 


TP  I  ‘Vk  |Jt(Ek  -  -  h  u0). 


where 


with 


-c(n0*l  )1/2 


6k  k,qf79z> 0Z ^ 


c  -  i( 


e*hw0 


cV. 


(10) 


(ID 


(12) 


is  the  matrix  element  for  the  Frohlich  electron-POP  interaction.  For  simplicity  we  have  only  considered 
phonon  emission  process.  Q  and  Vs  denote  the  three  dimensional  phonon  wavevector  (q,qz)  and  the  sample 
volume,  respectively.  The  6  is  a  modified  Kronecker  function  defined  as 


4a-b  "  f  ir  a»b 
(0  if  a-b 

f. notion  F  in  Eq.  (11)  is  the  form  factor  in  z-direction  which  is  defined  as: 


(13) 


F(qz,Lz)  -  2. 

Lz 


Lz 


e"1(1zZ  sin2  dz. 


(1H) 


here  is  to  realize  that  the  matrix  element  Mkik  is  the  function  of  a  r.v.  qz. 
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Therefore,  the  Wk'k  also  becomes  a  r.v..  The  most  likely  value  of  the  wk.  k  that  can  be  experimentally 
measured  is  given  by  its  expected  value  using  the  statistics  associated  with  the  r.v.  q,.  This  is. 


r 


wk ' k  *  p(Qz)wk'k  (QzMqz. 


05) 


The  total  emission  rate  is  then  given  by 

wk'k  “  >  wk ' k  d  q  -/  /  I^iq.'-z'^Ek'Ek'-h  w0)V-k»q  Tddd  s 

where 


0  6) 


IQ2D(q.U)  -  lo 


g(x,y)dydx, 


07) 


and  0  denotes  the  angle  between  the  k  and  q.  The  factor  I„  is  defined  3S: 

'  •d)e2q(Jwo 


Io  * 


0  8) 


2ue 


The  function  g(x,y)  is  given  by: 


g(x,y  =  0.1 


d7** 


sin2(  —■) 


1  -(xz  /  2ir)2"j2(xz  /  2) 2 

r  ^ 

COS2  (x*y)z  /  2  cos2(yz  /2) 
-(x*y  )2Z2  /"^l-y^2  /  n2j2 


09) 


with  z*q„Lz,  which  includes  1  /  Q2  from  Frohlich  electron-POP  interaction  Hamiltonian,  the  form  factor  in  z- 
direction  and  the  probability  density  function  p(q2)  for  r.v.  qz.  All  the  wave  vector  are  normalized  by 


/  2m 


do* 


hz 


L„~‘  which  is  about  1  /  UoX  for  GaAs. 5 


The  integral  lQ2D  given  gq.  (17)  for  Q2D  structures  represents  the  electron-POP  interaction  strength  as 
function  of  q.  It  contains  all  the  information  about  confinement  effects  as  well  as  size  effects'  on  the 
scattering  rates.  In  order  to  further  understand  how  electron  confinement  makes  the  difference  in  the 
electron  scattering  rate  between  Q2D  structures  and  its  bulk  counterpart,  a  plot  of  the  electron-POP 
interaction  strength  l3D  versus  q  for  3D  electron  system  is  needed  in  such  a  form  that  can  be  compared  with 
I2D.  It  can  be  shown  that 


l3D(q) 


q 2  *  Ak  z  ( ~ ) 2  q 2  *  Ak  z  ( ' ) 2 


(20) 


where  Akz(')-(kz*k’z).  The  first  and  second  terms  on  the  right  side  of  Eq.  (20)  represent  the  scattering 
rocesses  without  and  with  change  of  initial  direction  of  electron  momentum,  respectively.  Since  Ak,(*)  in 
q.  (20)  depends  on  both  the  magnitudes  as  well  as  the  scattering  angle  6,  one  cannot  find  an  explicit 
xpression  for  for  Akz^  )  in  terms  of  q  and  8  unless  a  variational  parameter  for  the  electron  kinetic  energy 


I! 


ratio  (r-  K2  /  kz2)  between  the  total  kinetic  energy  and  the  energy  due  to  the  z-direction  motion  is 


introduced . 


Using  the  momentum  and  the  energy  conservation  laws  for  the  3D  scattering  process,  the  Ak,(+) 


expressed  as 


can  be 


Ak 


.O 


(q2*£.)»  /(q2*E)!-6‘)q2cos28 


2qcos  0/ r- 1 


(21) 


where  OSOSarccos 


.  /  hiuo 


£"r~r  with  1  <  r  <  ®.  The  value  of  r  =  3  implies  that  the  kinetic  electron 
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energy  is  equally  distributed  in  all  directions,  while  r<3  more  electron  kinetic  energy  is  in  the  z-direction. 
The  2D  electron  motion  in  x-y  plane  corresponds  to  the  case  when  r  approaches  infinite. 


We  now  present  the  numerical  results  in  the  following  sections. 

The  two  solid  curves  in  Fig.  1  show  the  probability  density  profiles  for  the  r.v.  qz  obtained  from  Eq.  (9) 
for  well  width  Lz  *  2L0  and  5L0.  The  larger  well  width  Lz,  the  smaller  is  tne  variance  of  the  density  curve 
resulting  in  the  less  fuzzy  momentum  conservation  in  z-direction.  The  magnitudes  of  F(qz,Lz)  versus  q,  for 
Lz  =  2L0  and  5L0  are  also  displayed  in  Fig.  1  as  two  dotted  curves.  The  two  features  from  the  plots  should 
be  noticed.  First,  in  the  limit  qz-*0,  |F(qz,Lz)|  approaches  unity  independent  of  Lz  reflecting  a  scattering 
event  with  an  exact  momentum  conservation  in  z-direction.  This  is  also  the  most  favorable  process  according 
to  the  density  curve  for  qz.  Second,  for  very  thick  well  the  |F(qz,Lz)|  can  be  written  as  6k  _k,^  .  which 
implies  the  momentum  conservation  law  for  any  scattering  event  and  p(qz)  approaches  a  i-f urTc'tiih  .or  the 
bulk  situation.  The  fact  that  not  every  scattering  event  follows  the  exact  momentum  conservation  law  in  z- 
directron  is  due  to  the  electron  confinement.  This  procuces  the  weakening  of  electron-PCP  interaction  in  Q2D 
structures. 


The  solid  curves  in  Fig.  2  show  the  integral  lQ2D  as  function  of  the  parallel  component  of  phonon  wave 
vector  q  with  Lz  as  a  parameter.  Curves  1,  2,  3,  and  k  are  plotted  in  Fig.  2  for  the  value  of  L,  of  C.5L0, 
L0,  3L„  and  5L0,  respectively.  Our  calculation  procedure  using  the  probability  density  function  p(qO  for  r.v. 
qz  yields  a  significant  reduction  by  a  factor  2  on  the  value  of  the  integral  I'3-D  obtained  by  Leburton*  in  the 
range  from  0.1  to  k  for  q.  In  spite  of  this  key  difference  other  features  of  our  integral  r=2D  are  very 
similar  to  what  Leburton  has  shown."  Namely,  the  integral  lG2D  is  a  monotonic  decreasing  function  of  q  for 
Q2D  electron  system  and  weakly  dependent  on  Lz.  This  similarity  occurs  just  because  the  profile  of 
probability  density  function  for  qz  happens  to  have  similar  shape  as  |F(qz,Lz)|. 

The  interaction  strengths  1 3 D  for  30  system  with  various  values  of  9  and  r  are  also  displayed  in  Fig.  2  by 
the  dotted  curves  5  to  8.  The  values  of  r  for  curves  5,  6,  and  8  are  10,  3,  and  1.2,  respectively.  The 

value  of  0  for  those  three  curves  is  set  to  be  zero.  The  curve  7  with  r  =  l.l  and  6  =  71/6  is  used  to  see  the 
variation  of  l3 D(q )  by  the  change  of  the  scattering  angle  6.  A  comparison  of  the  dotted  curves  5  to  3  for  3D 
with  the  solid  curves  for  Q2D  provides  the  physical  insight  into  electron  confinement  in  quantum  wells.  The 
following  three  points  are  discussed:  first,  the  value  of  l3D  increases  with  increase  of  the  value  of  r.  It 
indicates  that  the  electron-POP  interaction  strength  for  given  value  of  q  in  the  system  with  r>3  is  stronger 
than  the  system  with  r<3.  The  reason  for  this  is  the-weaxness  of  electron  motion  m  z-direction  resulting  in 
small  amount  of  change  of  the  electron  momentum  &k,(  K  It  should  be  emphasized  that  the  value  of  I'S2-1  at  q 
as  Lz->0  do  not  approach  the  value  of  with  r-*«  which  corresponds  to  tne  case.  In  other  words,  as  far 
as  the  electron  scattering  process  is  concerned  the  C2D  quantum  well  structure  with  Lz  -*0  is  not  equivalent 
to  a  2D  system  due  to  the  potential  barriers  which  confine  electrons  inside  the  well.  The  values  of  l3°  and 
lQ2D  are  close  to  each  other  for  the  value  of  the  ratio  r  between  1.2  to  1,5.  This  implies  the  origin  of  the 
difference  in  electron  scattering  between  Q2D  and  3D  systems  comes  from  tne  fact  that  the  electron  kinetic 
energy  is  not  equally  distributed  in  all  directions  arising  from  electrons  in  the  well  obtain  extra  kinetic 

energy  in  z-direction  from  the  potential  barriers.  These  are  important  findings  which  were  overlooked  in 

previous  works. 1_* 

Second,  the  larger  scattering  angle  0  yields  a  larger  value  of  I^.  For  q >  1  ,  i.e.,  for  the  scattering 

processes  in  which  electrons  change  their  initial  momentum  direction,  the  value  of  l33  iS  almost  independent 

of  0  and  can  be  distinguished  from  various  values  of  r  by  the  second  term  of  the  right  hand  side  of  Eq.  (20). 
The  value  of  lQ2D  is  independent  of  0  and  insensitive  to  the  Lz  for  q<i  giving  rise  to  a  weak  Lz  dependence 
of  total  electron  scattering  rates.  It  also  appears  that  the  value  of  I'*2D  for  q>l  is  more  sensitive  to  L 
than  for  q<1.  z 


Third,  a  dent  occurs  in  each  dotted  curve  which  becomes  smaller  when  r  increases.  The  position  of  the 
dent  moves  towards  the  smaller  q  when  r  decreases.  It  should  be  noticed  that  the  dent  does  not  appear  in 
solid  curves  for  Q2D  structures.  This  is  because  the  non-commutation  relationship  given  by  Eq.  (1)  employed 
for  Q2D  systems  which  gives  rise  to  the  momentum  change  in  the  z-direction  does  not  correlate  to  the 
momentum  change  in  the  direction  parallel  to  the  quantum  well  interfaces. 

The  total  normalized  phonon  emission  rate  (Wk.  «/•»„)  m  Q2D  quantum  wells  versus  normalized  energy  Ek/huj0 
Is  displayed  In  Fig.  3  by  two  solid  curves.  The  normalization  factor  W„  is  given  by  w,=2r I0/hJq0!.  For  GaAs 
W»-1.2k  ps".  For  comparison,  the  calculated  results  for  Q2D  by  Leburton"  (dash  curves)  and  for  3D  (dotted 
curves)  are  also  plotted  in  Fig.  3.  The  results  of  our  approach  clearly  differ  from  the  others1**  in  two 
aspects.  First,  the  Q2D  phonon  emission  rates  for  various  Lz  are  3  times  smaller  at  the  onset  of  phonon 
emission  (Ek/h<j0-l)  than  previously  obtained.1  "  Second,  the  phonon  emission  rate  is  a  slowly  increasing 
function  of  Initial  electron  energy  instead  of  being  a  decreasing  function.  The  non-zero  value  of  the  phonon 
emission  rate  for  Ek=hiu()  stems  Trom  the  nature  of  Q2D  density  of  states.  For  thick  quantum  wells  (L  >5L0) 
the  phonon  emission  rate  Is  about  a  factor  of  2  smaller  than  the  3D  rate  for  Ek/hw„>2  resulting  rrom  the 
weakening  of  the  electron-POP  Interaction  in  the  z-direction  for  Q2D  case  and  neglecting  higher  intra-subband 
and  Inter-subband  contributions.  The  Q2D  emission  rate  will  approach  the  3D  result  when  Lz->-  if  all  the 
contributions  mentioned  above  are  included.  The  phonon  emission  rate  is  comparable  to  the  3D  rate  for  very 
thin  quantum  wells  arising  from  the  competition  between  Q2D  density  of  states  and  the  weakening  of  electron- 
POP  Interaction  strength  In  z-dlrectlon. 
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For  completeness,  it  should  be  pointed  out  that  the  phonon  absorption  scattering  rate  can  be  treated  in 
the  same  way  as  outlined  in  this  paper.  One  should  expect  that  the  emitted  phonon  by  Q2D  electrons 
traveling  in  the  direction  perpendicular  to  the  quantum  well  interface  may  be  easily  reabsorbed.  This  is 
because  those  phonons  are  not  plane-wave-iike  phor.on  mode  that  is  extended  to  the  whole  sample.  The 
phonons  are  localized  inside  the  well  the  same  as  the  localized  electrons. 


In  conclusion,  the  confinement  effects  on  the  phonon  emission  rate  of  the  electron  in  Q2D  structures 
through  tne  unscreened  electron-POP  interaction  are  shown  by  taxing  into  account  the  non-commutation 
relationship  between  the  Hamiltonian  of  the  system  and  the  electron  momentum  operator  in  the  z-direction. 
It  is  explicitly  remonstrated  how  one  can  exactly  treat  the  fuzzy  momentum  conservation  in  z-direction  using 
neither  "non-momentum-selection  rule"  nor  "momentum  conservation  approximation"  by  considering  the  phonon 
wave  vector  qz  as  a  r.v..  The  phonon  emission  rate  at  the  onset  is  found  3  times  smaller  than  previously 
reported  and  increases  with  initial  electron  energy. 


This  work  was  supported  by  the  Air  Force  Office  of  Scientific  Research  grant  #86-003'.  The  authors  would 
like  to  thank  Jianping  uu  and  Jimmy  Zheng  for  their  help  in  computer  analysis. 
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Fig.  2  The  electron-PCP  interaction  strength  for  Q2D  and  3D  systems  versus  the  parallel  component  of  phonon 
wave  vector  q.  The  solid  curves  from  1  to  »  are  for  the  32D  structures  with  well  width  0.5  L0,  L0,  3L0l  and 
5L0  accordingly.  The  main  contribution  of  the  integral  IQ2D  to  the  total  phonon  emission  rate  occurs  at  q-1, 
where  lQ2D  presents  a  weak  Lz  dependence.  The  dotted  curves  from  5  to  7  are  for  the  3D  system  with  ratio 
r-10,  3,  and  1.2,  respectively.  The  scattering  angle  9  for  those  three  curves  are  set  to  be  all  zero.  The 
dotted  curve  8  is  also  for  3D  system  with  r-1.2  and  9-ir/6.  The  larger  scattering  angle  6,  the  larger  the 
value  of  l3D  in  the  range  from  0.1  to  1  for  q. 
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Fig.  3  The  normalized  phonon  emission  rates  as  function  of  Ek/hu0  in  semiconductor  structures  at 
temperature.  The  solid  curves  are  from  this  work  and  the  dash  curves  were  generated  from 
result."  The  well  thicknesses  are  indicated  on  the  curves.  The  dotted  curve  is  the  3D  result. 
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Abstract 

The  recombination  lifetimes  for  the  radial  and  angular  quantum  number  conserving  1S-1S  and  1P-»1P  transi¬ 
tions  from  three-dimensionally  confined  electrons  in  CdSxSe,_x  were  measured  by  time-resolved  photolumines¬ 
cence.  The  assignment  of  the  observed  transitions  is  supported  by  calculations  of  eigen  energy  levels, 
squared  matrix  element  ratio  for  those  transitions,  steady-state  and  picosecond  spectroscopic  studies. 

Intr o  d  uc tion 

Carriers  localized  in  semiconductor  microstructures  have  attracted  much  attention  both  for  fundamental 
reasons  and  for  their  potential  applications.  Semiconductor  crystallites  in  a  transparent  insulating  matrix 
can  confine  electrons  in  three  dimensions  because  their  diameters  range  from  30  X  to  800  H. 1-1  The  envelope 
wave  functions  in  spherical  coordinates  [r,  6,  0]  and  eigen  energies  of  conduction  electrons  and  valence 
holes  localized  in  an  infinite-spherical-well  within  the  effective  mass  approximation  are  given  by: 

|nLm  >  -  CnL  jL(XnLr>  *Lm  (6.*).  (U 


and 


EnL 


2me,ha  1  ’ 


(2) 


respectively.  The  subscripts  n,L,  and  m  are  effective  radial,  angular  and  magnetic  quantum  numbers,  respec¬ 
tively.  j  and  Y  are  the  spherical  Bessel  and  spherical  Harmonic  functions,  respectively,  a  is  the  radius  of 
crystallite  and  me  p  is  the  effective  mass  of  electron  or  isotropic  hole  mass.5  The  value  of  xnL  *’or  the 
lowest  two  states  'of  either  conduction  electron  or  valence  hole  IS  and  IP  are  n  and  b . h 9 ,  respectively.  S 
stands  for  L*0  and  P  for  L*1.  The  allowed  transitions  which  conserve  angular  and  radial  quantum  numbers  are 
1S-»1S,  1P-»1P,  and  higher  transitions.  The  physical  picture  of  the  quasi-OD  electron  system  described  above 
has  been  experimentally  verified  by  several  groups. 1-5  Large  and  fast  optical  nonlinearities  arising  from 
photogenerated  electron  hole  pairs  have  been  also  reported  in  quasi-OD  electron  systems. 6-7  Pioneering  sub¬ 
nanosecond  time-resolved  photoluminescence"  and  picosecond  pump-probe  experiments’  have  been  performed.  It 
is  important  to  study  the  dynamics  of  various  excitations  in  quasi-OD  electron  systems  in  detail  in  order  to 
understand  the  physical  origins  underlying  these  various  novel  optical  properties. 


In  this  paper,  we  report  on  measurements  of  ultrafast  recombination  lifetimes  for  the  radial  and  angular 
quantum  number  conserving  1  S-»  1 S  and  1 P -» 1 P  transitions  in  quasi-OD  electron  systems  in  CdSxSe,_x  at  *).3K 
using  a  streak  camera  detection  system,  and  on  the  observation  of  optical  transitions  using  photoluminescence 
spectroscopy  between  quantized  levels  (IS,  IP)  in  the  conduction  and  valence  bands  in  quasi-zero  dimensional 
electron  system  in  CdSxSe,_x.  The  time-resolved  photoluminescence  detected  at  various  emission  ener¬ 
gies  allows  us  to  unambiguously  identify  the  1S-1S  and  1 P -» 1 P  transitions.  The  recombination  lifetime 
of  1P-*1P  transition  was  measured  to  be  3-5  times  shorter  than  1S-1S  transition.  The  ultrafast  decay  of  IP 
excitation  may  have  practical  importance  for  the  construction  of  ultrafast  reversible  optical  switches.  The 
assignment  of  the  observed  S  and  P  transitions  were  supported  by  calculations  of  eigen  energy  levels 
and  squared  matrix  element  ratios  for  these  transitions,  as  well  as  from  3teady-state  photolumines¬ 
cence  spectra  of  the  same  samples  used  for  the  picosecond  kinetic  3tudie3. 
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Ex£e  r  i_m  ent a_l  _Method 

The  samples  Investigated  were  four  optical  glass  filters  2-6! ,  2-59,  2-58  and  2 -6 h  manufactured  by  Corn¬ 
ing  years  ago.  The  samples  were  mounted  on  a  cold  finger  of  a  variable  temperature  helium  optical  cryostat 
(2K  to  300K).  The  x  value  for  each  sample  was  accurately  obtained  from  chemical  analysis.10  The  picosecond 
experimental  setup  used  in  this  research  has  been  described  in  detail  elsewhere.11  A  second  harmonic  (530 
nm)  of  a  Nd-gla3s  laser  pulse  of  8  ps  duration  was  used  to  excite  the  samples  on  the  front  surface.  The 
maximum  optical  energy  Incident  onto  the  front  surface  was  about  90  pj.  This  corresponds  to  101*  photons 
per  single  shot.  The  spot  size  was  8x10~5  cm2.  The  exact  photoexclted  carrier  density  in  quasi-OD  electron 
system  was  difficult  to  estimate  due  to  the  complexity  of  the  surface  of  the  crystallites.  Assuming  all  the 
photons  were  absorbed  with  10pm,  a  value  on  the  order  of  10*"  cm-5  was  estimated.  The  photoluminescence 
emitte  :  from  the  samples  wa3  collected  by  a  combination  of  lenses  and  imaged  onto  a  30-pm  slit  of  a  Hama¬ 
matsu  streak  camera  system.  Various  narrow-band  filters  centered  at  different  wavelengths  were  placed  In 
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front  of  the  streak  camera  to  select  different  transition  energies.  Two  color  filters  3*67  were  placed  in 
front  of  the  streak  camera  to  block  the  light  with  wavelengths  shorter  than  5*iOoJL  The  photoluminescence 
profiles  were  corrected  for  the  nonlinearity  of  streak  rate  and  spectral  response  of  the  detector  (S-20). 
For  steady-state  photoluminescence  experiments,  an  Argon-ion  laser,  a  double  Spex  spectrometer, 

S-20  photomultiplier  and  a  lock-in  amplifier  were  used. 


Picosecond  Luminescence  Studies 


The  photoluminescence  profiles  in  time  domain  were  measured  at  different  energies  covering  the  corres¬ 
ponding  spectral  range  of  each  sample.  The  expression  for  number  of  photons  detected  at  energy  E  at  time  t 
is  given  by: 

wp  ti  ,  u.r^-t/Thfb>  .  -t/tr-tE),  fn 


where  E  is  the  central  energy  of  the  i. arrow-band  filter  used;  a d ( E)  and  tr(E)  are  the  decay  and  the  rise 
times  at  energy  E;  and  N0  is  the  proportionality  constant.  The  value  of  tr  gives  information  about  how  fast 
electrons  and  holes  scatter  out  from  initially  excited  states  to  iess  energetic  states  where  they  recombine 
while,  gives  a  measure  of  the  recombination  lifetime  of  various  excited  states  when  recombination  domi¬ 
nates  the  depopulation  of  photo-excited  carriers. 


The  time-resolved  photoluminescence  profiles  obtained  at  various  emission  energies  at  H.2K  for  sample  2-58 
are  shown  in  Fig.  1(a-e).  The  left  curve  in  Fig.  la  shows  the  temporal  profile  of  exciting  laser  pulse  which 
reflects  the  time  resolution  of  detection  system  (10  ps).  The  dotted  curves  in  Fig.  l  are  the  fit  to  the 
data  using  Eq.  (3)  with  a  value  of  12  ps  for  the  risetime  tr.  The  decay  times  of  these  time-resolved  lumi¬ 
nescence  profiles  show  only  two  distinct  values.  For  emission  energies  ranging  from  2.213  to  2.175  eV  and 
from  2.1  to  1.967  eV  the  decay  times  are  29  and  100  ps,  respectively. 


The  decay  time  as  a  function  of  emitted  photon  energy  E  for  four  quasi-OD  electron  systems  are  plotted  in 
Fig.  2.  For  a  comparison  (Fig.  2e),  the  exciton  lifetime  vs  exciton  energy  in  quasi-JD  system  in  bulk 
CdS0, 5iSe0. . 7  while  the  other  two  samples  do  not  show  tr.is  feat  ure  (the  reason,  the  P  levels  have  not  been 
reached  in  2-61  and  2-59  samples).  It  is  expected  that  the  decay  time  should  be  nearly  invariable  over 
the  energy  range  of  emission  if  only  the  lowest  confined  state  13  for  electrons  and  holes  is  occupied.  The 
wide  spectral  range  most  likely  reflects  the  crystallite  size  distribution  and  the  fluctuations  in  the  value 
of  x  from  crystallite  to  crystallite.  When  the  two  lowest  states  13  and  IP  are  substantially  occupied  and 
the  recombination  lifetimes  associated  with  these  two  states  are  considerably  different,  a  steplike  change  of 
decay  time  will  be  observed  in  two  distinct  energy  regions.  The  exact  lifetime  ratio  in  the  two  energy 
regions  depends  on  the  transition  matrix  elements.  The  above  argument  should  in  principle  explain  what  we 
observed  in  Fig.  2.  However,  the  energy  dependence  of  idealized  exciton  lifetime  in  the  bulk  compound 
CdS„  sjSe0  ,i  as  displayed  in  Fig.  2e  also  exhibits  a  steplike  feature.  The  explanation  given  by  Kash  et  al.’ 
was  "based"  on  the  model  suggested  by  Cohen  and  Sturge"  where  the  exciton  migration  from  a  site  with  higher 
energy  into  another  site  with  lower  energy.  This  raises  a  question  whether  the  exciton  nigration  mechanism 
can  apply  to  the  present  case.  Since  exciton  migration  must  involve  a  transport  over  a  relatively  large  dis¬ 
tance  due  to  the  nature  of  the  exciton-phonon  interaction,"  such  migration  may  occur  in  relatively  large 
crystallites.  This  apparently  may  explain  the  results  observed  in  the  sample  2-58  and  2-6H.  In  smaller 
crystallites  the  migration  does  not  occur  in  correspondence  with  the  samples  2—6 1 .  However,  the  exciton 
migration  picture  cannot  explain  the  steplike  energy  dependence  of  lifetime  as  well  as  the  ratio  with  the 
value  of  about  3.5  because  the  depth  of  a  potential  well12  caused  by  compositional  fluctuations  are  entirely 
random  in  nature.  Furthermore,  the  exciton  lifetime  would  be  increased  monotonically  with  a  decrease  of 
exciton  energy  as  recently  observed  for  CdTe-ZnTe  interfaces. ' ’ 


A  two-state  (IS,  IP)  model  is  introduced  here  which  consistently  interprets  all  the  results  shown  in  Fig. 
2  as  well  as  other  results  observed  by  Cohen  and  Sturge."  It  is  possible  that  the  excitons  are  confined  by 
the  potential  wells  in  the  bulk  alloy  compound  which  have  two  lowest  states  designated  by  IS  and  IP.  The 
exciton  in  higher-lying  state  IP  can  be  scattered  to  the  lower-lying  state  IS  by  emission  of  phonons  instead 
of  migrating  from  site  to  site.  The  two-state  model  can  explain  not  only  the  results  of  Cohen  and  Sturge" 
but  also  our  time-resolved  dynamic  results  presented  here  as  well  as  the  data  reported  by  Kash  et  al..* 


In  order  to  further  verify  the  two-states  (IS,  IP)  model,  we  have  calculated  the  eigen  energy  levels  of 
the  qua3l-0D  electron  systems  and  the  recombination  lifetime  ratio  for  1S-1S  and  1 P -*  1  p  transitions.  The  con¬ 
finement  energy  for  the  ground  states  is  given  by 


E.S  -  E?S  *  eI?S 


for  the  qua3l-0D  electron  systems,  where  Efg  and  E^g  are  the  lowest  confinement  energies  for  electrons  and 
holes,  respectively.  These  energies  were  calculated  from  the  bulk  bandgaps'"  at  300K  and  the  measured  peak 
energies  of  the  fir3t  derivative  of  room  temperature  reflectance  (dR/dA).1!  Using  the  electron  effective 
mass  and  the  isotropic  hole  effective  mass  given  in  Ref.  R,  the  effective  diameter1*  of  crystallite  for  each 
sample  and  eigen  energies  were  calculated.  These  results  are  displayed  in  Table  1.  The  spin-orbit  split 
states'"  were  neglected  since  they  were  not  excited  in  our  experiments.  The  emitted  photon  energies  at  9.2K 
for  1  S*1  S  and  1P-HP  transitions  were  obtained  by  adding  corresponding  confinement  energies  to  the  values  of 
bulk  gap17  at  H.2K.  The  binding  energies  of  excitons  were  not  included.  The  calculated  energy  positions  are 
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in  good  agreement  with  what  we  measured.  These  positions  are  located  in  Fig.  2  as  arrows  labeled  by  S  and  P 
for  clarity. 


In  Table  1,  the  values  of  lifetime  for  1S-1S  transitions  denoted  by  1,3  were  measured  at  the  luminescence 
intensity  peaks.  A  value  of  85  ps  for  the  1S-»1S  ground  state  recombination  time  for  sample  2-59  is  in  exact 
agreement  with  the  value  obtained  by  Warnock  and  Awschalom’  using  pump-probe  experiment.  Since  the  energy 
per  pulse  in  our  experiment  is  about  four  orders  of  magnitude  larger  than  they  used,  the  Auger  mechanism  of 
carrier  depopulation  is  then  ruled  out  in  these  quasi-OD  electron  systems.  Furthermore,  we  did  not  observe 
the  change  in  lifetimes  By  varying  the  excitation  intensity  over  a  factor  of  33.  One  concludes  that  the  dep¬ 
opulation  of  photoexcited  carriers  must  be  dominated  by  radiative  recombination  path  in  the  present  study. 
The  recombination  lifetimes  for  1 P -*  1 P  transitions  denoted  by  T,p  are  also  entered  in  Table  1.  Since  the  pre¬ 
sent  experimental  condition  cannot  detect  the  1P-*1P  transitions  in  samples  2-61  and  2-59  the  numbers  inside 
the  circles  are  the  anticipated  lifetimes  of  1  P-»  1 P  transitions  assuming  the  same  ratio  of  T,s/r,p  for  the 
samples  2-58  and  2-69.  Note  that  a  value  of  18  ps  has  been  also  reported  by  Warnock  and  Qwschalom’  for 
2-61  although  they  did  not  identify  this  to  be  the  lifetime  of  1 P  —  1 P  transition.  It  should  be  emphasised 
that  strongly  enhanced  transition  probability  which  leads  to  shortened  lifetimes  of  various  excitations  in 
the  quasi-zero  dimensional  electron  system  is  not  a  surprising  result.  The  recombination  enhancement  due  to 
carrier  confinement  in  quantum  well  structure  (quasi-two  dimensional  electron  system)  was  clearly  demon¬ 
strated  by  Gobel,  Jung,  Kuhl,  and  Ploog.'*  These  authors'*  attributed  their  measured  recombination  lifetime 
to  be  the  spontaneous  carrier  lifetime  in  quasi-two  dimensional  the  electron  system.  The  measured  diameter 
dependence  of  recombination  lifetime  supports  their  arguments.'* 


A  further  test  to  support  the  two-states  picture  is  whether  a  reasonable  prediction  can  be  made  to  locate 
the  energy  positions  of  exciton  localized  in  IS  and  IP  states  for  the  quasi-OD  electron  system  in  the  bulk 
compound  CdS0t  s  ,Se0.  ..  7  studied  by  Kash  et  al..*  When  electrons  are  confined  in  a  one-dimensional  semiconduc¬ 
tor  quantum  well  with  thickness  Lz,  the  electron-hole  recombination  time  for  ground  states  may  be  estimated 
to  be  proportional  to  the  Lz  since  the  wave  function  is  "squeezed"  by  the  boundaries.’-'"  The  evidence  for 
this  has  already  been  provided  by  the  measurements  of  Cooel  et  al..‘*  Extending  the  above  concept,  the 
effective  diameter  (d)  for  the  spherical  potential  well  may  De  estimated  using  a  relation  1,5  »  c  d1.  A 
value  of  about  250  A  for  the  diameter  of  the  potential  well  in  the  bulk  was  estimated  using  the  lifetime  for 
energies  less  than  2.09  eV  assigned  by  us  to  be  the  recombination  lifetime  for  1S-1S  transition.  Using  this 
value  the  confinement  eigen  energies  for  exciton  (electron  and  hole  Dound  By  their  mutual  Coulomb  potential) 
were  calculated.  These  results  are  also  entered  in  TaDle  1  under  Duik  a.nd  displayed  in  Fig.  2e.  The  experi¬ 
mental  results  of  Kash  et  al."  are  in  excellent  agreement  with  our  two-states  model.  Furthermore,  when  the 
temperature  of  their*  sample  was  raised  above  3oK  (3  new)  the  steplixe  feature  of  lifetime  in  two  distinct 
energy  regions  disappeared.  Since  exciton  binding  energy  in  this  material  is  about  20  meV,  the  exciton  disso¬ 
ciation  is  not  responsible  for  this  disappea-ance.  However,  the  energy  difference  between  ground  state  IS 
and  first  excited  state  IP  for  hole  is  just  3  meV  (Table  1);  therefore,  the  strong  hole  scattering  between  IS 
and  IP  3tates  assisted  by  thermal  energy  is  responsible  for  the  absence  of  the  two  distinct  regions  above 
36k.  This  temperature  dependence  of  exciton  lifetime  further  supports  our  model. 


The  lifetime  ratio  for  1 S -» 1 S  and  IP-lP  can  be  predicted  following  tne  approach  given  by  Casey  and  Panish” 
and  Efros  and  Efros.20  The  matrix  element  for  3  transition  from  localized  state  in  conduction  band  to  the 
localized  state  in  valence  can  be  expressed  as: 


where  is  the  average  matrix  element  for  the  31och  states  for  bands  in  absence  of  new  eigen  states  due  to 
the  confinement  and  Menv  is  the  envelope  part  of  the  matrix  element  which  is  given  by: 


<  1  L  |  1L>. 


Since  IP  state  is  a  threefold  degeneracy  state  with  magnetic  quantum  number  -1,  0,  and  1,  a  value  of  5  for 
the  ratio  of  |Menv  and  |Menv  |,2s+1s  can  be  readily  calculated.  This  value  should  be  compared  with 
experimentally  measured  value  of  3.5  for  1,3/1  ,p.  The  agreement  is  reasonable  in  view  of  assumptions  made 
in  the  theoretical  approach  and  the  uncertainty  in  data  fitting. 


St  e  a  dy_j  s  1 3  te_  _Photol  u  m  i_nes_cence_  _St  u  die  s 


Following  work  of  Warnock  and  Awschalom’  steady-state  photoluminescence  studies  was  used  to  confirm  our 
findings.  In  Fig.  3a,  photoluminescence  spectra  of  sample  2-6]  at  various  excitation  levels  taken  at  room 
temperature  is  plotted.  The  excitation  wavelength  was  the  988  nm  line  of  the  Argon-ion  laser.  The  excita¬ 
tion  levels  were  2.5  and  10  times  higher  for  the  broken  curve  and  dot-dashed  curve  than  the  solid  curve 
respectively.  The  scale  of  luminescence  Intensity  was  2  times  and  10  times  larger  for  the  broken  curve  and 
dot-dashed  curve  than  the  solid  curve,  respectively.  The  broken  vertical  line  indicates  the  peak  position  of 
dR/dA.  The  most  salient  feature  we  would  like  to  stress  is  the  appearance  of  the  two-peak  structure  at  the 
high  energy  side  of  the  peak  position  of  the  first  derivative  of  room  temperature  reflectance  spectrum. 
Based  on  results  of  picosecond  studies  listed  in  Table  1  we  attribute  these  two  peaks  to  the  1S-1S  and  1P-.1P 
optical  transitions’  as  Indicated  by  S  and  P  in  the  figure.  The  energy  separation  between  P  and  S  is  about 
220  meV.  This  rules  out  the  possibility  of  transition  from  IS  electron  state  to  IS  spin-split-hole  state, 
since  the  energy  separation  of  spin-orbit  splitting  is  about  350  meV  for  Se  mole  composition  of  0.73.  Two 
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other  features  should  also  be  mentioned.  The  excitation  dependence  of  luminescence  intensity  did  not  show 
the  saturation  of  available  electronic  states  up  to  the  maximum  excitation  level  ICGW/cm2;  md  the  lumines¬ 
cence  peak  shifts  to  the  lower  energy  when  the  excitation  level  increases  reflecting  the  renormalization 
process  in  quasi-zero  dimensional  electron  system. 


The  photoluminescence  spectrum  of  2-61  measured  at  U.3K  using  the  H57.9  r.m  line  of  the  Argon-ion  laser  as 
excitation  wavelength  is  shown  in  Fig.  lb.  As  can  be  seen  from  the  spectrum  the  position  of  1S-1S  peak 
snifts  to  the  higher  energy  and  its  FWHM  (58  meV)  decreases  by  a  factor  of  2  in  comparison  with  room  temper¬ 
ature  data.  The  latter  indicates  the  broadening  of  IS-is  transition  at  room  temperature  is  not  entirely 
dominated  by  size  fluctuation  of  the  microcrystallites  under  the  photoexcitation.  The  line  broadening 
mechanisms  may  apply  for  the  quasi-zero  electron  system.  We  assume  the  broadening  of  the  IS  peak  is  domi¬ 
nated  by  the  fluctuation  of  the  diameter  of  the  microcrystallites  at  A.  3  K.  The  inset  of  Fig.  lb  shows  a 
spectrum  fit  by  assuming  the  Gaussian  shape21  probability  density  function  for  the  random  variable  d.21  The 
fit  gives  a  mean  value  of  7k  \  for  the  effective  diameter  and  a  variance  a2  =  25  w  which  characterizes  the 
size  fluctuation  of  microcrystallites.  The  structure  beside  the  main  IS  peak  is  from  the  1P-1P  transition 
since  its  energy  separation  with  the  IP  peak  is  about  210  meV.  The  broad  peak  centered  at  507  nm  may  arise 
from  the  transition  from  the  split  orbit  band. 


In  Fig.  9,  the  photoluminescence  soectra  of  sample  2-59  measured  both  at  room  temperature  (a)  and  k.JK  (b) 
are  displayed.  The  most  important  feature  is  the  appearance  of  both  1 S -*  1 S  and  1 P -» 1 P  transitions.  The 
energy  separation  between  IS  and  IP  peaks  is  21 k  meV.  The  calibrated  value  for  the  confinement  energy  is 
213  meV.  The  corresponding  peak  positions  are  shifted  to  lower  energy  in  comparison  with  2-61  due  to  the 
larger  diameter  of  crystallites  and  smaller  mole  percent  of  the  composition  of  sulfur. 

The  broad  peak  centered  at  500  nm  in  Fig.  kb  at  the  k.3K  spectrum  appears  in  the  k.3K  spectrum  of  the 
2-6l  sample.  The  origin  of  this  emission  peak  at  500  nm  is  unlikely  from  the  ID-ID  transition.  This  emis¬ 
sion  may  arise  from  the  optical  transition  from  the  IS  electron  state  to  the  IS  spin-orbit  split  hole  state. 
The  energy  separations  between  the  weak  peak  and  the  main  IS  peak  are  about  370  meV  and  -50  meV  for  2-ol 
and  2-59,  respectively.  These  energies  are  consistent  with  the  measured  selenium-composition-percentage  dep¬ 
endence  of  spin-orbit  splitting.1" 

We  have  also  measured  the  room  temperature  spectra  of  2-58  and  2-6k  samples.  The  peaws  for  1-5S  and 
2-6k  are  at  1.9537  eV  and  1.911k  eV,  respectively.  These  peaks  were  broad  so  that  the  IP  structures  were 
not  clearly  resolved. 

Sumjnar  y 

In  summary,  we  have  reported  on  the  measurements  of  recombination  lifetimes  of  15-lS  and  )F*ip  transi¬ 
tions  in  quasi-OD  electron  systems  by  time  resolved  photoluminescence  measur ements.  The  assignment  of  these 
observed  transitions  are  supported  by  calculations  of  eigen  energy  levels  and  the  corresponding  matrix  ele¬ 
ment  ratio  as  well  as  the  direct  observation  of  optical  transitions  between  quantized  energy  levels  (IS,  IP) 
In  the  conduction  and  valence  bands  by  steady-state  photoluminescence  measurements. 
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Note:  After  the  completion  of  our  work,  we  recieved  a  preprint  from  Dr.  Hall,  Corning  glass  works  (to  be 
published  in  J.  Appl.  Phys.).  According  to  them,  the  wet  chemical  method  for  determining  x  value  from  0.3  to 
0.8  was  uncertain  by  15  J  from  the  x-ray  method  measured  corresponding  values.  We  should  state  that  the 
energy  peaks  we  measured  can  not  be  attributed  to  the  changes  of  x  determined  by  them  for  "a  given  filter 
number.  The  IBM  group  does  not  agree  with  Hall  et  al  x  values.  Therefore,  we  conclude  that  the  each  filter 
does  not  have  a  unique  x  even  though  the  designated  filter  number  may  be  same.  Our  response  to  Hall's 
statement,  "that  the  crystallite  size  in  filter  glasses  are  such  that  quantum  confinement'  effects  do  not 
exist,"  is:  even  though  there  is  a  large  distribution  of  particle  sizes,  quantum  confinement  does  exist  in 
each  particle. 
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Table  1:  The  measured  energy  peak  of  first  derivative  of  reflection  at  300K,  bulk  energy  gap3  at  4.2k  and 
JOOK,  the  calculated  confinement  energies,  emission  energies  for  IS  *  IS  and  IP  -►  IP  transitions  and  effec¬ 
tive  diameters  (d),  the  measured  lifetimes  for  IS  -»  IS  and  IP  *  IP  transitions  as  well  as  their  ratio  are 
lisplayed.  The  t  are  anticipated  values  assuming  a  same  ratio  t13/t,p  for  the  samples  2-58  and  2-64. 
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Fig.  1  Time-resolved  photoluminescence  profiles  detected  at  different  energies  (a  to  e)  for  the  sample  2-58 
at  n.2K.  The  dotted  curves  are  the  theoretical  fit  with  12  ps  for  rlsetlme  and  decay  times  are  indicated  on 
the  corresponding  profiles.  The  left  curve  of  (a)  is  the  temporal  response  profile  to  the  laser  pulse  used 
to  excite  the  samples. 
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Fig.  2  The  lifetimes  are  plotted  as  function  of  emitted  photon  energy.  The  sample  temperature  for  (a),  (b), 
(c),  (d),  13  4.3K  and  for  (e)  is  2K.  The  arrows  indicated  the  energy  positions  for  1S->1S  transition  (S)  and 
1P-MP  transition  (P). 
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Fig.  3(a)  Photolumlnescence  spectra  of  sample  2-61  at  various  excitation  levels  taken  at  room  temperature 
using  «88  nm  line  of  the  Argon-Ion  laser.  The  excitation  levels  (Intensity  scale)  were  2.5  (2)  and  10  (10) 
times  higher  for  the  broken  curve  and  dot-dashed  curve  than  the  solid  curve,  respectively.  The  broken  verti¬ 
cal  line  Indicates  the  peak  position  of  the  first  derivative  of  reflectance  at  room  temperature.  S  and  P 
stand  for  the  1  S-*l S  and  1P-1P  transitions,  respectively.  (b)  Photolumlnescence  spectra  of  2-61  at  a. 3  k. 
The  Inset  shows  that  the  IS  peak  broadening  Is  dominated  by  the  diameter  fluctuation  of  mlcrocrystallltes. 
The  solid  dots  are  calculated  by  assuming  a  Gaussian  shape  of  the  probability  density  function  for  the  random 
variable  d  with  a  variance  of  25  A. 


158  /  SPtE  Vol.  793  Ultrafast  Las  ft  Probe  Phenomena  in  Bulk  and  Microstructure  Semiconductor  t  (19871 


531 


m 


'.i.rfPi 


PHYSICAL  REVIEW  B 


VOLUME  37,  NUMBER  3 


15  JANUARY  1988-1! 


Reply  to  “Comment  on  ‘Determination  of  valence-band  discontinuity  via  optical  transitions 

in  ultrathin  quantum  wells’  ” 

» 

■’  '  ’  Kai  Shum,  C.  Zhang,  P.  P.  Ho,  and  R.  R.  Alfano 

Institute  for  Ultrafast  Spectroscopy  and  Lasers,  Departments  of  Electrical  Engineering  and  Physics, 

The  City  College  of  Sew  York,  Sew  York,  Sew  York  10031 
(Received  16  January  1987;  revised  manuscript  received  20  August  19871 

The  salient  points  of  our  previous  paper  [Phys.  Rev.  B  33,  7259  (1986)]  are  reinforced  in  this 
Reply  to  the  preceding  Comment  by  Miller. 


The  essential  points  of  our  work1  were  to  demonstrate 
what  experimental  data  from  optical  transitions  should 
be  used  and  how  sensitive  they  are  to  the  Q  value  of  band 
offset.  Contrary  to  Miller’s  comment.2  Dingle  et  a! . J  did 
not  emphasize  how  sensitive  the  Q  value  depended  upon 
the  value  of  the  energy  separation,  A E(L.),  between  the 
n  =1  heavy-  and  light-hole  subbands.  We  found1  that 
the  values  of  A E  are  most  sensitive  to  the  band  offset  for 
L,  ranging  from  15  to  80  A.  Based  on  this  key  finding  a 
currently  fashionable  choice  of  Q,,  =0.40  proposed  by 
Miller  et  al.A  was  ruled  out  using  Dingle’s  connection 
rule  and  the  then  available  experimental  data  in  this 
range  of  L..  Another  key  point  made  in  our  paper1  was 
that  systematic  measurements  should  be  performed  in 
the  sensitive  zone  of  well  width  (15  A  to  80  A)  to  pre¬ 
cisely  determine  the  Q  value.  These  essential  points  are 
still  overlooked  in  depth  by  Miller's  comment2  to  our  pa¬ 
per.1 

It  should  be  pointed  out  here  that  the  extent  of  data 
used  in  our  paper1  was  not  the  main  issue  of  the  paper 
and  more  data  could  not  alter  the  essential  spirit  of  the 
paper.  Most  recently,  Miller  et  al.  have  performed  new 
measurements  providing  more  data  in  the  sensitive  zone 
which  will  be  discussed  below  in  this  Comment. 

We  have  recognized  that  there  has  been  a  large  body 
of  data  pointing  to  a  larger  Q„~0.40,  which  was  the 
mainstream  of  thought  in  this  field.  However,  most  of 
the  experimental  probes  carry  their  own  source  of  uncer¬ 
tainties.6  The  electrical  measurements  are  often  plagued 
by  residual  doping  whereas  the  intersubband  optical 
transitions  may  prove  to  be  too  weakly  dependent  upon 
the  band  offset.  The  optical  method  described  in  our  pa¬ 
per1  provides  a  sensitive  optical  test  of  band  offsets  in 
GaAs/Alj.GaAs,.,  structures.  It  should  be  mentioned 
that  there  are  some  most  recent  determinations  of  Q 
value  for  GaAs/AI^Ga^^As  which  do  not  yield 
Qu=0.40  but  point  to  values  of  0.31  from  Raman 
scattering  data7  and  0.23  ±0.07  from  electrolyte 
electroreflectance  study.’ 

As  noticed  by  Miller,2  one  of  the  data  in  the  original 
paper1  to  support  lower  Qt,  was  obtained  in  a  circular 
fashion.  However,  the  other  data  quoted  by  us  do  not 
suffer  from  this  problem. 

It  is  well  known  that  the  mismatch  of  effective  mass  in 
well  and  barrier  is  responsible  for  the  appearance  of 


various  connection  rules''9'10  for  the  envelope  wave  func¬ 
tion  and  its  derivative.  In  general,  using  a  different  con¬ 
nection  rule  will  result  in  different  eigen  energy  levels 
and  hence  differences  in  the  spacing  between  the  levels. 
Therefore,  the  Q  value  determined  by  A E  versus  L.  de¬ 
pends  much  on  what  connection  rule  is  employed.  This 
is  where  some  of  the  problems  come  from. 

In  the  spirit  of  our  past  paper,  the  light-  and  heavy- 
hole  energy  splitting  is  plotted  in  Fig.  1  as  function  of 
well  thickness  Lz  using  various  connection  rules. 3'9, 11 
The  Q  value  is  adjusted  in  such  a  wav  the  A£s  for 
Lz>  40  A  calculated12  using  different  connection  rules 
coincide  with  each  other  within  —2  meV.  D,  B ,  and  S 
in  the  Fig.  1  refer  to  the  Dingle’s  connection  rule, 
Bastard's  connection  rule,  and  new  connection  rule.11 
The  material  parameters  used  for  the  calculation  of  the 


Lz  (A) 


FIG.  1.  Calculated  light-  and  heavy-hole  energy  splitting  as 
function  of  L,.  D.  B,  and  .V  refer  to  Dingle’s.  Bastard's,  and 
new  connection  rules.  The  values  of  Q„  are  indicated  on  the 
corresponding  curves.  The  masses  used  for  the  dotted  curve  is 
Miller’s  masses.  The  masses  used  for  the  other  curves  are  the 
conventional  masses  given  in  the  text.  The  plus  and  triangle 
data  are  from  Refs.  5  and  15,  respectively. 
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curves  are  jc=0.37,  A£,.  =  1247x£)ll  (meV)  (the  value  of 
Q,,  is  indicated  next  to  the  corresponding  variable  in  the 
Fig.  1),  masses11,14  of  heavy  hole  (light  hole)  for  GaAs 
and  AlAs  are  0.403m,,  <0.087m(i)  and  0.487mfi 

(0.208mo),  respectively;  except  for  the  curve  where  the 
heavy-hole  mass  (0.34m0)  and  light-hole  mass  (0.094m„) 
proposed  by  Miller  et  al .*  were  used.  The  plus  and  tri¬ 
angle  data  points  in  Fig.  1  are  from  Refs.  5  and  15,  re¬ 
spectively. 

Figure  1  shows  two  important  features.  First,  as  ex¬ 
pected  the  energy  splitting  A £  is  sensitive  to  the  connec¬ 
tion  rule  in  the  sensitive  zone  of  Lz,  especially  for 
Lz  <40  A.  This  features  makes  it  possible  to  experimen¬ 
tally  demonstrate  which  connection  rule  is  appropriate 
by  systematically  measuring  A £  in  the  sensitive  zone. 
Based  on  the  data  given  by  Miller  et  al.,s,l>  Dingle’s 
connection  rule  which  was  used  by  us  to  show  how  to 
determine  Q  value  in  the  sensitive  zone  of  Lt  seems  to  be 
not  appropriate.  However,  the  method  itself  described 
in  our  paper1  is  still  very  useful.  Second,  for  the  given 
connection  rule  the  splitting  A £  is  very  sensitive  to  the 
effective  masses.  In  order  to  fit  simultaneously  the  data 
given  by  Miller  et  al. 5,15  it  seems  to  be  necessary  to  use 
the  Bastard’s  connection  rule  and  the  mass  parameter  set 
(mHn  =0.34/nn,  m LI|  =0.094/n0 )  and  £3,.  =0.40  pro¬ 
posed  by  Miller  et  al.*  The  use  of  new  connection  rule11 
with  Qu  =0.23  seems  to  fit  the  data  for  Lz>  30  A.  It 
should  be  pointed  out  here  that  in  our  calculation  the 
binding-energy  difference  between  the  light  hole  and  the 
heavy  hole  is  taken  to  be  a  constant  ( —  0.5  meV).  This 
may  be  a  good  approximation  for  Lz  >  40  A.  In  order  to 
comparethe  measured  A £  with  the  calculated  A £  for 
£,<40  A  an  exact  analysis  taking  the  proper  exciton 
binding  energy  into  account  must  be  performed.16 

In  Fig.  2,  the  energy  separation  A En  between  first 
and  second  conduction  subbands  is  plotted  as  function  of 
Lt  using  three  different  connection  rules  by  including17 
the  energy-dependent  mass  m„.(£„)  (Ref.  18)  as  well  as 
mblE„  ),  n  =  1,2,  for  well  material  GaAs  and  barrier  ma¬ 
terial  Al0 jGa0  7As,  respectively,  and  compared  A£,: 
with  the  published  experimental  data  by  West  and 
Eglash.l<>  This  work  is  not  affected  by  the  valence-band 
complexity,  the  exciton  binding  energy,  and  can  provide 
an  independent  test  of  Q  value. 

The  mb(E„)  is  obtained  by 


F 2 

fi2k  2£< 

1/2  £ 

d2 

ft2k2 

dk1 

4  + 

2m  cd 

2 

dk2 

2mb(En ) 

(I) 

The  expression  inside  of  the  large  brackets  on  the  left- 
hand  side  of  Eq.  (1)  is  the  energy-momentum  relation¬ 
ship  within  semiconductor  band  gap.20  After 
differentiation  the  mb(E„)  reduced  to 

mb(E„  )=mCB(  1  +4V„  ),/J,  n  =  l,2, ...  ,  (2) 

where  mCB  is  electron  effective  mass  at  conduction  band 
edge  of  AI0  JGa07As,  =  |  ff2fc„2„  /2mCB£t  | ,  and 
n  =  l,2  for  first  and  second  subbands.  Note  that  k2  is 


Lr  (A) 


FIG.  2.  Calculated  energy  separation  A£i;  between  first  and 
second  conduction  subbands  as  a  function  of  well  width  L.  us¬ 
ing  Dingle's  (Z>),  Bastard's  (B),  and  new  (,V)  connection  rules. 
The  curves  are  calculated  using  A£r=(?fA£,  (A£,  =  1247* 
meV,  jc  =  0. 30) .  The  values  of  Q,  arc  indicated  on  the  corre¬ 
sponding  curves.  The  data  of  West  and  Eglash  I  Ref.  |9>  are  in¬ 
dicated  as  two  plus  signs.  The  inset  shows  schematically  the 
well  and  the  energy-momentum  dispersion  curves.  The  elec¬ 
tron  effective  masses  at  band  edges  for  GaAs  and  Al0  ,Ga„  ,As 
are  also  indicated  in  the  inset. 

negative  and  kb  is  imaginary  reflecting  the  electron  wave 
function  within  the  gap  of  Al03Ga07As  is  exponentially 
damped.  The  dispersion  relations  dependent  upon  the 
real  and  imaginary  wave  vectors  in  the  well  and  the  bar¬ 
rier  are  schematically  shown  in  the  inset  of  Fig.  2.  The 
connection  rule  and  the  Q  value  used  are  indicated  on 
the  corresponding  curves.  The  value  of  .t  is  equal  to 
0.30  in  order  to  make  a  comparison  with  the  data  A£,2 
of  West  and  Eglash.19 

Several  features  are  very  apparent  in  Fig  2.  First,  for 
the  given  connection  rule  A£l2  is  very  sensitive  to  the  Q 
value  in  the  sensitive  zone  of  L.  giving  rise  to  a  possibili¬ 
ty  of  extracting  the  value  of  A£c  and  Q  value.  Second, 
the  value  of  A£l2  calculated  using  the  new  connection 
rule  approaches  the  values  of  A£]2  at  —45  A  and  at 
Lz  >  100  A  calculated  using  Bastard's  and  Dingle's  con¬ 
nection  rules,  respectively.  Third,  the  two  experimental 
points  support  neither  Bastard's  nor  Dingle’s  connection 
rules,  but  seems  to  be  more  favorable  to  the  new  connec¬ 
tion  rule  with  Qc  =0.77. 

In  conclusion,  we  have  shown  different  connection 
rules  and  hole  masses  yield  different  Q  values.  Further¬ 
more,  the  appropriate  connection  rule  can  be  discrim¬ 
inated  by  systematically  measuring  A  £  and  A£,2  in  the 
sensitive  zone  of  Lz.  Dingle’s  85-15  rule  is  based  on  his 
connection  rule  while  Miller’s  b0-40  rule  is  based  on 
Miller's  effective  masses  and  Bastard's  connection  rule. 
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The  essential' information  described  in  our  paper1  is  not 
compromised  by  the  comments  raised  by  Miller  but  rein¬ 
forced  here.  We  stress,  that  more  systematic  photo- 
luminescencer  PLE,  and  inter-conduction-subband  tran¬ 
sition  data  must  be  taken  in  the  sensitive  zone  (15  to  80 


A)  to  determine  the  proper  connection  rule  and  an  accu¬ 
rate  Q  value. 

This  work  was  supported  by  Air  Force  of  Scientific 
Research  under  Grant  No.  AFOSR-86-003 1 . 
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Gallium  arsenide  photoluminescence  under  picosecond-laser-driven  shock 
compression 
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A  pump-and-probe  technique  was  used  to  investigate  shock  effects  on  the  photoluminescence 
spectra  (  ~  833  nm)  at  T  =  80  K  due  to  the  direct  transition  £0  from  the  T6  conduction  band 
to  the  r„  fourfold  degenerate  top  valence  band  in  GaAs.  Under  the  shock  loading  condition, 
the  photoluminescence  peak  was  observed  to  blue  shift  and  split  into  two  components, 
corresponding  to  the  transitions  from  the  r6  conduction  band  to  the  valence  heavy-  and  light- 
hole  subbands,  because  of  symmetry  breaking  by  the  uniaxial  shock  compression  along  the 
{001  ]  direction.  From  the  blue  shift  of  the  photoluminescence  peaks,  we  deduced  our 
picosecond-laser-driven  shock  pressure  of  ~  10  kbar. 


Significant  changes  in  the  optical  and  transport  proper¬ 
ties  of  semiconductors  occur  with  the  application  of  a  uniax¬ 
ial  stress  due  to  changes  in  symmetry  and  lattice  param¬ 
eters. 1-1  Knowledge  of  stress  effects  on  band  structure  is 
important  for  operation  of  GaAs  in  high-speed  ultrasmall 
devices  and  switches.  In  this  letter  we  report  on  a  new  obser¬ 
vation  of  the  photoluminescence  spectra  from  gallium  arsen¬ 
ide  ( GaAs)  under  picosecond-laser-driven  shock  wave  load¬ 
ing.  The  shock  waves  were  generated  by  focusing  intense 
picosecond  laser  pulses  onto  an  aluminum  (Al)  foil5'"  at¬ 
tached  to  the  sample. 

GaAs  has  a  zinc-blende  structure9  of  the  space  group 


T l .  The  lowest  conduction  band  is  s-like  and  twofold  degen¬ 
erate,  having  a  T6  symmetry  at  the  Brillouin  zone  center. 
The  top  valence  band  is  a  sixfold  degenerate />-like  band  that 
splits  into  a  fourfold  Pir_  multiplet  ( J  =  3/2,  Mj  —  ±  3/2, 
+  1/2)  with  a  r„  symmetry  and  Ptn  doublet  ( J  =  1/2, 
Mj  =  +1/2)  with  a  T7  symmetry  due  to  the  spin-orbit  in¬ 
teraction.  The  band  structure  is  shown  in  Fig.  1  (the  left- 
hand  side).  The  lowest  energy  band-to-band  direct  transi¬ 
tion  (labeled  £0)  occurs  at  the  center  of  the  Brillouin  zone. 
The  total  Hamiltonian  matrix  for  the  valence  band  in  the 
presence  of  a  uniaxial  compression  P  along  the  {001 )  direc¬ 
tion  can  be  written9 


-6Eh-\8Es 


SEh  +  J  8ES 
2~'n8E*. 


—  A0  —  8EH 


where  8E„  =  —a(Su+2Su)P,  8EH  = 

+  2 SU)P,  SEs  =  -2b(Su -Sl2)P,  8ES 

=  —  2h'(Sn  — S^J/’.andA,,  is  the  spin-orbit  splitting,  Sn 
and  5, 2  are  components  of  the  elastic  compliance  constants, 
a  and  a'  are  the  deformation  potentials  of  the  Tg  and  T7 
valence  bands  for  hydrostatic  effects,  and  b  and  b '  are  the 
deformation  potentials  of  the  T8  and  T,  bands  for  shear  ef¬ 
fects.  If  the  pressure-induced  shift  and  splitting  are  much 
smaller  than  A„,  the  eigenvalues  of  Eq.  ( 1 )  can  be  written 


<5£,.  =  -8Eh-  1 5ES, 
8E,.  =  —  6Eh  +  ( SEs  + 


(8ES)- 


f>E,.  =  -  A0  -  8Eh  -  (<5£s—  .  (2c) 

2A0 

Including  the  effect  of  the  conduction-band  shift,  the  pres¬ 
sure  dependence  of  the  energy  gap  changes  is  given  by 

8(Et  -£„.)  =  A£„  +J  6ES,  (3a) 


8 (Ec  -£„,)  =  A£„  -  J  8ES  - 


(8ES  )2 


8(EC  —  £,,t )  =  A0  +  A£w  + 


2Ao 
(8ES  ) 2 


where  A£H  =  —  (c,  +  j)(5,,  +  2£u)£.  A EH 

=  —  (c,  +a')  (S,i  +  2SU)P,  andc,  is  the  deformation  po¬ 
tential  of  the  conduction  band.  For  light  polarized  perpen¬ 
dicularly  to  the  stress  axis  c,  the  transitions  between  the  s- 
like  T6  conduction  band  and  all  three  valence  bands,  v,,v2, 
and  vy,  are  allowed.9  From  Eqs.  (3a)  and  (3b),  we  deduce 
that  under  shock  loading  condition,  the  photoluminescence 
peak  corresponding  to  the  £„  transition  will  be  subject  to  a 
blue  shift  because  of  the  hydrostatic  component  of  the  pres¬ 
sure  [represented  by  the  term  A EH  in  Eqs.  (3a)  and  (3b) ). 
It  will  also  split  into  two  components  corresponding  to  the 
transitions  £„(  l )  and£0(2  ;  (see  Fig.  1,  the  right-hand  side) 
because  of  the  reduced  symmetry  by  the  uniaxial  nature  of 
the  shock  pressure  [  represented  by  the  term  8ES  in  Eqs.  ( 3a) 
and  (3b)  J.  Since  the  heavy-hole  band  in  has  a  lower  energy 
under  uniaxial  compression,10  the  corresponding  transition 
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UNSHOCKED  SHOCKED 

FIG.  1.  Band  structure  of  GaAs  at  T  =  80  K  for  unshocked  (left)  and 
shocked  (right)  crystals. 


£0(2)  will  have  a  larger  blue  shift.  In  addition,  the  blue  shift 
of  the  emission  peak  of  the  electron-light-hole  recombina¬ 
tion  l  the  transition  £„(  1 )  1  will  increase  sublinearly  with  the 
increasing  pressure  because  of  the  pressure-induced  cou¬ 
pling  between  the  sets  of  Mj  ~  +  1  bands  (v,  and  v,). 

A  schematic  diagram  of  the  experimental  geometry  is 
shown  in  Fig.  2.  The  experimental  setup  was  described  in 
Ref.  11.  A  Quantel  NdrYAG  laser  of  30  ps  pulse  width  and 
25  mJ  pulse  energy  at  l  .064  (im  was  utilized  as  the  pump 
beam  to  generate  shock  waves  in  the  aluminum  foil.  The 
shock  waves  were  launched  into  the  GaAs  sample  located 
behind  the  aluminum  foil.  The  probe  beam  of  27  ps  at  532 
nm  was  delayed  through  a  white  cell  approximately  by  46  ns 
relative  to  the  arrival  of  the  pump  beam  at  the  aluminum  foil 
to  ensure  that  the  probe  pulses  reach  the  sample  surface  area 
right  after  the  arrival  of  the  shock  front.  This  was  necessary 
to  observe  the  maximum  effect  of  shock  pressure. 

The  GaAs  sample  was  n  type  with  a  carrier  concentra¬ 
tion  of  ~  10"1  cm~\  The  thickness  of  the  sample  was  about 
50 /tm  with  the  cubic  [001 J  axis  perpendicular  to  the  sample 
surface.  The  aluminum  foil  had  a  thickness  of  20  /im.  The 
sample  and  aluminum  foil  were  loaded  between  two  quartz 
plates  and  placed  in  the  center  of  a  liquid-nitrogen  Dewar. 


APIEZON  N-grease  was  added  to  the  interface  of  the  sam¬ 
ple  and  aluminum  foil  to  ensure  efficient  mechanical  cou¬ 
pling  at  T  —  80  K.  The  pump  beam  was  focused  by  a  1 5-cm 
focal  length  lens  to  a  500-/im-diam  spot  on  the  aluminum 
foil,  which  ensured  planarity  of  the  shock  wave  and  a  peak 
power  density  of  —  4.2X  10"  W/cm:.  The  probe  pulse  was 
focused  to  a  400-/im-diam  spot  by  a  15-cm  focal  length  lens 
on  the  sample,  giving  a  peak  power  density  of  —  1.6  X  107  W/ 
cm:. 

The  emitted  light  was  dispersed  by  a  1/4  m  spectrom¬ 
eter,  detected  by  a  silicon-intensified  target  coupled  to  an 
optical  multichannel  analyzer  and  stored  in  a  PDP  11/23  + 
host  computer.  The  resolution  of  the  detection  system  was 
about  2  nm.  All  spectra  were  recorded  at  T  —  80  K  by  single 
shots  of  the  laser. 

The  observed  photoluminescence  spectrum  of  GaAs 
without  shock  loading  is  dominated  by  a  single  peak  due  to 
the  band-to-band  recombination  at  the  lowest  energy  gap  £,„ 
as  shown  in  Fig.  3  ( left ) .  No  emission  lines  due  to  impurity 
transitions  were  observed.  The  effect  of  shock  waves  was 
most  pronounced  for  the  delay  time  between  the  Dump  and 
probe  beam,  r  —  40  —  60  ns.  At  r<  35  ns,  no  shock  effect  was 
observed.  A  typical  spectrum  showing  the  shock  effects  on 
the  photoluminescence  from  GaAs  is  displayed  in  Fig.  3  for 
r  =  46  ns  and  probe  intensity  1.6  X  10  W/cm\  Under  shock 
loading,  the  emission  peak  was  split  into  two,  with  blue  shifts 
of  20  and  85  meV.  respectively. 

The  lowest  energy  band  gap  in  GaAs  in  the  absence  of 
the  shock  pressure  can  be  written 

E,lT)  =  1.530eV  -  (50x  10  ‘eVK  ')T.  (4) 

The  peak  energy  of  the  emission  due  to  the  band-to-band 
transition  £„  at  T  =  80  K  is  Er  —  £„  ( 80  K )  =  1.4W  eV, 
i.e.,  A.,,  =  832  nm,  which  agrees  with  our  observed  un¬ 
shocked  emission  line  at  833  nm  (  Fig.  3 )  within  our  experi¬ 
mental  accuracy.  Since  the  energy  gap  E,  —  £,  expands 
linearly  with  the  increasing  pressure  P,  the  blue  shift  of  the 
transition  energy  £„( 2)  was  chosen  to  calibrate  the  shock 
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FIG  .V  Unshocked  (left)  and  shocked  (fight)  photoluminescence  spectra 
al  adelav  lime  r  =  46nsand  probe  iniensm  I  h  .  IO'W/cm:  Theshocked 
spectrum  consists  of  two  peaks,  corresponding  lo  ihe  transitions  £,,(  1 )  and 
£„(2>. 
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pressure.  From  Eq.  (3a)  the  pressure  coefficient  of  this  tran¬ 
sition  energy  is  given  by 

dE.,(  2) 

=  -  (c,  +aHS,t  +2Si2) -b(Sn  —5,,). 

(5) 

With  the  reported  values,13  c,+a=  —13.0  eV. 
6=  -1.66  eV,  and14  S,,  =  1.16X  10“ 3  kbar'\ 

=  —  3.67 X  10-4  kbar-',  we  obtain  dEt)(2)/dP  =  8.5 
meV/kbar.  Therefore,  the  experimentally  measured  blue 
shift  of  85  meV  (Fig.  3)  for  the  electron-heavy-hole  recom¬ 
bination  emission  line  gives  a  magnitude  of  the  laser-driven 
shock  pressure  in  GaAs  of  ~  10  kbar. 

From  Eq.  (3b)  and  the  reported  values,4  b '  =  —  2.47 
eV  and  A()  =  0.34  eV,  we  obtain  S(EC  -  E,: )  =  21.5  meV 
for  a  pressure  of  10  kbar,  which  consistently  agrees  with  our 
observed  blue  shift  of  ~  20  meV  under  the  same  shock  com¬ 
pression  condition.  Therefore,  the  shocked  system  is  in  qua¬ 
siequilibrium  state. 

In  conclusion,  the  picosecond-laser-drive  shock  waves 
offer  a  new  technique  to  investigate  properties  of  materials  at 
very  high  pressures  with  a  table-top  laser  facility.  Using  this 
technique,  we  have  observed  blue  shift  and  line  splitting  in 
the  photoluminescence  spectrum  of  GaAs  due  to  band-gap 
expansion  and  symmetry  breaking  by  the  uniaxial  shock 
compression  along  the  [001  ]  direction.  The  pressure  range 
of  the  shock  waves  generated  in  our  experiment  by  a  25-mJ 
picosecond  laser  pulse  is  —  10  kbar. 
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Shock-wave-induced  collision  broadening  of  the  photoluminescence  spectra 
in  GaSe 
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Significant  spectral  broadening  of  the  photoluminescence  in  GaSe  under  the  picosecond-laser- 
driven  shock  pressure  has  been  observed  for  the  first  time.  The  broadening  of  the  spontaneous 
emission  was  found  to  be  proportional  to  the  shock  pressure  and  attributed  to  a  shock-wave- 
induced  exciton  collision  mechanism  due  to  the  directional  motion  of  particles  in  the  shocked 
region. 


There  are  significant  differences'"*  between  the  ob¬ 
served  physical  and  chemical  properties  of  condensed  mat¬ 
ter  which  is  shocked  and  hydrostatically  compressed  at  the 
same  pressure  and  temperature.  Production  of  the  picose- 
cond-laser-driven  shock  waves  ofTers  a  new  powerful  tech¬ 
nique  to  study  the  properties  of  semiconductor  circuits  and 
devices  at  very  high  pressures  with  a  table-top  laser  facility. 
Recently,  we  reported3  on  the  spontaneous  and  stimulated 
emission  from  GaSe  at  different  excitation  intensities  and 
attributed  both  the  spontaneous  and  stimulated  emission  to 
the  same  origin — exciton-exciton  scattering  process.  In  this 
letter,  photoluminescence  spectral  line  broadening  observed 
under  picosecond-laser-driven  shock  loading  in  GaSe  is  re¬ 
ported  for  the  first  time  which  is  attributed  to  enhanced  exci- 
tonic  collision  processes  associated  with  the  shock  wave. 

A  pump-and-probe  technique  was  used  to  observe  the 
shock  effect  on  the  photoluminescence  emission  in  GaSe. 
The  GaSe  sample  of  50/rm  thickness  with  the  c  axis  perpen¬ 
dicular  to  the  layers  was  attached  to  an  aluminum  foil  of 
thickness  20  fim.  A  Quantel  Nd:YAG  laser  of  30  ps  pulse 
width  and  35  mj  pulse  energy  at  1.064 /im  was  focused  to  a 
450-/ii.i-diam  spot  on  the  A1  foil  surface  to  generate  shock 
waves  which  propagated  through  the  A1  foil  into  the  GaSe 
sample.  A  532-nm  probe  beam  of  27  ps  pulse  width  was  fo¬ 
cused  to  a  350-/im-diam  spot  onto  the  GaSe  surface.  The 
photoluminescence  from  GaSe  was  dispersed  by  a  1/4-m 
spectrometer  and  detected  by  a  silicon-intensified  target 
coupled  to  an  optical  multichannel  analyzer  OMA  III,  and 
stored  in  a  PDP1 1/23  -|-  computer. 

Typical  spontaneous  emission  spectra  from  GaSe  with 
and  without  shock  loading  are  shown  in  Fig.  1 .  The  salient 
features  displayed  in  Fig.  1  are  red  shift  of  28  nm  and  broad¬ 
ening  of  the  shocked  spectrum  relative  to  the  unshocked  one. 

The  emission  due  to  exciton-exciton  scattering  can  be 
expressed  as  (£*,£*  )  —  (hv,e  —  h).  The  spontaneous 
emission  spectrum  due  to  the  exciton-exciton  scattering  pro¬ 
cess  is  given  by6 


(£  -  E,  )2  +  (va/e )£2  kT , 


« _ & _ 


wher  ep(E)dE  is  the  number  of  photon  modes  in  the  crystal 
between  energies  E  and  £  +  dE,  wherep(£)  <x  £ 2,  E,  is  the 
energy  of  the  Is  exciton  at  K  =  0,  Eb  is  the  binding  energy  of 
excitons,  Ex  =  Et  —  Eb,Tx  is  exciton  temperature,  a  and  e 
are  the  polarizability  of  excitons  and  dielectric  constant  of 
the  material,  respectively.  The  peak  energy  is  located  at 
Ep  —  Eg  —  2 Eb .  Figure  2  shows  a  comparison  between  the 
observed  spontaneous  emission  spectrum  (solid  dots)  and 
the  calculated  one  (solid  line)  from  Eq.  ( 1 )  with  Ex  =  2.0 
eV,  £„  =  0.020  eV,  Tx  =  500  K,  and  the  coupling  coeffi¬ 
cient  ira/e  =  3.3  X  10"3. 

It  is  a  well  known  fact  that  the  direct  energy  band  gap  of 
GaSe  shrinks  under  pressure  causing  the  red  shift  of  the  pho¬ 
toconductivity,  absorption,  and  luminescence  spectra.7'’ 
With  the  reported  static  values  of10  dEt/dP  =  —6.2 

meV(kbar) " 1  and  3Eb/3P  =  -  0.6  meV(kbar)  ” ’  we  ob¬ 
tain  the  ( hydrostatic )  pressure  coefficient  of  the  peak  energy 
3Ep/3P  =  —  5.0  meV(kbar)  ~ '.  There  is  evidence1 1  show¬ 
ing  that  the  stress  anisotropy  of  the  shock  compression  could 
be  neglected  in  many  cases.  The  red  shift  of  28  nm  ( 84  meV ) 
displayed  in  Fig.  1  deduces  a  shock  pressure  of  ~  1 7  kbar.  To 
see  if  this  approximation  is  reasonable,  we  also  performed  an 
independent  measurement  of  the  shock  pressure  in  the  same 
experimental  setup  by  an  jc-cut  quartz  transducer  tech¬ 
nique. 12  The  pressure  of  —  1 5  kbar  was  obtained  which  con¬ 
firmed  our  approximation. 

Significant  line  broadening  of  the  spontaneous  emission 
spectra  under  the  shock  loading  was  clearly  and  repeatedly 
observed.  The  magnitude  of  broadening  was  about  24  nm  for 
a  shock  pressure  of  17  kbar,  as  displayed  in  Fig.  1.  It  was 


WAVELENGTH  ( nm ) 

FIG.  I.  Observed  unshocked  and  shocked  spontaneous  emission  spectra 
from  GaSe  at  room  temperature  and  excitation  intensity  17  MW/cm2. 
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FIG.  2.  Calculated  spontaneous  emission  spectrum  from  the  exciton-exci- 
ton  scattering  model  of  Eq.  ( 1 )  with  E,  =  2.0  eV,  Eb  =  0  020  eV,  not/ 
f  =  3.3  X  10  and  T ,  =  500  K.  The  solid  dots  represent  the  experimental 
unshocked  spontaneous  spectrum. 

found  that  the  peak  amplitude  decreased  under  shock  load¬ 
ing  such  that  the  integrated  spontaneous  emission  intensity 
(area)  approximately  remains  the  same  as  that  of  the  un¬ 
shocked  spectrum. 

It  was  reported13  that  although  the  electronic  (or  exci- 
tonic )  temperature  in  the  probed  region  ( ~  434  K )  was  con¬ 
siderably  higher  than  the  room  temperature  due  to  heating 
effect  by  the  probe  beam,  no  significant  temperature  rise  due 
to  the  shock  wave  loading  was  detected.  The  Doppler 
linewidth  at  temperature  Tx  due  to  the  thermal  motion  of 
the  emitting  excitons  is  given  by14 


where  kp  is  the  peak  wavelength  of  the  emission  and  Mx  is 
the  exciton  mass.  If  this  mechanism  was  responsible  for  the 
observed  line  broadening  (of  a  factor  of  2,  approximately), 
the  temperature  would  increase  by  a  factor  of  4  due  to  the 
shock  loading,  i.e.,  the  exitonic  temperature  would  exceed 
1700  K  in  the  shocked  region  which  is  definitely  unrealistic. 
Also,  a  Doppler  broadened  line  should  have  a  Gaussian  pro¬ 
file14  which  disagrees  with  our  observed  Lorentzian  line 
shape.  Therefore,  the  possibility  of  temperature  effect  caus¬ 
ing  the  linewidth  increase  due  to  the  shock  loading  can  be 
safely  excluded. 

The  probed  region  in  our  experiment  formed  a  disk  of 
350 //m  in  diameter  and  5^im  in  thickness  (the  penetration 
depth  of  the  532-nm  probe  beam  in  GaSe).  Since  the  shock 
pressure  is  proportional  to  the  pump  beam  fluence,12  trans¬ 
verse  inhomogeneity  of  the  shock  pressure  over  the  disk  area 
(350/im  in  diameter)  originated  from  the  Gaussian  profile 
of  the  pump  beam  (450/rm  in  diameter)  intensity  should  be 
less  than  20%.  There  was  also  a  longitudinal  inhomogeneity 
of  the  shock  pressure  in  the  direction  of  shock  wave  propaga¬ 
tion  due  to  pressure  decay  in  space.  In  our  picosecond  laser- 
driven  shock  wave  case,  the  shock  decay  time  was  mea¬ 
sured13  to  be  —70  ns,  corresponding  to  a  decay  length  of 
—  140/rm.  Thus,  the  longitudinal  variation  of  pressure  in  the 
shocked  region  was  less  than  0.5  kbar.  However,  the  ob¬ 
served  line  broadening  corresponded  to  a  pressure  variation 
from  zero  to  twice  the  mean  value  of  the  shock  pressure  as 
deduced  from  the  observed  red  shift.  Therefore,  neither 
transverse  nor  longitudinal  inhomogeneity  of  the  shock 
pressure  can  explain  the  observed  large  line  broadening. 
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The  observed  line  broadening  of  the  spontaneous  emis¬ 
sion  is  attributed  to  the  shock-wave-induced  exciton  colli¬ 
sions.  When  a  shock  front  propagated  through  the  sample, 
all  molecules  behind  it  would  gain  a  particle  velocity  up  via 
the  intermolecular  bonding,  but  excitons  would  not,  since 
there  was  no  such  tight  bonding  between  excitons  or  between 
exciton  and  molecules.  Therefore,  the  emitting  excitons  in 
the  shocked  region  all  suffer  from  additional  collisions  with 
the  entire  array  with  a  directional  particle  velocity  up .  The 
collision  frequency  for  each  exciton  is  fc  — a:up  Nmo[,  where 
a  is  the  Bohr  radius  of  the  exciton  in  GaSe  which  is 1 5  —  32  A, 
Nmol  is  the  molecular  density  in  GaSe  which  is  ~2x  10” 
cm-3.  The  particle  velocity  at  17  kbar  is  about  2X104 
cms-1.  Therefore,  the  collision  broadening  at  17  kbar  is 
Av  =  2fc/ir~ 2.5x  1013  Hz  or  A/l  —  33  nm  which  agrees 
well  with  the  observed  value.  It  should  be  noted  that  colli- 
sional  broadening  is  homogeneous  and  should  have  a  Lor¬ 
entzian  line  shape. 14  This  is  also  clearly  shown  in  Fig.  1. 
Since  the  collision  frequency  is  proportional  to  up  which  is 
related  to  the  shock  pressure  P  by  the  jump  condition": 
P  —  P0  =  ppU,  up ,  where  Us  is  the  shock  velocity,  p0  and  P0 
are  density  and  pressure  in  unshocked  region.  The  shock¬ 
wave-induced  collision  broadening  is  expected  to  be  propor¬ 
tional  to  the  shock  pressure.  This  is  confirmed  by  our  experi¬ 
mental  data  as  shown  in  Fig.  3.  The  observed  line  broadening 
from  different  pairs  of  unshocked-shocked  data  plotted 
against  the  applied  shock  pressure  displays  a  linear  depen¬ 
dence.  A  pressure  coefficient  of  line  broadening  for  the  spon¬ 
taneous  emission  of  GaSe  —  1.3  nm/kbar  is  obtained. 

The  observation  that  the  shocked  spontaneous  emission 
has  the  same  integrated  intensity  as  the  unshocked  one  sug¬ 
gests  that  the  elastic  collision  (phase-perturbing  collision) 
plays  a  major  role  in  our  line  broadening  mechanism  because 
inelastic  collision  (quenching  collision)  will  decrease  the 
number  of  emitting  excitons  and  reduce  the  integrated  inten¬ 
sity  significantly.  Note  that  the  usual  pressure  broadening 
which  is  crucial  for  gaseous  systems  essentially  comes  from 
an  increase  in  particle  density  and  a  decrease  in  mean  free 
path  due  to  high  pressure.  This  latter  effect  cannot  be  signifi¬ 
cant  in  condensed  matter  because  of  its  extremely  small 


FIG.  3  Observed  wavelength  broadening  of  the  spontaneous  emission  spec¬ 
trum  from  GaSe  vs  shock  pressure. 
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compressibility  (typically  p/p„~  1.2  at  P  =  200  kbar).  The 
shocked-wave-induced  exciton  collision  broadening  is  a  new 
mechanism  different  from  the  usual  pressure  broadening 
and  directly  related  to  the  unique  nature  of  the  shock  waves 
and  excitons  in  condensed  matter. 
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